level 1
君一临
楼主
s = NDSolve[{r[t] \[Theta]''[t] +
2 r'[t] \[Theta]'[t] == -5*10^12 \[Theta]'[t]/r[t],
r''[t] - r[t] \[Theta]'[t]^2 == -4/3*10^20/r[t]^2,
r[0] == 10^11, \[Theta][0] == 0,
r'[0] == 0, \[Theta]'[0] == 3.65*10^-7}, {r, \[Theta]}, {t, 45.8}]###求微分方程数值解
Plot[Evaluate[r[t] /. s], {t, 0, 45}, PlotRange -> All]###上限为45s时作图
Plot[Evaluate[r[t] /. s], {t, 0, 46}, PlotRange -> All]###上限为46s时作图
今天第一次用mathmetica，麻烦大家了。按理说改变数值计算时间上限不会改变函数关系啊，这是怎么回事呢？
2019年12月25日 08点12分
1
2 r'[t] \[Theta]'[t] == -5*10^12 \[Theta]'[t]/r[t],
r''[t] - r[t] \[Theta]'[t]^2 == -4/3*10^20/r[t]^2,
r[0] == 10^11, \[Theta][0] == 0,
r'[0] == 0, \[Theta]'[0] == 3.65*10^-7}, {r, \[Theta]}, {t, 45.8}]###求微分方程数值解
Plot[Evaluate[r[t] /. s], {t, 0, 45}, PlotRange -> All]###上限为45s时作图
Plot[Evaluate[r[t] /. s], {t, 0, 46}, PlotRange -> All]###上限为46s时作图
今天第一次用mathmetica，麻烦大家了。按理说改变数值计算时间上限不会改变函数关系啊，这是怎么回事呢？
