I am stuck in a question

my code is

syms y(t);

ode = diff(y,t) == (1-y)*exp(-0.1*y)

ysol(t) = dsolve(ode)

cond = y(0) == .5;

ysol(t) = dsolve(ode,cond)

t=[0:5]

plot(t,cell2mat(double(ysol(t))))

how would i determine the stationary solution in matlab when t approaches to infinity

Answer 1

[David Goodmanson]

Hello Arslan, (I checked this by ode45 but the principle should be the same). You are solving dy/dt = (1-y)*exp(-.1*y). For a stationary solution, dy/dt = 0. So either exp(-.1*y) = 0 or (1-y) = 0. The first case corresponds to y = inf which is probably not what you are looking for, and the second case should be the solution.