Программа расчета для функции - Free Pascal

PROGRAM z1_3;
USES crt;
VAR a,b,x1,x2,y,dx:real;
n1:string;
f1:text;
BEGIN
clrscr;
write('Файл исходных данных: ');
readln(n1);
n1:=n1+'.txt';
assign(f1,n1);
reset(f1);
read(f1,a,b,x1,x2,dx);
close(f1);
TextBackGround(5);
Clrscr;
Textcolor(10);
Nosound;
Delay(2);
while x1<=x2 do
begin
if x1<=2.5 then y:=ln(x1*x1+b*b*b)/ln(10)-(a+b*b*x1)/PI
else if (x1>2.5) and (x1<5) then
y:=exp(a*x1)+(10*sqrt(5*exp(5*ln(x1))))
else if (x1>=5) then y:=ln(x1*x1)-arctan(1/exp(3*ln(a+x1)));
writeln('x=', x1:5:2,' y=',y:8:2);
x1:=x1+dx;
end;
readln
END.