Function as an integrator

1 次查看（过去 30 天）
Hugo Hernández Hernández
评论： Hugo Hernández Hernández 2021-1-25
Hi everybody,
I'm computing a numerical solution by Euler's method for a second order differential equation, I have defined the function for the differential equation but I also need to define an additional one as an integrator (like ode45), but by Euler's method. I've tried some computations and agorithms but it does not work. Please, could you help me?
This is my code:
% Sampling interval
h = 0.1;
%h = 0.05;
%h = 0.01;
t0 = 0; % Initial time
tf = 4; % Final time
t = [t0:h:tf]; % Time vector
% Initial value problem
th0 = pi/40; % theta(t0)
om0 = 0; % ommega(t0)
yu = [th0 om0];
e = euler(@yprime,t,h,yu); % Want to call euler's function as an integrator but does not work and has an issue with * multiplication
function [yp] = yprime(t,y)
L = 0.61;
G = 9.81;
yp = zeros(2,1);
yp(1) = y(2); % Angle theta [rad]
yp(2) = (-G/L)*sin(y(1)); % Angular velocity omega [rad/s]
end
% and my Euler's function as a script:
function eu = euler(yp,t,h,y)
% Function that computes the values of the angle theta followinf the Euler
% method
eul = y + yp*h;
eu = eul;

请先登录，再回答此问题。

采纳的回答

Walter Roberson
Walter Roberson 2021-1-24
eul = y + yp(t,y)*h;
  6 个评论
显示 4更早的评论
Walter Roberson
Walter Roberson 2021-1-24
eu1(i+1) = y(i)+h*yp(t(i));
That line of code does not appear in anything you have posted so far. It is failing to pass y values into yp. It is also trying to work with one y value at a time, which is incorrect as yp needs both y values and returns two outputs.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

产品


版本

R2020a

Community Treasure Hunt

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

Start Hunting!

Translated by