daytime

on 8 Oct 2022

% Set up axes

c=[1 1 0];

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

hold

Current plot held

% 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;

surf(x,s.*rand(26),x')

shading interp

% Set equal aspect ratio

axis equal

% Sun

t=0:.1:6;

fill3(2*cos(t),2*sin(t)+1,0*t+9,'.y')

% Horiz sun lines

plot3([-2,2],[1;1]*(-1:.1:3),[9,9],Col=c);

% Add lighting

light(Po=[0,9,30],Col='y')

camlight