用MATLAB求微分方程的奇异性

Question

我无法通过ode45解决这个微分方程，因为它具有奇点。

xy"=3xcos(x)+sin(x) ; x(0)=0 , x'(0)=0

你可以帮我写ode45函数吗？

Answer 1

您可以使用sinc(x)函数，该函数定义为sin(π*x)/(π*x)，但x=0除外，其值为1.因此，您可以将ODE重写为：

y'' = 3*cos(x) + sinc(x/π)

ode45不应该解决任何问题。