need plot a graph using Euler's method for both numerical and analytical
조회 수: 4 (최근 30일)
이전 댓글 표시
% Euler's Method
% Initial conditions and setup
h = (10); % step size
x = (0):h:(60); % the range of x
y = zeros(size(x)); % allocate the result y
y(1) = (2); % the initial y value
n = numel(y); % the number of y values
for i=1:n-1
pi = 3.1416
a = 2
r = 1
f = a + (2 - 0.4*a) * x / r^2 * pi
y(i+1) = y(i) + h * f;
end
plot (x,y,'b')
grid on
댓글 수: 4
Jan
2022년 10월 30일
@Mpho Mokgotlo: If you mention an error, post a copy of the complete error message. Seeing some parts of the message only is less useful, e.g. the readers do not know, in which line the error occurs.
답변 (2개)
Alan Stevens
2022년 10월 30일
편집: Jan
2022년 10월 30일
f = a + (2 - 0.4*a) * x(i) / r^2 * pi; % x(i) not just x.
% Also take note of the other comments above.
댓글 수: 1
Jan
2022년 10월 30일
h = 10; % step size
x = 0:h:60; % the range of x
y = zeros(size(x)); % allocate the result y
y(1) = 2; % the initial y value
n = numel(y); % the number of y values
a = 2;
r = 1;
for i = 1:n-1
f = a + (2 - 0.4 * a) * x(i) / r^2 * pi;
% ^^^
y(i+1) = y(i) + h * f; % I guess this line fails
end
plot (x,y,'b')
grid on
A cleaned version of your code. Alan Stevens' hint is included (please accept his answer).
You can solve such problems using the debugger. Type in the command window:
dbstop if error
and run the code again. Matlab stops at the problem and you can check the values of the used variables.
size(y(i+1)) % A scalar, of course
size(y(i) + h * f) % Not a scalar
size(f) % Same size as x!
You do not want x in the line to get the value of f, but only the element x(i).
Do you really mean: ... / r^2 * pi or .../ (r^2 * pi) ?
댓글 수: 0
참고 항목
카테고리
Help Center 및 File Exchange에서 Graphics Performance에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!