Classes of stable three-layer schemes 4470
Let H = SL h be the space of all complex functions defined on the grid
wh and vanishing at the points x = 0 and x = 1. Also, it will be sensible
to introduce
N-1
(y, v) = L Yk vk h,
k=l
where vk is the complex conjugate function of vk. The next step is to write
the initial equation in the formOU. 82 u.
- = -z --+ zqu
ot 8x^2
and introduce in the space H the operatorA y = i Yxx - i q 1l.
The operator A is skew-symmetric: (Ay, v) = -(y, Av), since
N-1 N-1
(Ay,v) =I: (iYxx,k)vkh+ I: (-iqyhvkh
k=l k=l
N-1 N-1
=-2.::: Yk(ivfa-,k)h-2.::: Yk(-iqv)kh=-(y,Av)
k=l k=lThe norm of the operator A admits the estin1ate4
llAll < h 2 +q.The explicit scheme (82) for the Schrodinger equation takes the fonn0
, yo t + i Yxx. - i q y = 0 1 y E Ah.It is stable for rll A II < 1, that is, under the constraints
orh2
T <
4 + q h^2As the second possible example we look at the schemeyo t +yo x = 0