読者です 読者をやめる 読者になる 読者になる

Wolfram / ローレンツアトラクター

sol=NDSolve[{x'[t]==10(y[t]-x[t]),y'[t]==x[t](28-z[t])-y[t],z'[t]==x[t]y[t]-(8/3)z[t],x[0]==1,y[0]==0,z[0]==0},{x[t],y[t],z[t]},{t,0,50}];
p[t_]={x[t],y[t],z[t]}/.sol[[1]];
lorenz=ParametricPlot3D[p[t],{t,0,50},PlotRange->All]

f:id:ti-nspire:20170403105325p:plain
参考: レクチャーズオンMathematica, p.141-142

広告を非表示にする