Краевая задача разностным методом

program odu2test;
{Решение тестовой задачи с известным точным решением yt(t)}
const n=10; {Число точек}
      a=0;  {Границы интервала}
      b=1;
      {Краевые условия}
      k1=1; k2=0; a1=1;
      l1=1; l2=0; b1=10;
type vec=array [0..n] of real;
var
 h,eps,epsmax:real;
 al,bet,t,u,aa,bb,cc,ff:vec;
 i,j:integer;
 
{Функция коэффициентов при 1-й производной u'(t)}
function p(t:real):real;
begin
 p:=sin(t)/cos(t);
end;
{Функция коэффициентов при u(t)}
function g(t:real):real;
begin
 g:=sqr(cos(t));
end;
{Функция правой части уравнения}
function f(t:real):real;
begin
 f:=0;
end;
{Точное решение yt(t)}
function yt(t:real):real;
var c1:real;
begin
 c1:=(10-cos(sin(1)))/sin(sin(1));
 yt:=cos(sin(t))+c1*sin(sin(t));
end;
 
begin
 {Вычисление шага сетки}
 h:=(b-a)/n;
 {Формирование векторов коэффициентов для 3-диагональной СЛАУ}
 for i:=0 to n do begin
  t[i]:=a+h*i;
  aa[i]:=1-p(t[i])*h/2;
  bb[i]:=1+p(t[i])*h/2;
  cc[i]:=2-g(t[i])*h*h;
  ff[i]:=h*h*f(t[i]);
 end;
 {Краевое условие при t=a}
 al[1]:=k2/(k2-k1*h);
 bet[1]:=-(a1*h)/(k2-k1*h);
 {Прямой ход метода прогонки}
 for i:=1 to n-1 do begin
  al[i+1]:=bb[i]/(CC[i]-al[i]*aa[i]);
  bet[i+1]:=(aa[i]*bet[i]-ff[i])/(cc[i]-al[i]*aa[i]);
 end;
 {Решение на правой границе отрезка}
 u[n]:=(l2*bet[n]+b1*h)/(l2+h*l1-l2*al[n]);
 {Обратный ход метода прогонки}
 for i:=n-1 downto 0 do
  u[i]:=al[i+1]*u[i+1]+bet[i+1];
 {Вывод таблицы и определение макс.погрешности}
 writeln('i':4,'t[i]':10,'y[i]':10,'yt(t[i])':10,'Погрешн.':10);
 for i:=0 to n do begin
  eps:=abs(u[i]-yt(t[i]));
  if eps>epsmax then epsmax:=eps;
  writeln (i:4,t[i]:10:4,u[i]:10:4,yt(t[i]):10:4,eps);
 end;
 writeln ('Максимальная погрешность=',epsmax);
 reset (input); readln;
end.