حرکت پرتابه در متیمتیکا
Cear[x,v,t,a1,a2,f1,f2
S=Dsolve[{x''[t]
Š-x'[t]-x'[t]^2,y''[t]Š-y'[t]+y'[t]^2,z''[t]Š-z'[t]-z'[t]^2-g,x[0]Š0,y[0]Š0,z[0]Š0,x'[0]Švx0,y'[0]Švy0,z'[0]Švz0},{x[t],y[t],z[t]},t]D=Solve[{5050*t^2+(5-(Sin[f])*Cos[a])*t+1
Š0,50.55*t-Cos[f]Š0,49.50*t+Sin[f]*Sin[a]Š0/.s,},{t,f,a}]V={vxo
®100,vy0®100,vz0®100,g®10};ListPlot[Evaluate[{x[t],y[t],z[t]}/.s/.v],{t,0,1},Frame
®True,FrameLabel®{"x[t]","y[t]","z[t]"},RotateLabel®False]
+ نوشته شده در چهارشنبه نوزدهم تیر ۱۳۸۷ ساعت 3:11 توسط
|