Решение дифференциальных уравнений параболического типа.
       3. Разностная схема Кранка-Николсона

     3.1. Вывод разностной схемы Кранка-Николсона

     Рассмотрим одномерное дифференциальное уравнение параболического типа с начальным и граничными условиями:
  (4.17)
Запишем для уравнения (4.17) разностную схему, представляя вторую производную функции u по координате х в виде двух слагаемых
 
и аппроксимируя первое слагаемое на n-ом шаге по времени, а второе - на (n + 1)-ом шаге по времени:
  (4.18)
Данная разностная схема называется разностной схемой Кранка-Николсона в честь авторов, создавших её.
     Разложение второй производной функции u по координате х на две равноценные составляющие, одна из которых аппроксимируется на n-ом шаге по времени, а другая - на (n + 1)-ом шаге по времени, указывает на то, что аппроксимацию этой производной в целом следует рассматривать на шаге по времени (n + 1/2). В то же время конечная разность, аппроксимирующая производную функции u по времени, по отношению к точке (n + 1/2) является центральной конечной разностью, имеющей, как известно, второй порядок аппроксимации. Следовательно, разностная схема Кранка-Николсона аппроксимирует уравнение (4.17) со вторым порядком и по времени, и по координате:
 
Таким образом, порядок аппроксимации разностной схемы Кранка-Николсона выше, чем порядок аппроксимации явной и неявной разностных схем; то есть, результаты, получаемые при использовании разностной схемы Кранка-Николсона будут более точными.