2017-04-03から1日間の記事一覧

Wolfram / 三次元勾配ベクトル場

f=x^2+y^2+z^2; v=VectorPlot3D[gradf,{x,-1,1},{y,-1,1},{z,-1,1}, VectorStyle -> "Segment"]; c=ContourPlot3D[f,{x,-1,1},{y,-1,1},{z,-1,1},Mesh->None,ContourStyle->Opacity[0.5,Blue]]; Show[v,c] 参考: レクチャーズオンMathematica, pp.144-145

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] 参考: レクチ…