y «+i y
hf + ^h2Df + ^h2(D2f + fyDf)+ (2)
Rozwinięcie prawej strony (1): korzystamy z ogólnego wzoru na szereg Taylora funkcji dwóch zmiennych
f^xn+aih,yn+^^i^hfn^ = ^Ui‘D‘f(xn,y„),
d_
dy
Stąd
=hf\xn+ ath,yn + £ =
¥
X, + <*ih>yn + Z [Piihfn Pij - )]J =
j=i
/=0 n 7=2 ^