Error "Not enough Input Arguments" in predator-prey model

Question

summyia qamar 2018-1-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/375483-error-not-enough-input-arguments-in-predator-prey-model

评论： Star Strider 2018-1-4

I'm trying to understand the working of ODE45 to solve multiple variable equations. I copied a code from internet but when I run this code, matlab gives error

Not enough input arguments.
Error in predPrey (line 2)
f=[-x(1)+3*x(3);-x(2)+2*x(3);x(1)^2-2*x(3)];
function f=predPrey(t,x);

the code is given here:

        f=[-x(1)+3*x(3);-x(2)+2*x(3);x(1)^2-2*x(3)];
          [t,xa] = ode45(@predPrey,[0 1.5],[0 1/2 3]);
          plot(t,xa(:,2))
          title('y(t)')
          xlabel('t'), ylabel('y')

can someone please help me in this error so I can continue modifying model for my system.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2018-1-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/375483-error-not-enough-input-arguments-in-predator-prey-model#answer_298617

在 MATLAB Online 中打开

You did not copy it correctly or implement it correctly.

Try this:

predPrey =  @(t,x) [-x(1)+3*x(3);-x(2)+2*x(3);x(1)^2-2*x(3)];
[t,xa] = ode45(predPrey,[0 1.5],[0 1/2 3]);
plot(t,xa(:,2))
title('y(t)')
xlabel('t'), ylabel('y')

3 个评论
显示 1更早的评论隐藏 1更早的评论

Walter Roberson 2018-1-4

在 MATLAB Online 中打开

plot(t, xa)

Star Strider 2018-1-4

@summyia qamar — My pleasure.

@Walter — Thank you.

请先登录，再进行评论。

Answer 2

Walter Roberson 2018-1-3

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/375483-error-not-enough-input-arguments-in-predator-prey-model#answer_298618

在 MATLAB Online 中打开

The way you give your code makes it appear as if all of your code is within predPrey.m and that your call to ode45 is on line 3 of predPrey.m and that you try to run the ode by executing the predPrey file. If so then you need to break out the line starting from [t,xa] into a different file. Or, if you are running R2016b or later, you could use one file that is not named prepPrey.m that contains

          [t,xa] = ode45(@predPrey,[0 1.5],[0 1/2 3]);
          plot(t,xa(:,2))
          title('y(t)')
          xlabel('t'), ylabel('y')
          function f=predPrey(t,x)
          f=[-x(1)+3*x(3);-x(2)+2*x(3);x(1)^2-2*x(3)];

The alternative to all of this is to use just one file containing:

          prePrey = @(t, x) [-x(1)+3*x(3);-x(2)+2*x(3);x(1)^2-2*x(3)];
          [t,xa] = ode45(predPrey,[0 1.5],[0 1/2 3]);
          plot(t,xa(:,2))
          title('y(t)')
          xlabel('t'), ylabel('y')