如何解決代碼

Question

我有一個可以解決dy/dt = t/y和y（0）= 1的matlab代碼，我該如何概括這個來解決dy/dt = f（t，y），y o）= y_0對於任何給定的f（t，y）？

Answer 1

您將不得不閱讀一些有關ODE的數值方法的書籍。我認爲即使對於這個相對簡單的情況，你也不能製作通用的代碼。你在這裏y' + F(y,t) = 0這只是一階ODE。在不瞭解F的結構的情況下，我們無法談論尋找解決方案的單一方式。但是如果F至少是Lipschitz連續的，你可以嘗試任何Euler方法。他們並不困難。我可以想象你的代碼已經使用了歐拉方法之一，所以如果你不關心F是一些「壞」函數，你可以修改你的代碼來使用F和y（0）= y_0而不是y/t和（0,1）。