Physics Forums

Physics Forums (http://www.physicsforums.com/index.php)
-   General Math (http://www.physicsforums.com/forumdisplay.php?f=73)
-   -   Euler Integration (http://www.physicsforums.com/showthread.php?t=2365)

Jimbo May22-03 02:39 PM

Euler Integration
 
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 :)

Tom Mattson May22-03 04:51 PM

Re: Euler Integration
 
Hi Jimbo, and welcome to PF.

Quote:

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,

Jimbo May22-03 05:14 PM

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 :)

Tom Mattson May22-03 05:39 PM

Re: Thanks
 
Quote:

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!

Quote:

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

edit: typo

Integral May22-03 09:55 PM

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

Integral May23-03 01:40 AM

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.

gnome May23-03 07:37 AM

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.

Jimbo May23-03 04:31 PM

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 :)


All times are GMT -5. The time now is 05:26 PM.

Powered by vBulletin Copyright ©2000 - 2014, Jelsoft Enterprises Ltd.
© 2014 Physics Forums