用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應該沒有任何問題解決。