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.

2 次查看（过去 30 天）
sv
sv 2014-8-17
评论： Star Strider 2014-8-18
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.

请先登录，再回答此问题。

采纳的回答

Star Strider
Star Strider 2014-8-17
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

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by