ode45 not enough input arguments error

trying to solve a 2nd order ODE using ode45, i looked up several tutorials to help me as i was working through it. My code looks very similar to the tutorials with the exception of different equations and initial values.

function rk=ode45bessel(t,v)
rk=[v(2);-v(1)-(1/t)*v(2)];

t=[0:50];

initialvalues=[-4.0690e-05,-8.1380e-05];

[t,v]=ode45(@ode45bessel,t,initialvalues);

plot(t,v)

thats my code and i get an error on line 4 saying not enough input arguments

Answer 1

[Walter Roberson]

You need to break that into two files

function rk=ode45bessel(t,v)
rk=[v(2);-v(1)-(1/t)*v(2)];

which you store into ode45bessel.m

and the rest of the lines in a second file (with a name other than 'ode45bessel.m' or 'ode45.m') . Then you would run the second file.

[Daniel Hunt]

thanks, this is now running correctly but giving me the ouput NaN (not a number) for most of my values of v, any idea why this could be?

[Walter Roberson]

You start at t = 0, and you have 1/t in your ode45bessel, which is leading to an infinite value there that forces the second element of the first response to be infinite. That infinite value then comes back to the very next iteration and gets carried over the first output, so you get [inf -inf] for the second iteration. That gives you -inf+inf to be integrated for the second element, which gives nan. You are lost from there.

[Daniel Hunt]

didnt even think of that, all sorted now