Euler's Method

24 visualizaciones (últimos 30 días)
Matthew Russell
Matthew Russell el 10 de Jun. de 2019
Respondida: Shagufta Perveen el 24 de Mayo de 2021
Hello,
New Matlab user here and I am stuck trying to figure out how to set up Euler's Method for the following problem:
?′ =sin(?)∗(1−?) with ?(0)=?0 and ?≥0
The teacher for the class I am taking provided us with the following code to use for Euler's Method. The code has been modified with my most recent attempt to solve the problem above. His template worked fine for the problem listed at the top and several others, but when I changed out the variables for the problem above and I get the error message listed at the bottom. I have tried various ways of inputing the code that I found on YouTube or Google searches and none have been able to help. What am I doing wrong?
%This code solves the differential equation y' = 2x - 3y + 1 with an
%initial condition y(1) = 5. The code uses
%the Euler method, the Improved Euler method, and the Runge-Kutta method.
%The function f(x,y) = 2x - 3y + 1 is evaluated at different points in each
%method.
h = 1/16; %Time Step
a = 0; %Starting x
b = 20; %Ending x
n = 321; %Number of Iterations
x = zeros(n,1);
y = zeros(n,1);
x = linspace(a,b,n)'; %Array of x values where evaluate the function
y(1) = 6; %Initial Condition
for i = 1:n-1
y(i+1) = y(i) + h *((sin(x) * ( 1 - y(i)) ; %Euler's Method
end
[x y] %Table showing x and y values
Error: File: Euler_Method.m Line: 21 Column:
47
Invalid expression. When calling a function
or indexing a variable, use parentheses.
Otherwise, check for mismatched delimiters.

Respuesta aceptada

Debasish Samal
Debasish Samal el 10 de Jun. de 2019
There is a parentheses mismatch in your code for the euler's method.
Replace that line with:
y(i+1) = y(i) + h *(sin(x) * ( 1 - y(i))) ; %Euler's Method
  1 comentario
Matthew Russell
Matthew Russell el 11 de Jun. de 2019
Thank you so much. I had to change sin(x) to sin(x(i)) for it to work, but it worked perfect after that.

Iniciar sesión para comentar.

Más respuestas (2)

Raghunandan V
Raghunandan V el 10 de Jun. de 2019
Hi,
Simple code error. Check this line
y(i+1) = y(i) + h *((sin(x(i)) * ( 1 - y(i)))) ; %Euler's Method

Shagufta Perveen
Shagufta Perveen el 24 de Mayo de 2021
y = te3t − 2y, 0 ≤ t ≤ 1, y(0) = 0, with h = 0.5
if true
% code
end

Categorías

Más información sobre Error Functions en Help Center y File Exchange.

Etiquetas

Productos


Versión

R2019a

Community Treasure Hunt

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

Start Hunting!

Translated by