Euler's Method plotting
Afficher commentaires plus anciens
Hi,
i found this code for the euler's method on the internet:
function [wi, ti] = euler ( RHS, t0, x0, tf, N )
neqn = length ( x0 );
ti = linspace ( t0, tf, N+1 );
wi = zeros( neqn, N+1 ) ;
wi(1:neqn, 1) = x0';
h = ( tf - t0 ) / N;
for i = 1:N
x0 = x0 + h * feval ( RHS, t0, x0 );
t0 = t0 + h;
wi(1:neqn,i+1) = x0';
end
% RHS string containing name of m-file defining the
% right-hand side of the differential equation; the
% m-file must take two inputs - first, the value of
% the independent variable; second, the value of the
% dependent variable
% t0 initial value of the independent variable
% x0 initial value of the dependent variable(s)
% if solving a system of equations, this should be a
% row vector containing all initial values
% tf final value of the independent variable
% N number of uniformly sized time steps to be taken to
% advance the solution from t = t0 to t = tf
%
% output:
% wi vector / matrix containing values of the approximate
% solution to the differential equation
% ti vector containing the values of the independent
% variable at which an approximate solution has been
% obtained
%
Now, I want to plot the result of the function. Unfortunately I wasn't able to solve this problem by my own.
I tried: plot(xo,to) and many other variations of this, but the result didn't make any sense.
Many thanks in advance, Fabian
2 commentaires
Youssef Khmou
le 15 Déc 2014
can you briefly describe the input variables.
Amit
le 15 Déc 2014
First question is: Do you know Euler's method? (if yes, it is a rather simple code and I think you should try it on your own - it will be a good exercise)
Second question is: where is your attempt or what have you done to understand this code?
Réponses (0)
Catégories
En savoir plus sur Mathematics dans Centre d'aide et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!