Problem with Backward Euler Results
 

Hi,

I have tried to solve a system of equations in the form:
dy1/dt = y1 +dt*f(y1,y2)
dy2/dt = y2 +dt*f(y1,y2)

using different schemes such as forward Euler and backward (implicit) Euler and Runge Kutta order 4. My results seem to be OK for the forward Euler scheme and Runge Kutta above a certain number of time steps (A). But for implicit Euler using A time steps, I get an error after a certain number of time steps (say A-100). I was wondering if anyone had any ideas of what the problem might be. I can't think of anything at all.

Thanks!

[QUOTE=Hooman;336756]Hi,

I have tried to solve a system of equations in the form:
dy1/dt = y1 +dt*f(y1,y2)
dy2/dt = y2 +dt*f(y1,y2)

Do you perhaps mean f( ... ) for the first equation

and g( ... ) for the second equation?

yes, that's what I meant. Sorry!

Any ideas?

I just realized that if I use really small time steps (much smaller than that of explicit euler) my reuslts become OK. I am now thinking it could be because I used forward euler to work out the initial RHS's at the new time step, to calculate the y1 and y2 at the new time steps.

ooh I just realize that I have written the equations wrong here. They should be:
dy1/dt = f(y1,y2)
dy2/dt = g(y1,y2)