function z = myfun(x,y,c)
%evaluates f(x) = cx^2*y
%x and y are coordinates from meshgrid
%PC 02/03 MA302

z = c*(x.^2).*y;
end