# Euler's method code

조회 수: 1(최근 30일)
Mike Asmanis 2021년 6월 18일
댓글: Mike Asmanis 2021년 6월 18일
Hello , how can I code this : y'(t)=cos(t + y) y(0)=0 t[0,3] (exact solution y(t)=-t + 2arctan(t)) using Euler's method like the code below?
clc
clear all
close all
h=0.01;
N=100; % number of steps
y(1)=1;
for n=1:N
y(n+1)=y(n)+h*(-6*y(n));
x(n+1)=n*h;
end
plot(x,y)
hold on
x=0:0.01:1;
y=exp(-6*x); % exact solution
plot(x,y)

댓글을 달려면 로그인하십시오.

### 답변(1개)

Paul 2021년 6월 18일
The code that is posted solves a differential equation of the form
y'(x) = f(x,y), y(0) = y0
and compares it to the exact solution
y(x) = some function of x
over the interval
x [0 1]
As it happens, your problem has the exact same form, with the independent variable being t instead of x and the endpoint of the interval is different.
##### 댓글 수: 2표시숨기기 이전 댓글 수: 1
Mike Asmanis 2021년 6월 18일
I'm getting Array indices must be positive integers or logical values.

댓글을 달려면 로그인하십시오.

### Community Treasure Hunt

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

Start Hunting!