clear; h = 0:0.05:1; [x y] = meshgrid(h); z = fun1(x,y); mesh(x,y,z);
%print -deps draw_ex.eps
