Euler Integration

  1. 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...


    Jimbo :)
  2. jcsd
  3. Tom Mattson

    Tom Mattson 5,526
    Staff Emeritus
    Science Advisor
    Gold Member

    Hi Jimbo, and welcome to PF.

    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:


  4. Thanks


    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 :)
  5. Tom Mattson

    Tom Mattson 5,526
    Staff Emeritus
    Science Advisor
    Gold Member

    Re: Thanks

    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!

    We're getting there, but I don't know what f(x,t) is. [?]

    edit: typo
  6. Integral

    Integral 7,288
    Staff Emeritus
    Science Advisor
    Gold Member

    Here is what I have on Euler's Method,

    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: May 22, 2003
  7. Integral

    Integral 7,288
    Staff Emeritus
    Science Advisor
    Gold Member

    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.
  8. 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.
  9. Thanks


    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 :)
Know someone interested in this topic? Share this thead via email, Google+, Twitter, or Facebook

Have something to add?