clear all
global r b s c N

N=6;
r=28;
b=8/3;
s=10;
c=[1.5,0,0];


ee=1.e-03;
 psi0=[1;0.3;0.5];
 psi00=[-0.4252;-1.7460;19.2357];
 psi01=randn(3,1);
 %psi01=[-0.4252+ee;-1.7460+2*ee;19.2357-ee/10];
 psi0=[psi00;psi01];
  
  
  
  options=odeset('RelTol',1e-5,'AbsTol',1.e-06*ones(N,1));
  dt=0.01;
  tmax=50;
  Nt=floor(tmax/dt);
  Psi=zeros(N,Nt);
  Psi(:,1)=psi0;
  tt=zeros(1,Nt);
  
  for it=2:Nt;
  tt(it)=(it-1)*dt;
  tspan=[tt(it-1),tt(it)];
  [TT,UUt] = ode45(@dlorenz,tspan,Psi(:,it-1),options);
  Psi(:,it)=(UUt(length(TT),:))';
  end;
  
  xt=Psi(1,:);yt=Psi(2,:);zt=Psi(3,:);
  x1t=Psi(4,:);y1t=Psi(5,:);z1t=Psi(6,:);
  figure(3);
  plot3(xt,yt,zt,'b')
  hold on
  plot3(x1t,y1t,z1t,'r')
  xlabel('x');ylabel('y'),zlabel('z');
  hold off
  
  
  figure(4);
  plot(tt,xt,'b',tt,yt,'r',tt,zt,'k',tt,0*tt)
  hold on
  title('x : blue , y: red , z: red')
  plot(tt,x1t,'--b',tt,y1t,'--r',tt,z1t,'--k',tt,0*tt)
  title('x : blue , y: red , z: red')
  hold off