Lunar Shadows

on 11 Oct 2022

%colormap(bone) % Uncomment for Pluto :)

p=@rand;

q='no';

a=500;

u=ones(a);

u(p(a)>.98)=132;

image(u);

camva(4)

axes(Color=q,CameraP=[0,2,-7],CameraT=[0,1,8],Pr='p')

hold;

Current plot held

% Surface code from Adam Danz (orig. remix)

% Rocky surface

% This is an inverted super-gaussian + noise

% www.mathworks.com/matlabcentral/answers/575647#answer_475645

x=meshgrid(-9:.7:9);

s=-2*exp(-(x.^2/9).^3)+2.1;

% This magic from the legend Jenny Bosten

d=-200:.801:200;

surf(x,s.*p(26),x',FaceC='texturemap',EdgeC=q,CData=abs(ifft2(abs(d+d'*i).^-1.6.*cos(7*p(a)))))

% Set equal aspect ratio

axis equal