i am trying to solve this differential equation and the solution i get it is wrong so where is the error? and this my function dy/dt=k*y*ln(450/y) y(0)=40 and y(15)=95 and k is constant and this my MATLAB
tspan = [0 300];
y0 = 40;y15=95;
[t,y] = ode45(@(t,y) 'k'*y*exp(450/y), tspan, y0,y15);
plot(t,y,'-o')

 Réponse acceptée

Takey Asaad
Takey Asaad le 14 Sep 2018

0 votes

Hi , the plot is still wrong , k should be get from the initial condition and the solution is y(t)=450/( e^e^(-0.02948t+0.8839) so how can i get this by MATLAB

Plus de réponses (1)

Torsten
Torsten le 14 Sep 2018
Modifié(e) : Torsten le 14 Sep 2018

0 votes

This is a boundary value problem - use "bvp4c" instead of "ode45".
Or try "dsolve".
Best wishes
Torsten.

Community Treasure Hunt

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

Start Hunting!

Translated by