I am working at a mathematical problem and I wanted to solve it numerically through MATLAB. The equation has the form:
epsilon*y'(t) =t*y(t)(y(t)-a)(y(t)-b)+c
where a<0 and b>0 and c>0

 채택된 답변

Torsten
Torsten 2019년 2월 7일
편집: Torsten 2019년 2월 7일

0 개 추천

a = -0.2;
b = 0.1;
c = 0.1;
epsilon = 1.0;
fun = @(t,y)(t*y*(y-a)*(y-b)+c)/epsilon;
tspan = [0 3.5];
y0 = 0;
[T,Y] = ode45(fun,tspan,y0);
plot(T,Y)

댓글 수: 4

Thank you for your help, but are you sure for you answer, could you explain more the method of your work, please.
Torsten
Torsten 2019년 2월 7일
Take a look at how ODE integrators are invoked:
help ode45
Would you change this algorithm to Runge Kutta algorithm, pleas?!
a = -0.2;
b = 0.1;
c = 0.1;
epsilon = 1.0;
fun = @(t,y)(t*y*(y-a)*(y-b)+c)/epsilon;
tspan = [0 3.5];
y0 = 0;
[T,Y] = ode45(fun,tspan,y0);
plot(T,Y)
Torsten
Torsten 2019년 3월 4일
ODE45 is a Runge-Kutta integrator ...

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

도움말 센터File Exchange에서 Numerical Integration and Differential Equations에 대해 자세히 알아보기

제품

질문:

2019년 2월 6일

댓글:

2019년 3월 4일

Community Treasure Hunt

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

Start Hunting!

Translated by