MATLAB: Error in ODE45 when using …………………. >> [t,x] = ode45(@odefcn,tspan,[0 0]);

ajmal r sameerluca amerio

function dxdt = odefcn(t,x,a,b,c,u)
dxdt = zeros(2,1);
dxdt(1) = a*x(1)+2*x(2)+u;
dxdt(2) = c*x(1)+ b*x(2)-2*u;
end
clc;clear all; close all;
a=-1; b=-2; c=0; u=3;
tspan = [0 20];
[t,x] = ode45(@odefcn,tspan,[0 0]);
figure(1)
plot(t,x(:,1),'r',t,x(:,2),'b')
i am getting an error (attached screenshot)
while i am not getting error when i am using 
a=-5; b=-3; c=1; u=3;
tspan = [0 20];
[t,x] = ode45(@(t,x) odefcn(t,x,a,b,c,u),tspan,[0 0]);
figure(2)
plot(t,x(:,1),'r',t,x(:,2),'b')
can anyone please explain

Best Answer

[t,x] = ode45(@odefcn,tspan,[0 0]);
That is valid syntax.
However, your function is relying on the variables a, b, c, and u. MATLAB never automatically passes variables from the calling environment into a function. If you define
x = 42;
y = f();
and you have defined
function result = f(x)
  result = x.^2;
end
then when you call f without passing any argument into it, MATLAB does not look and see that the internal parameter name for the first argument is x and then look into the script and see there is an x there and use the caller's x as the value for x in the function. That never happens for named parameters. (There are cases where MATLAB will look for a variable in an outer environment, but that is never done for named parameters.)
So... at the time that your function odefcn first needs to use a, b, c, or u, the function would fail, because you would not have passed the values into the function.
The approach you used with the anonymous function 
@(t,x) odefcn(t, x, a, b, c, u)
is the recommended approach. http://www.mathworks.com/help/matlab/math/parameterizing-functions.html
The approach Alan suggested, of instead passing extra parameters, has not been documented since about MATLAB 5.2 or so, and has some subtle internal incompatibilities so it is not completely reliable.
Related Question