
h=1/100;

f=inline('3*x+y','x','y');

% DEFINIAMO IL QUADRATO [a,b] x [c,d].
a=0; b=1;
c=0; d=1;

x=0:h:1;
y=0:h:1;

[X,Y]=meshgrid(x,y);
Z=feval(f,X,Y);

surf(X,Y,Z);

