Question

求大神指点！matlab编写了一段龙格库塔法求解微分方程，可是为什么运行了会显示输入的参数不足啊

function dy=ofun(t,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=-(y(1)-6)*y(2)+2*y(1);
t0=0;
tp=20;
y0=[0;0.25];
[t,Y]=ode45(@ofun,[t0,tp],y0);
plot(t,Y(:,1),'r')
hold on;
plot(t,Y(:,2),'g')
hold off;
legend('y','y‘');
xlabel('time t');
ylabel('Y');

保存运行后弹出警告：
输入参数的数目不足。
dy(1)=y(2);

Answer 1

微分在数学中的定义：由函数B=f(A)，得到A、B两个数集，在A中当dx靠近自己时，函数在dx处的极限叫作函数在dx处的微分，微分的中心思想是无穷分割。微分是函数改变量的线性主要部分。微积分的基本概念之一。
设函数y = f(x)在x的邻域内有定义，x及x + Δx在此区间内。如果函数的增量Δy = f(x + Δx) - f(x)可表示为 Δy = AΔx + o(Δx)（其中A是不随Δx改变的常量，但A可以随x改变），而o(Δx)是比Δx高阶的无穷小（注：o读作omicron，希腊字母）那么称函数f(x)在点x是可微的，且AΔx称作函数在点x相应于因变量增量Δy的微分，记作dy，即dy = AΔx。函数的微分是函数增量的主要部分，且是Δx的线性函数，故说函数的微分是函数增量的线性主部（△x→0）。
通常把自变量x的增量 Δx称为自变量的微分，记作dx，即dx = Δx。于是函数y = f(x)的微分又可记作dy = f'(x)dx。函数因变量的微分与自变量的微分之商等于该函数的导数。因此，导数也叫做微商。
当自变量X改变为X+△X时，相应地函数值由f(X)改变为f(X+△X)，如果存在一个与△X无关的常数A，使f(X+△X)-f(X)和A·△X之差是△X→0关于△X的高阶无穷小量，则称A·△X是f(X)在X的微分，记为dy，并称f(X)在X可微。一元微积分中，可微可导等价。记A·△X=dy，则dy=f′(X)dX。例如：d(sinX)=cosXdX。
微分概念是在解决直与曲的矛盾中产生的，在微小局部可以用直线去近似替代曲线，它的直接应用就是函数的线性化。微分具有双重意义：它表示一个微小的量，因此就可以把线性函数的数值计算结果作为本来函数的数值近似值，这就是运用微分方法进行近似计算的基本思想。
希望我能帮助你解疑释惑。