如何解决颂歌

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代替of y / t和（0,1）。