Mathematica求非线性微分方程，用Python求解还算容易MMA却怎么都出不来，请教如何解决？

首先你Python用的是数值求解，但是你在MMA中却使用DSolve进行符号求解计算，这怎么看都不公平。

首先MMA用数值求解NDSolve 100%能算

Block[{g=98/10,L=1,sol},
sol=NDSolve[{(\[Theta]^\[DoublePrime])[t]+(3*g)/(2*L) Sin[\[Theta][t]]==0,(\[Theta]^\[Prime])[0]==0,\[Theta][0]==\[Pi]/6},\[Theta],{t,0,10}];
Plot[Evaluate[{\[Theta][t],-Integrate[\[Theta][x],{x,0,t}]}/.sol],{t,0,10}]
]

然后这个符号解其实也是可以算的(用mma的特殊函数表示)

Block[{g=98/10,L=1},DSolve[{(\[Theta]^\[DoublePrime])[t]+(3*g)/(2*L) Sin[\[Theta][t]]==0,(\[Theta]^\[Prime])[0]==0,\[Theta][0]==\[Pi]/6},\[Theta][t],t]//FullSimplify]

但是出现了反函数求解的通病问题，就是间断，我们直接对上面这个解画图，会是这样的

然后查询了下帮助JacobiAmplitude这个函数的帮助文档，在可能存在问题中，文档里给出的解决方案

于是，我们修改一下，就可以得到