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')

 채택된 답변

Takey Asaad
Takey Asaad 2018년 9월 14일

0 개 추천

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

추가 답변 (1개)

Torsten
Torsten 2018년 9월 14일
편집: Torsten 2018년 9월 14일

0 개 추천

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

카테고리

질문:

2018년 9월 14일

댓글:

2018년 9월 14일

Community Treasure Hunt

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

Start Hunting!

Translated by