Back to Gallery
  • Remix
  • Share
  • New Entry
on 7 Nov 2023
  • 12
  • 47
  • 1
  • 0
  • 802
drawframe(1);
Write your drawframe function below
function drawframe(f)
h = 0.01;
a = 10;
b = 28;
c = 8/3;
persistent x y z ColorData
if isempty(x)
[X,Y,Z] = meshgrid(-10:1:10);
X2 = []; Y2 = []; Z2 = [];
for n = 1:length(X)
for m = 1:length(X)
X2 = [X2 X(n,:,m)];
Y2 = [Y2 Y(n,:,m)];
Z2 = [Z2 Z(n,:,m)];
end
end
ColorData = zeros(size(X2));
for i = 1:length(X2)
ColorData(i) = norm([X2(i) Y2(i) Z2(i)]);
end
x = X2;
y = Y2;
z = Z2;
end
for j = 1:5
dx = -a*x + a*y;
dy = -x.*z + b*x -y;
dz = x.*y - c.*z;
kx1 = dx;
kx2 = dx + h*kx1/2;
kx3 = dx + h*kx2/2;
kx4 = dx + h*kx3;
x = x + h*(1/6)*(kx1+2*kx2+2*kx3+kx4);
ky1 = dy;
ky2 = dy + h*ky1/2;
ky3 = dy + h*ky2/2;
ky4 = dy + h*ky3;
y = y + h*(1/6)*(ky1+2*ky2+2*ky3+ky4);
kz1 = dz;
kz2 = dz + h*kz1/2;
kz3 = dz + h*kz2/2;
kz4 = dz + h*kz3;
z = z + h*(1/6)*(kz1+2*kz2+2*kz3+kz4);
end
colormap jet
scatter3(x,y,z,5,'filled','CData',ColorData);
daspect([1 1 1])
axis([-40 40 -40 40 -10 60])
view(25,15)
axis off
end
Animation
Remix Tree