% ex33.m
n=1:1000;
r=sqrt(n);
theta=137.51*pi*n/180;
x=r.*cos(theta);
y=r.*sin(theta);
plot(x,y,'o');
axis('equal');