如何解决代码

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）。