Euler's Method plotting
Show older comments
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 Comments
Youssef Khmou
on 15 Dec 2014
can you briefly describe the input variables.
Amit
on 15 Dec 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?
Answers (0)
Categories
Find more on Mathematics in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!