解決在ODE45 -Matlab

Question

我有一個ODE，y'=y^2+y;所以，我寫了一個腳本：

foo=@(y)(y.^2+y); 
[x y]=ode45(foo,[1 4],1); 

但它返回以下錯誤：

Error using @(y)(y.^2+y) 
Too many input arguments. 

Error in odearguments (line 88) 
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} 
to yp0. 

Error in ode45 (line 114) 
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, 
odeArgs, odeFcn, ... 

我必須做出定義函數的一些錯誤。 我很欣賞任何意見和建議。

Answer 1

引述official documentation：

All solvers solve systems of equations in the form y′ = f(t,y) or problems that involve a mass matrix, M(t,y)y′ = f(t,y).

你的函數只接受一個變量y，而它必須接受兩個：t和y。所以，如果您的ODE是y′ = y2 + y，按以下方式確定foo：

foo = @(t, y)(y .^ 2 + y); 

，它應該工作。