Euler Integration

  • Thread starter Jimbo
  • Start date
  • #1
10
0

Main Question or Discussion Point

Hi

Im a bit stuck on the method for Euler Integration. I have the following first order differential equation:

dx/dt = (x-at) / (x / a+t)

where constant a = 1.0V/s, and initial condition x = 1.0V at t=0s

I have a time step of 0.02 and I need to calculate the output voltage at a time t = 0.2s

So I have the Euler Integration formula: x(t + ^t) = x(t) + ^t.f(x,t)
and I've started putting the values into an Excel spreadsheet, with columns as follows:

t | x(t) | f(x,t) | ^x = ^t.f(x,t)

I have values for all the first row for these columns but an example in a textbook has an additional column called x(exact) which appears to calculate an error value? I am not sure how to calculate this value, and am also unsure about how to move on with the iteration? What values should I move down for further rows?

I have searched the web but can't find any sites I can make sense of - can anybody help me? I hope my question makes sense...

Thanks

Jimbo :)
 

Answers and Replies

  • #2
Tom Mattson
Staff Emeritus
Science Advisor
Gold Member
5,500
7
Hi Jimbo, and welcome to PF.

Originally posted by Jimbo
So I have the Euler Integration formula: x(t + ^t) = x(t) + ^t.f(x,t)
and I've started putting the values into an Excel spreadsheet, with columns as follows:
This bit I don't understand. Could you try to clear that up? If you are trying to do superscripts, do this:

x[ sup ]2[ /sup ]

with no spaces in the brackets. The result will be:

x2

Thanks,
 
  • #3
10
0
Thanks

Hi

Thanks for your reply

Sorry about the unclear formula - the ^ were meant to be deltas to signify the small change in t

x(t + delta(t)) = x(t) + delta(t).f(x,t)

Hope that makes more sense now

Jimbo :)
 
  • #4
Tom Mattson
Staff Emeritus
Science Advisor
Gold Member
5,500
7


Originally posted by Jimbo
x(t + delta(t)) = x(t) + delta(t).f(x,t)
OK, that helps. If you want to make it look prettier, check the Announcement at the top of the forum Howto: Making Math Symbols.

& Delta ;

without the spaces gives you:

Δ

Oooohhhh! Pretty!

Hope that makes more sense now
We're getting there, but I don't know what f(x,t) is. [?]

edit: typo
 
  • #5
Integral
Staff Emeritus
Science Advisor
Gold Member
7,198
55
Here is what I have on Euler's Method,

Given
y' = f(x,y)
y(a)= y0

For a solution over the interval [a,b]
Choose a step h=(b-a)/N. Set xn= a + nh for n= 0,1,2...N

yn+1 = yn+ hf(xn,yn)

Error = h2y''(ξ)/2

Where xn < &xi; < xn+ h

I am haveing trouble intrepeting your notation, perhaps you can apply my notation to your problem.

Edit: Typo+ some content
 
Last edited:
  • #6
Integral
Staff Emeritus
Science Advisor
Gold Member
7,198
55
I just reread your initial post, you did say EXACT and not error term. The ONLY way to generate the exact solution is to solve the DE. This solution is what Eulers, or any other numerical method, is approximating.
 
  • #7
1,036
1
My textbook has an "exact" column too. That's not part of the Euler's method computation. It's there just for purposes of the example, to show you the amount of the error in the Euler approximation.

As Integral said, that value was obtained by solving the DE (i.e. by integration), then evaluating the solution for various values of t so that the error could be determined for each iteration.
 
  • #8
10
0
Thanks

Hi

Thanks very much everyone for your help

I have it sorted now. gnome was correct about the exact column not being part of the calculation.

To move onto the next iteration I needed to add the value &Delta;(t).f(x,t) to x(t) to get a new value for x(t)

Thanks again

Jimbo :)
 

Related Threads on Euler Integration

Replies
1
Views
519
Replies
3
Views
4K
Replies
2
Views
2K
  • Last Post
Replies
5
Views
743
Replies
36
Views
15K
  • Last Post
Replies
3
Views
881
  • Last Post
Replies
5
Views
601
  • Last Post
Replies
11
Views
2K
  • Last Post
Replies
5
Views
1K
Top