Question

使用mathematica解常微分方程（数值解）

方程是
t为自变量，c、k、h为常数。取正号。

Answer 1

这个方程可以求符号解啊，虽然答案里用到了特殊函数，但是稍微带入下具体数值求个数值解也是可以的：
eqn = r'[t]^2 == c r[t]^2 + 2 k^2 r[t] - h^2/r[t];
DSolve[eqn, r[t], t]

假如想直接来数值解就要用NDSolve，当然还要补相应的边界条件和具体参数，语法与上面的大同小异，看看帮助吧。