home *** CD-ROM | disk | FTP | other *** search
- for i=1:20,for j=1:20, x(i,j)=0;
- for i=1:20 , x(i,1)=12.55-i*1.26;
- for i=2:20, for j=1:20,x(i,j)=x(i,1);
- for i=1:20, for j=1:20, y(i,j)=0;
- for i=1:20, y(1,i)=9.7-i*.875;
- for i=1:20, for j=2:20, y(i,j)=y(1,j);
- for i=1:20, for j=1:20, if x(i,j)<>0, x(i,j)=sin(x(i,j))/x(i,j);
- for i=1:20, for j=1:20, if x(i,j)=0, x(i,j)=1;
- for i=1:20, for j=1:20, if y(i,j)<>0, y(i,j)=sin(y(i,j))/y(i,j);
- for i=1:20, for j=1:20, if y(i,j)=0, y(i,j)=1;
- for i=1:20, for j=1:20, z(i,j)=x(i,j)*y(i,j);
- plot(z,'sinc2.plt')
-