- /
-
Scarf
on 16 Oct 2021
- 2
- 17
- 0
- 0
- 279
L=3e5;
rng(14193)
p=-1.2+2.4*rand(1,30);
a=0.05;
b=a;
c=a;
t=zeros(3,L);
for i=1:L
z=[1 a b c a*b a*c b*c a^2 b^2 c^2];
v=[0,0,0];
for k=1:3
v(k)=sum(p((k-1)*10+1:k*10).*z);
end
t(:,i)=[v(1),v(2),v(3)];
a=v(1);
b=v(2);
c=v(3);
end
scatter3(t(1,:),t(2,:),t(3,:),1,t(3,:),'filled')
set(gca,'Color','k')
view(160,0)