
clear a;

A1 = 20.;
B1 = 2.5;
A2 = 30.;
B2 = 3.5;

SX = 25;
SY = 25;

DX = (A2-A1)/SX;
DY = (B2-B1)/SY;

for i=1:SX;
 for j = 1:SY;
  a(i,j) = func(A1+DX*i,B1+DY*j);
end;
end;

clf;
mesh (a);
