function z = myfun(x,y,c)
%function z = myfun(x,y,c)
%evaluates f(x) = cx^2*y
%x and y are coordinates from meshgrid
%PC MA302
r=sqrt(x.^2+y.^2);
z = c*sin(r)./r;
end