что делает эта программа распишите пожалуйста? это ведь метод прогонки?

#include <cstdio>
#include <math.h>
#include <iostream>
using namespace std;
 
int n,i;
double ax,bx,h,x,a[201],b[201],c[201],f[201],alf[201],bet[201],u[201];
void Show (int);
 
int _tmain(int argc, _TCHAR* argv[])
{
    setlocale (0,"");
    Show(n);
    ax=1; bx=3; n=20;
    h=(bx-ax)/n;
    x=1+h;
    for (i=1;i<n;i++)
    {
        a[i]=2+h/x;
        c[i]=2-h/x;
        b[i]=2*pow(h,2)/pow(x,2)-4;
        f[i]=16*pow(h,2)*x;
        x=x+h;
    }
    alf[1]=0; bet[1]=3;
    for (i=1;i<n;i++)
    {
        bet[i+1]=(f[i]-a[i]*bet[i])/(a[i]*alf[i]+b[i]);
        alf[i+1]=-c[i]/(a[i]*alf[i]+b[i]);
    }
    cin.get();
    u[n]=60.2958;
    for (i=n;i>0;i--)
    {
        u[i-1]=alf[i]*u[i]+bet[i];
        cout << "u[" << i << "]=" << u[i] << endl;
    }
    cout << "u[0]=" << u[0];
    cin.get();
}
 
void Show (int n)
{
    cout << "Исходная краевая задача:\n";
    cout << "        y'      y\n";
    cout << "y'' - ----- + ----- = 8x,   x=[1;3]\n";
    cout << "        x      x^2\n";
    cout << "y(1)=3\n";
    cout << "y(3)=60.2958\n";
    cout << endl;
    cout << "Искомая функция: y = x + xlnx + 2x^3\n";
    cout << endl;
}