Function as an integrator
조회 수: 1 (최근 30일)
이전 댓글 표시
Hugo Hernández Hernández
2021년 1월 24일
댓글: 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;
댓글 수: 0
채택된 답변
Walter Roberson
2021년 1월 24일
eul = y + yp(t,y)*h;
댓글 수: 6
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개)
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!