function lorenz
  g=10; // γの値
  b=8/3; // bの値
  mu=28; // μの値
  dt=0.01; // 時間刻み幅 Δt
  nstep=10000; // 実行するステップ数
  
  x=zeros(nstep+1);
  y=zeros(nstep+1);
  z=zeros(nstep+1);
  
  x(1)=0; // x の初期値
  y(1)=1; // y の初期値
  z(1)=0; // z の初期値

  for n=2:nstep+1
    x(n)=x(n-1)-dt*g*(x(n-1)-y(n-1));
    y(n)=y(n-1)+dt*(mu*x(n-1)-y(n-1)-x(n-1)*z(n-1));
    z(n)=z(n-1)+dt*(-b*z(n-1)+x(n-1)*y(n-1));
  end

  param3d(x,y,z)
  param3d1(x(nstep),y(nstep),list(z(nstep),-4))
  
