## Is there more than one Crank-Nicolson scheme?

Hi everybody...

I want to solve the diffusion equation in 1D using the Crank-Nicolson scheme. I have two books about numerical methods, and the problem is that in "Numerical Analysis" from Burden and Faires, the differences equation for the diffusion equations is:

$\frac{w_{i,j+1}-w_{i,j}}{k}-\frac{\alpha^2}{2h^2}\Big[w_{i+1,j}-2w_{i,j}+w_{i-1,j}+w_{i+1,j+1}-2w_{i,j+1}+w_{i-1,j+1}\Big]=0$

On the other hand, in "Numerical and analytical methods for scientists and engineers using mathematica", the same equation is expressed as:

$\frac{w_{i,j}-w_{i,j-1}}{k}-\frac{\alpha^2}{2h^2}\Big[w_{i+1,j}-2w_{i,j}+w_{i-1,j}+w_{i+1,j-1}-2w_{i,j-1}+w_{i-1,j-1}\Big]=0$

$i$ represents the space steps, $j$ the time steps, $k$ is $\Delta t$, $h$ is $\Delta x$

Should this schemes yield the same results? Why the differences?

I mean, in the first term of the first scheme, the numerator is $w_{i,j+1}-w_{i,j}$, but in the second scheme is $w_{i,j}-w_{i,j-1}$.

In addition to this, the last 3 terms of the equations (inside the brackets) are $w_{i+1,j+1}-2w_{i,j+1}+w_{i-1,j+1}$ and $w_{i+1,j-1}-2w_{i,j-1}+w_{i-1,j-1}$.

Are both schemes named Crank-Nicolson?

Can somebody help me with this?? Thanks!!

 PhysOrg.com science news on PhysOrg.com >> Hong Kong launches first electric taxis>> Morocco to harness the wind in energy hunt>> Galaxy's Ring of Fire
 Recognitions: Science Advisor It is just a difference in notation. If you replace j+1 by j and j by j-1 in the first equation, you get the second equation (but with the terms in the [ ] written in a different order). The method described in the first book is going to solve for the j+1 terms using the j terms. The second book is going to solve for the j terms using the j-1 terms.

 Tags crank nicolson, diffusion, heat, numerical, pde