Для заданной функции вывести на экран значения функции

program functia;
uses crt,graph;
function D(f:real):real;
begin
D:=0.6*sqrt(f*f-2*f+3)+3-1.1*sin(2*f/(f+1))/cos(2*f/(f+1))
end;
var a,c,f,x,h,mx,my:real;
    x0,y0,n,i:integer;
    s:string;
begin
clrscr;
a:=-0.8; {начало интервала}
repeat
write('Введите конец интервала -0.8<c<=2 c=');
readln(c); {конец интервала}
until(c>-0.8)and(c<=2);
h:=0.2;{шаг}
{вычисляем и выводим значения функции на интервале [-0.8;c] с шагом 0.2}
f:=a;
repeat
if cos(2*f/(f+1))=0 then writeln('f=',f:4:1,' функция не определена')
else writeln('f=',f:4:1,' D(f)=',D(f):7:3);
f:=f+h;
until f>c;
writeln;
write('Нажмите Enter для просмотра графика');
readln;
initgraph(x0,y0,''); {переходим в графический режим}
{определим положение начала координат на экране
и масштабы по осям для перевода реальных координат в экранные}
if c<=0 then {если все значения <=0}
 begin
  x0:=getmaxX-50;{ось Y справа}
  n:=trunc(abs(a*10));{удвоенное количество засечек по оси Х}
  mx:=(x0-40)/abs(a);{масштаб по Х}
 end
else {если c>0}
 begin
  x0:=round(-a*(getmaxX-30)/(c-a));{ось Y пропорционально отрицательной и
                                    положительной частям интервала}
  if abs(a)>abs(c) then{если отрицательная часть больше}
   begin
    n:=trunc(abs(a*10));
    mx:=(x0-10)/abs(a);
   end
  else {если положительная часть больше}
   begin
    n:=trunc(abs(c*10));
    mx:=(getmaxX-x0-40)/abs(c);
   end;
 end;
y0:=getmaxY div 2;{положение оси Х}
my:=(y0-20)/50;{масштаб по оси Y,
                достаточно поскольку значения функции не ограничены}
{проведем оси}
line(0,y0,getmaxX,y0);
line(x0,0,x0,getmaxY);
{нарисуем на них засечки и подпишем}
for i:=1 to n do
if i mod 2=0 then{шаг 0.2}
 begin
  line(x0+round(i*mx/10),y0-3,x0+round(i*mx/10),y0+3);
  str(i/10:0:1,s);
  outtextXY(x0+round(i*mx/10),y0+10,s);
  line(x0-round(i*mx/10),y0-3,x0-round(i*mx/10),y0+3);
  outtextXY(x0-round(i*mx/10),y0+10,'-'+s);
 end;
for i:=1 to 5 do
 begin
  line(x0-3,y0-round(i*my*10),x0+3,y0-round(i*my*10));
  str(i*10,s);
  if c<=0 then outtextXY(x0+10,y0-round(i*my*10),s)
  else outtextXY(x0-25,y0-round(i*my*10),s);
  line(x0-3,y0+round(i*my*10),x0+3,y0+round(i*my*10));
  if c<=0 then outtextXY(x0+10,y0+round(i*my*10),'-'+s)
  else outtextXY(x0-30,y0+round(i*my*10),'-'+s);
 end;
outtextXY(x0+5,y0+10,'0');
outtextXY(x0+5,10,'Y');
outtextXY(getmaxX-30,y0-20,'X');
{график точками, поскольку есть разрывы}
f:=a;
repeat
if cos(2*f/(1+f))<>0 then putpixel(x0+round(f*mx),y0-round(D(f)*my),12);
f:=f+0.0001;
until f>c;
readln
end.

program grafik;
uses crt,graph;
function Y(x:real):real;
begin
Y:=2.1*cos((x*x-4)/(x+3.4))+x/sqrt(abs(x))
end;
var a,b,x,h,maxy,mx,my:real;
    x0,y0,n,i:integer;
    s:string;
begin
clrscr;
a:=-3; {начало интервала}
repeat
write('Введите конец интервала -3<c<=6 b=');
readln(b); {конец интервала}
until(b>-3)and(b<=6);
h:=0.2;{шаг}
{вычисляем и выводим положительные значения функции на интервале [-3;c] с шагом 0.2}
x:=a;
n:=0;
maxy:=abs(Y(x));
repeat
if abs(Y(x))>maxY then maxY:=abs(Y(x));
if Y(x)>0 then
 begin
  writeln('x=',x:4:1,' Y(x)=',Y(x):7:3);
  n:=n+1;
 end;
if n mod 20=0 then
 begin
  write('Press Enter');
  readln;
 end;
x:=x+h;
until x>b;
writeln;
write('Нажмите Enter для просмотра графика');
readln;
initgraph(x0,y0,''); {переходим в графический режим}
{определим положение начала координат на экране
и масштабы по осям для перевоа реальных координат в экранные}
if b<=0 then {если все значения <=0}
 begin
  x0:=getmaxX-50;{ось Y справа}
  n:=trunc(abs(a*2));{удвоенное количество засечек по оси Х}
  mx:=(x0-40)/abs(a);{масштаб по Х}
 end
else {если c>0}
 begin
  x0:=round(-a*(getmaxX-30)/(b-a));{ось Y пропорционально отрицательной и
                                    положительной частям интервала}
  if abs(a)>abs(b) then{если отрицательная чать больше}
   begin
    n:=trunc(abs(a*2));
    mx:=(x0-10)/abs(a);
   end
  else {если положительная чать больше}
   begin
    n:=trunc(abs(b*2));
    mx:=(getmaxX-x0-40)/abs(b);
   end;
 end;
y0:=getmaxY div 2;{положение оси Х}
my:=(y0-40)/maxy;{масштаб по оси Y}
{проведем оси}
line(0,y0,getmaxX,y0);
line(x0,0,x0,getmaxY);
{нарисуем на них засечки и подпишем}
for i:=1 to n do
 begin
  line(x0+round(i*mx/2),y0-3,x0+round(i*mx/2),y0+3);
 // str(i/2:0:1,s);
  outtextXY(x0+round(i*mx/2),y0+10,s);
  line(x0-round(i*mx/2),y0-3,x0-round(i*mx/2),y0+3);
  outtextXY(x0-round(i*mx/2),y0+10,'-'+s);
 end;
for i:=1 to 2*(trunc(maxy)+1) do
 begin
  line(x0-3,y0-round(i*my/2),x0+3,y0-round(i*my/2));
 // str(i/2:0:1,s);
  outtextXY(x0-30,y0-round(i*my/2),s);
  line(x0-3,y0+round(i*my/2),x0+3,y0+round(i*my/2));
  outtextXY(x0-35,y0+round(i*my/2),'-'+s);
 end;
outtextXY(x0+5,y0+10,'0');
outtextXY(x0+10,10,'Y');
outtextXY(getmaxX-20,y0-20,'X');
{график}
setcolor(12);
x:=a;
repeat
if x=a then moveto(x0+round(x*mx),y0-round(Y(x)*my))
else lineto(x0+round(x*mx),y0-round(Y(x)*my));
x:=x+0.01;
until x>b;
readln
end.

uses crt;
function D(f:real):real;
begin
D:=0.6*sqrt(f*f-2*f+3)+3-1.1*sin(2*f/(f+1))/cos(2*f/(f+1))
end;
var a,c,f,x,h:real;
begin
clrscr;
a:=-0.8; {начало интервала}
repeat
write('Введите конец интервала -0.8<c<=2 c=');
readln(c); {конец интервала}
until(c>-0.8)and(c<=2);
h:=0.2;{шаг}
{вычисляем и выводим значения функции на интервале [-0.8;c] с шагом 0.2}
f:=a;
repeat
if cos(2*f/(f+1))=0 then writeln('f=',f:4:1,' функция не определена')
else writeln('f=',f:4:1,' D(f)=',D(f):7:3);
f:=f+h;
until f>c;
readln
end.

uses crt;
function Y(x:real):real;
begin
Y:=2.1*cos((x*x-4)/(x+3.4))+x/sqrt(abs(x))
end;
var a,b,x,h:real;
    n:integer;
begin
clrscr;
a:=-3; {начало интервала}
repeat
write('Введите конец интервала -3<c<=6 b=');
readln(b); {конец интервала}
until(b>-3)and(b<=6);
h:=0.2;{шаг}
{вычисляем и выводим положительные значения функции на интервале [-3;c] с шагом 0.2}
x:=a;
n:=0;
writeln('Положительные значения функции');
repeat
if Y(x)>0 then
 begin
  writeln('x=',x:4:1,' Y(x)=',Y(x):7:3);
  n:=n+1;
 end;
if n mod 20=0 then
 begin
  write('Press Enter');
  readln;
 end;
x:=x+h;
until x>b+h/2;
readln
end.