The su1n1narized approximation method 631grid in T on the segment 0 < t < T and a grid wh in G remains unchanged
(for more detail see Chapter 5).
If G is a parallelepiped in the space RP, problem (63)-(65) can be
solved through the use of an economical factorized scheme with accuracy
O(r^2 + lhl^2 ). The design of such a scheme was made in Section 2 and it
was investigated there in full details. Applying the same procedure serves
ti 1not.ivate that, first, the operators(66)with fa still subject the normalization conditionare approximated successively with step T /p.
Second, the difference expresions(67)are ai111ed at approxin1ating to the derivative 82 tt/8t^2 with step r/p, where
u(u)=tt - j-l+a/2 , ll(o)-tt-tl - j-1 ,
( 68)u(a) - ll(c"-1) - tl(a-2) + il(o:)
rv - --
T2 9 fJt2 'C\'=l,2,3 for p=:3,
with the members u(-l) = ii( 2 ) = u(j-l)+^2!^3 and u(- 2 ) = u(l) = ttj-^2 /^3.
Third, the operator Lau+ fa is approximated to second order by the
homogeneous difference operator Ao:Y +'Po: on the grid wh in the space RP.
The coefficients at the member Ao: and the righ-hand side 'Po: are taken at
the moment