Borrar filtros
Borrar filtros

Solve differential equation using ode23 and ode45

2 visualizaciones (últimos 30 días)
emonhossain roy
emonhossain roy el 24 de Jun. de 2019
Comentada: Pullak Barik el 25 de Jun. de 2019
Question: solve the differential equation regarding y as the independent variable rather than x. Also find solution for x(0)=1. Finally plot y vs x(y). Also use ode23 and ode45 to solve the equation and plot again.
My approach:
code:
syms x(y);
eqn=diff(y,x)==(x+y*exp(y));
cond=x(0)==1;
dsolve(eqn,cond)
But now how to use ode23 or ode45 and how to plot? Any solution will be appreciated .
Thanks in advance .
  3 comentarios
emonhossain roy
emonhossain roy el 24 de Jun. de 2019
@Star Strider I edited my post
Torsten
Torsten el 25 de Jun. de 2019
Editada: Torsten el 25 de Jun. de 2019
No, you didn't since you still write
eqn=diff(y,x)==(x+y*exp(y));

Iniciar sesión para comentar.

Respuesta aceptada

Pullak Barik
Pullak Barik el 25 de Jun. de 2019
Firstly, as correctly mentioned in the comments, the eqn variable should be written as-
eqn = diff(x, y) == (x + y*exp(y));
Now, for plotting the solution, you can use fplot- fplot documentation
xSol = simplify(dsolve(eqn, cond)); %simplify, as the name suggests, simplifies the expression obtained after dsolve
fplot(xSol);
Now, about using ode23 or ode45, I suggest you to refer to the documentation to gain an understanding on the same. I am posting a code below which uses ode45 (referred directly from the documentation)-
[y, x] = ode45(@(y, x) (x + y*exp(y)), [-10 10], 1)
Link to ode45 documentation- ODE45
  1 comentario
Pullak Barik
Pullak Barik el 25 de Jun. de 2019
You can use plot(y, x) after using the ode45 code above to plot the result

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Ordinary Differential Equations en Help Center y File Exchange.

Community Treasure Hunt

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

Start Hunting!

Translated by