matlab用泰勒展开解微分方程,mathematica的解微分方程的能力让人大失所望啊

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

Clear["Global`*"]

c = 299792458*10^2(*光速，单位cm/s*)

G = 6.67259*10^-8(*gravitational constant，引力常数，单位cm^3/g*s^2*)

Msun = 1.9891*10^33(*Subscript[M, \[CircleDot]]，太阳质量，单位g*)

Itilder = 0.283(*Overscript[I, ~]，(7)式下面*)

Jtilder = 1.81*10^-2(*Overscript[J, ~]，(23)式下面*)

Mbi = 1.4*Msun(*Subscript[M, b,i]，(29)式下面*)

\[Eta] = 0.01(*\[Eta]，fig1的不同情况*)

t0 = 3000(*Subscript[t, 0]，fig1的不同情况*)

B = 2*10^14(*(29)式下面，单位G*)

R115 = 1(*Subscript[R, 11.5]*)

R = 11.5*10^5(*单位cm，R=11.5km，(29)式下面*)

M14 = 1(*Subscript[M, 1.4]*)

B15 = 1(*Subscript[B, 15]*)

Bt14 = 1(*Subscript[B, t,14]*)

T9 = 1(*Subscript[T, 9]*)

u = B*R^3(*\[Mu]，磁偶极