Sunset in the Matrix

on 7 Oct 2021

% Create a 1000x1000 figure

fig1=figure('Position',[10,10,1000,1000]);

% Specify the basic plot parameters

plot(0,0);

xlim([-1,1]);

ylim([-1,1]);

axis off

hold on

% Draw the image

for i=0.5:.5:360 % For cyclical functions, you'll want to iterate across a

% full circle, which is 360 degrees.

% For these four lines, choose any trigonometric functions you like

xt=sind(i);

yt=cosd(i);

xt2=.9*sind(3*i);

yt2=.9*cosd(3*i);

xt3=.3*sind(5*i);

yt3=.3*cosd(5*i);

%p01=plot(sind(0:.5:360),cosd(0:0.5:360),'b-');

%p02=plot(.7*sind((0:.5:360)*2),.7*cosd((0:0.5:360)*2),'r-');

% Plot each line of the drawing

p=plot([xt,xt2],[yt,yt2],'Color',[0 0 1]);

p2=plot([xt2,xt3],[yt2,yt3],'Color',[1 0 0]);

end