Euler Method without using ODE solvers

Question

Marko Kovacevic 2015-6-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/224319-euler-method-without-using-ode-solvers

评论： James Tursa 2016-4-11

I am trying to write a code that will solve a first order differential equation using Euler's method(Improved Euler's, Modified Euler's, and Euler-Cauchy). I don't want to use an ode solver, rather would like to use numerical methods which will return values for (x,y) and f(x,y) and plot of function f. I'm not sure how to begin to write this in MATLAB. I have solved the equation by hand and I'm now trying to write a code that solves that equation. Formulas for Improved Euler's, Modified Euler's, and Euler-Cauchy which I need to use are:

Improved Euler's:

Modified Euler's:

Euler-Cauchy's:

The equation to be used is y'=xy^2+y, y(0)=1, x∈[0,0.5], with step h=0.05

Any help is appreciated. Thank you in advance.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Marko Kovacevic 2015-6-22

Anyone?

请先登录，再进行评论。

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

Answer 1

James Tursa 2015-6-17

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/224319-euler-method-without-using-ode-solvers#answer_183060

编辑：James Tursa 2015-6-17

在 MATLAB Online 中打开

What have you done so far? Can you code up the simple Euler's Method? E.g., if you started with this generic expression of Euler's Method:

y(i+1) = y(i) + h * f(x(i),y(i))

can you put that inside of a MATLAB m-code for loop to generate results for some range of x with a given h?

Then, using that code, modify it for the other methods you are interested in.

Your "Modified Euler's" method doesn't look like Modified Euler to me ... it just looks like Euler's method using a h/2 step instead of a h step. (Maybe there is more to these methods in your doc that you didn't post?). E.g., see this link:

https://en.wikipedia.org/wiki/Heun%27s_method

6 个评论
显示 4更早的评论隐藏 4更早的评论

James Tursa 2015-6-17

编辑：James Tursa 2015-6-17

在 MATLAB Online 中打开

OK, here is a simple outline to get you started:

% Euler's Method
% Initial conditions and setup
h = 0.05;  % step size
x = 0:h:0.5;  % the range of x
y = zeros(size(x));  % allocate the result y
y(1) = 1;  % the initial y value
n = numel(y);  % the number of y values
% The loop to solve the DE
for i=1:n-1
    f = the expression for y' in your DE
    y(i+1) = y(i) + h * f;
end

You need to figure out how to write the code for f and insert it into the place I have indicated. Then you will have Euler's Method completed.

For your other methods, you will need to expand the code to calculate the multiple f values and y values and then combine them via the formula that is given.

(I used the variable i in the indexing above so that the code looks like the formulas in your text. Normally I would have used k for an index since MATLAB uses i for the imaginary value sqrt(-1), but I didn't want to confuse you with that since k has a specific meaning in some of the Runge Kutta links you may read)

Sanjida Ahmed 2016-4-11

Running this code,it gives nothing to display.whats the problem?I am new in Matlab.please help me.

James Tursa 2016-4-11

What code are you running? If you have a different problem to solve, it would be better to open up a new question.

请先登录，再进行评论。

Euler Method without using ODE solvers

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

回答（1 个）

6 个评论
显示 4更早的评论隐藏 4更早的评论

另请参阅

类别

标签

Community Treasure Hunt

Euler Method without using ODE solvers

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

回答（1 个）

6 个评论 显示 4更早的评论隐藏 4更早的评论

另请参阅

类别

标签

Community Treasure Hunt

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

6 个评论
显示 4更早的评论隐藏 4更早的评论