Проверить программу нахождения корня уравнения

const a=0;
      b=0.8;  // начало и конец промежутка
var   e, x2, x1, c :real; // точность, корень уравнения
 
function tg(x:real): real;
begin
 result:= sin(x)/cos(x);
end;
 
function sec(x:double): double;
begin
 result:= 1/cos(x);
end;
 
 //Значение функции
 function f(x:real):real; //  функция
   begin
    f:=sin(x)/cos(x)+1/3*power(sin(x)/cos(x),3)+1/5*power(sin(x)/cos(x),5)-1/3;
   end;
 
 //Производная функции
 // первая производная функции
 function f1(x:real):real;
    begin
     f1:=1/8*(cos(4*x)+7)*power(1/cos(x),6);
    end;
 
 // вторая производная функции
  function f2(x:real):real;
     begin
      f2:=3/4*(cos(4*x)+7)*sin(x)/cos(x)*power(1/cos(x),6)-1/2*sin(4*x)*power(1/cos(x),6);
     end;
     begin
      if(f(a)*f1(a)>0)
then x1:=a
else x1:=b;
 begin
 
  if(f(a)*f2(a)>0)
  then c:=b
  else c:=a;
 end;
  begin
   write ('Точность=');
   read(e);
   repeat
    x1:=x1-f(x1)*(x1-c)/(f(x1)-f(c));
    c:=x2;
    x2:=x1;
   until abs(c-x1)<e;
  writeln(' ',x1);
   readln;
end;
end.

var a,b,e,x,c,x1:real;
Function f(x:real): real;
 
begin
f:=tan(x)+1/3*tan(x)*tan(x)*tan(x)+1/5*tan(x)*tan(x)*tan(x)*tan(x)*tan(x)-1/3
End;
 
Function f1(x:real): real;
Begin
 f1:=1/8*(Cos(4*x)+7)*sqr(1/cos(x))*sqr(1/cos(x))*sqr(1/cos(x));
End;
 
Function f2(x:real): real;
Begin
 f2:=3/4*(Cos(4*x)+7)*tan(x)*sqr(1/cos(x))*sqr(1/cos(x))*sqr(1/cos(x))-1/2*sin(4*x)*sqr(1/cos(x))*sqr(1/cos(x))*sqr(1/cos(x));
End;
 
 
Begin
  read(a,b,e);
if(f(a)*f2(a)>0)
  then  x:=b
  else  x:=a;
 
if(f(a)*f2(a)>0)
 then   c:=a
 else   c:=b;
 
repeat
   x:=x-f(x)*(x-c)/(f(x)-f(c));
   c:=x1;
   x1:=x;
 
until abs (f(x))<e;
writeln(x);
end.