Converting Euler's method calculation into matlab

Question

Dai Nguyen 2020-9-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/599560-converting-euler-s-method-calculation-into-matlab

评论： Dai Nguyen 2020-9-25

Hi I wanna use Matlab to solver for an Euler calculation with the range of t from 0 to 10 with 0.5 step. however for some reason my code isn't working. The initial y=0

this is the equation dy/dt=(3*Q/A*(sin(x))^2))-(alpha*(1+y)^1.5)/Q

% Euler's Method

% Initial conditions and setup

Q=400

A=1200

alph=160

h = 0.5; % step size

t = 0:0.5:10; % the range of t

y = zeros(size(t)); % allocate the result y

y(1) = 0; % the initial y value

n = numel(y); % the number of y values

% The loop to solve the DE

for i=1:n-1

f = (3*(Q/A)*(sin(t))^2)-((alph*(1+y(i))^1.5)/A)

y(i+1) = y(i) + h * f;

end

plot(x,y); grid on

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Rafael Hernandez-Walls 2020-9-25

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/599560-converting-euler-s-method-calculation-into-matlab#answer_500188

在 MATLAB Online 中打开

tray this:

% Initial conditions and setup
Q=400;
A=1200;
alph=160;
h = 0.5;  % step size
t = 0:0.5:10;  % the range of t
y = zeros(size(t));  % allocate the result y
y(1) = 0;  % the initial y value
n = numel(y);  % the number of y values
% The loop to solve the DE
for i=1:n-1         %   v---- t fot t(i)
    f = (3*(Q/A)*(sin(t(i))).^2)-((alph*(1+y(i)).^1.5)/A);
    y(i+1) = y(i) + h * f;
end
plot(t,y); grid on