
I chose a finitedifference method which is based
on the CranckNicholson implicit time differencing:

You can write the equation as follows :
where .

The term is
expanded as :
where can
be replaced by :

Finally if we employ forwardtime differences
and centeredspatial differences ,
the final equation is obtained :