work2.m

v1=V(:,1);v2=V(:,2);v3=10*V(:,3);
plot3(u0(1),u0(2),u0(3),'.r','MarkerSize',10);box
grid;view(-117,20);
xlabel('x');ylabel('y');zlabel('z');
axis([0 20000 0 10000 0 20000]);hold on;
quiver3(u0(1),u0(2),u0(3),v1(1),v1(2),v1(3),1000);
quiver3(u0(1),u0(2),u0(3),v2(1),v2(2),v2(3),1000);
quiver3(u0(1),u0(2),u0(3),-v2(1),-v2(2),-v2(3),1000);
n=64;bas=0.01;M=30;
for i=0:M
    tt=i*pi/M;ttt=(1-cos(tt))/2*pi;
    v=(cos(ttt)*v2+sin(ttt)*v1)*1000;
    alpha=v(1);beta=v(2);gamma=v(3);
    [x,coef,sol]=leuclag1(n,bas,Ps,u0,alpha,beta,gamma);
    plot3(u0(1)+v(1),u0(2)+v(2),u0(3)+v(3),'.g','MarkerSize',10);
    z(1)=u0(1)+sum(coef(1:n));z(2)=u0(2)+sum(coef(n+1:2*n));z(3)=u0(3)+sum(coef(2*n+1:3*n));
    plot3(z(1),z(2),z(3),'.k');
    rez(i+1,:)=z;
    options=odeset('Events',@events_leuc);
    [t,u]=ode45(@leucemia3d,[0,-15000],z,options);
    plot3(u(:,1),u(:,2),u(:,3),'g');
    plot3(u(end,1),u(end,2),u(end,3),'.g');
    i,pause(0.1);
end
[t,u]=ode45(@leucemia3d,[0,25000],[u0(1)+v3(1),u0(2)+v3(2),u0(3)+v3(3)],options);
plot3(u(:,1),u(:,2),u(:,3),'r');
[t,u]=ode45(@leucemia3d,[0,25000],[u0(1)-v3(1),u0(2)-v3(2),u0(3)-v3(3)],options);
plot3(u(:,1),u(:,2),u(:,3),'r');