//stacksize('max')

x(1)=8030      
x(2)=8040       
x(3)=7998     
x(4)=8008   
x(5)=8037       
x(6)=7956      
x(7)=7937       
x(8)=7963     
x(9)=7919     
x(10)=7920   


y(1)=9478       
y(2)=9533       
y(3)=9578      
y(4)=9590     
y(5)=9597      
y(6)=9626      
y(7)=9627      
y(8)=9723     
y(9)=9766       
y(10)=9768   



z(1)=-2.39965177    
z(2)=1.81709373    
z(3)=-1.14068067    
z(4)=-2.95249772    
z(5)=-3.25881290    
z(6)=1.36802530    
z(7)=-1.46287727    
z(8)=-6.40911865    
z(9)=1.03929353    
z(10)=1.49115121   

n=length(z)

for i=1:n
rx(i)=round(x(i))
ry(i)=round(y(i))
end

write(%io(2),'length(rx)')
write(%io(2),length(rx))

x1=min(rx);
x2=max(rx);

ilen=round(x2-x1+1)
write(%io(2),'ilen')
write(%io(2),ilen,'(i5)')


y1=min(ry);
y2=max(ry);


jlen=round(y2-y1+1)

write(%io(2),'jlen')

write(%io(2),jlen,'(i5)')



for i=1:ilen
    for j=1:jlen
        newmat(i,j)=0.;
    end
end

sn=size(newmat)

write(%io(2),'size of newmat')
write(%io(2),sn)

for k=1:n
    i=rx(k)-x1+1;
    j=ry(k)-y1+1;
    newmat(i,j)=z(k);
end


snewmat=size(newmat);
write(%io(2),'size of newmat again')
write(%io(2),snewmat)


for i=1:ilen
    xaxis(i)=i;
end


for j=1:jlen
    yaxis(j)=j;
end



fig5=scf(5)
dmin=min(newmat)
dmax=max(newmat)

admin=abs(dmin)
admax=abs(dmax)

amax=max(admin,admax)
dmin=-amax
dmax=amax

//write(%io(2),dmin)
//write(%io(2),dmax)

xtitle("data")

xset("colormap",jetcolormap(64))
colorbar(dmin,dmax)  

grayplot(xaxis,yaxis,newmat)

h = gcf();
h.color_map = jetcolormap(64);

h.color_map(32,1)=1
h.color_map(32,2)=1
h.color_map(32,3)=1

h.color_map(33,1)=1
h.color_map(33,2)=1
h.color_map(33,3)=1