Euler method for vectors?

5 次查看（过去 30 天）

qwerty ed 2015-4-8

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/196675-euler-method-for-vectors

编辑： James Tursa 2015-4-8

Here is code I wrote for computing Euler's method. My question is, how can I make this work for vector inputs? (I have another function that returns a vector equation). Thank you.

EDIT: I'm not sure why the formatting is messed up and I'm not quite sure how to fix it.

_______________________________________

function euler(func, y0, t0, tf)
% y0 = initial y value
% m = number of spatial discretization points
% t0 = initial time (t0 = 0)
% tf = final time (tf = 4)
h = (2*pi)/32; % define step size
t=t0:h:tf; % time interval
% setting initial y value
y(1) = y0;
% loop using euler's method
for i = 1:length(t)-1
    y(i+1) = y(i) + h*(feval(func,t(i)));
end
% print column vectors of t and y
t = t'
y = y'
plot(t,y)
xlabel('time')
ylabel('y')