Euler Integration

  • Thread starter Jimbo
  • Start date
  • #1
Jimbo
10
0
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
quantumdude
Staff Emeritus
Science Advisor
Gold Member
5,566
23
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
Jimbo
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
quantumdude
Staff Emeritus
Science Advisor
Gold Member
5,566
23


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,224
63
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,224
63
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
gnome
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
Jimbo
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 :)
 

Suggested for: Euler Integration

  • Last Post
Replies
1
Views
758
Replies
3
Views
4K
Replies
9
Views
2K
Replies
2
Views
2K
  • Last Post
Replies
5
Views
910
  • Last Post
Replies
5
Views
4K
  • Last Post
Replies
5
Views
2K
Replies
36
Views
16K
  • Last Post
Replies
5
Views
1K
  • Last Post
Replies
12
Views
4K
Top