vectorize a four parameter function
see vectorize3 for lower parameter case