"iam49erfreak " <firstname.lastname@example.org> wrote in message news:email@example.com... > I am trying to solve a system of ODEs and when I follow an example on how > to use ode45 I keep getting an error saying that I don't have enough input > arguments and that 't' isnt being used. What am I doing wrong? Any help is > greatly appreciated! > > function dy = rigid(t,y) > dy = zeros(3,1); > dy(1) = (16.38124)*(y(2) - ((y(1))*(y(2))) + y(1) - (.008376*((y(1))^2))); > dy(2) = (.0027436)*(y(3) - y(2) - (x*(y(2)))); > dy(3) = (.034132)*(y(1) - y(3)); > > options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]); > [T,Y] = ode45(@rigid,[0 610],[4 1.1 4],options);
Move the ODE45 call with @rigid as its first input OUT of the rigid function itself. It should be in a separate function or script.