MATLAB: 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.

differential transcendental equation

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.

Best Answer

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
Related Question