i want to plot a equation dy/dt=Sqrt ((1-1/y^4)+log(y)). What type of ode solver should i use. The initial value of y at t=0 is 1.005. Also the limits of t are 0 to 0.05.

function bb [t,y]=ode(@odefun,[0 0.05],1.005); plot(t,y(:,1),'-o') end

function dy =odefun(t,y) a=250; b= 1e5; dy=(a*(1-1/y(1)^4)+b*log(y(1))); end

This code is not solving my problem.

Answer 1

[Star Strider]

Your ODE function .m file ‘odefun’ will work. (I used an anonymous function for your ODE, because it’s more convenient for me.) You can use ode45 to integrate your ODE:

odefun = @(t,y) (250*(1-1./y(1).^4)+1E+5.*log(y(1)));

[t,y]=ode45(odefun,[0 0.05],1.005); 

figure(1)
plot(t,y(:,1),'-o')
grid

[sv]

Thanks Sir for your answer. This was really helpful.

[Star Strider]

My pleasure!