Sunset in the savanna

on 22 Oct 2021

X=-10:.02:10;

s=600./abs(X+i*X'-3i);

colormap([pink.*autumn;1 1 .9]);

image(s);

hold

Current plot held

E=39*F(0:1);

I=800-E/i;

fill(E+510,I,0);

Warning: Imaginary parts of complex X and/or Y arguments ignored.

h=scatter(E+510,I,50,-14*abs(E));

Warning: Using only the real component of complex data.

alpha(h,.1);

camva(3)

function V=F(p)

V=[];

if abs(diff(p))>.01

Z=p+.5i*diff(p)*[4;

4-i+i^cos(431*norm(p))^.5];

V=[p(1);

F(Z(1:2));

F(Z(2:3));

p(2)];

end