- /
-
i
on 14 Oct 2021
- 19
- 50
- 2
- 0
- 264
X=-164:164;
[t,r]=cart2pol(X',X);
t=t+pi;
for k=0:.1:2*pi
t(t<k)=t(t<k)+rand;
end
t=rescale(t,20,50);%controls colour - indices into colormap
t(r>80)=258;
t(r<25)=0; %controls pupil size
colormap([[0 0 0];turbo;[1 1 1]]);
image(t)
hold on
k=[166 84 1;166 244 -1];
for j=1:2
plot(X+k(j,1),k(j,3)*.003*X.^2+k(j,2),'-k','LineWidth',4)
end
axis equal off