>> [x y]=meshgrid(0:100,0:60);

>> z=1635711.892-351.963775*x-2259.1306*y+42.74653*x.*y;

>> vpa(z,10);

zz=vpa(z,10);

>> fix(zz/100000);

>> (fix(zz/100)*100- fix(zz/100000)*100000)/100;

>> vpa(zz-fix(zz/100)*100,5);

>> zz1=fix(zz/100000);

>> zz2=(fix(zz/100)*100- fix(zz/100000)*100000)/100;

>> zz3=vpa(zz-fix(zz/100)*100,5);

plot3(zz)

hold on

plot3(zz)

hold on

plot3(zz1)

hold on

plot3(zz2)

hold on

plot3(zz3)

zzz(61,100)

a=1

b=1

>> for i=1:1:61

for j=1:1:101

if zz2(i,j)>150&10

zzz(a,b)=z(i,j)

a=a+1;

b=b+1;

end

end

end