Euler's method in matlab

4 visualizzazioni (ultimi 30 giorni)
Haya Ali
Haya Ali il 2 Set 2021
Modificato: Haya Ali il 8 Set 2021
If I write the model using Euler's method. Like this
clear all; close all; clc;
%value of constants
a_n=0.1;
omega_0=2;
dt=0.01; %step size
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x(1)=1;
y(1)=1;
for i=2:1000
x(i)=x(i-1)+((a-x(i-1)^2-y(i-1)^2)*x(i-1)-omega_0*y(i-1))*dt;
y(i)=y(i-1)+((a-x(i-1)^2-y(i-1)^2)*y(i-1)+omega_0*x(i-1))*dt;
end
Does that mean it will give me values of x and y instead of dx/dt and dy/dt?

Risposta accettata

Fabio Freschi
Fabio Freschi il 2 Set 2021
Modificato: Fabio Freschi il 2 Set 2021
Any ODE solver gives the unknown function(s) as result, not its derivative. So, yes, your results are x and y.
Note that in your solution you forgot the term

Più risposte (0)

Categorie

Scopri di più su Programming in Help Center e File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by