Фунция с построеним графика!!!

program test;
var m,n,i,j,p:integer;b,max,min,xj,s,ha,y,an:real;
 
function a(j:integer;an,ha:real):real;
begin
a:=an+j*ha;
end;
 
function c(i:integer;b:real):real;
begin
c:=exp(ln(3+b)*(1/i));
end;
 
function x(j:integer;b,an,ha:real):real;
begin
x:=(a(j,an,ha)*j+1)/(b*j-3);
end;
 
begin
write('vvedite m ');readln(m);
write('vvedite n ');readln(n);
write('vvedite an ');readln(an);
write('vvedite ha ');readln(ha);
write('vvedite b ');readln(b);
max:=0;min:=0;p:=0;
for j:=1 to m do
if j mod 4=0 then
begin
xj:=x(j,b,an,ha);
if xj<=b then
y:=((5*(xj*xj*xj))-6*sqr(xj)+j)/(sqr(j)+xj-2)+2*sqrt(abs(xj))
else
begin
s:=0;
for i:=1 to n do
s:=s+(sqr(c(i,b))+exp(ln(3)*(a(j,an,ha)/b))+2*b);
y:=s;
end;
if p=0 then begin max:=y; min:=y; p:=1;end;
if max<y then max:=y;
if min>y then min:=y;
end;
writeln('max ',max:10:2);
writeln('min ',min:6:2);
readln;
end.

uses crt,graph;
const m=29;
      n=5;
var j,i,jmn,jmx,gd,gm,x0,y0,xp,yv,yn:integer;
    an,h,a,b,s,c,sg,mn,mx,ym,xm:real;
    x,y:array[1..m] of real;
    f:text;
    st:string;
begin
clrscr;
an:=7.3;
h:=3.7;
b:=-1.5;
assign(f,'C:\Rezult\file1.txt');
rewrite(f);
writeln(f,'Znacenia argumenta X:');
for j:=1 to m do{массив аргументов}
 begin
  a:=an+j*h;
  x[j]:=(a*j+1)/(b*j-3);;
  write(f,x[j]:8:2);
 end;
writeln(f,'');
writeln(f,'');
writeln(f,'Znachenia functii Y:');
for j:=1 to m do{массив значений функции}
 begin
  if x[j]<=b then
  y[j]:=(5*x[j]*x[j]*x[j]-6*x[j]*x[j]+j)/(j*j+x[j]-2)+2*sqrt(abs(x[j]))
  else
   begin
    s:=0;a:=an+j*h;
    for i:=1 to n do{вычисляем сумму}
     begin
      c:=exp(ln(3+b)/i);
      s:=s+c;
     end;
    y[j]:=s+exp(ln(3)*(a/b))-5*a+2*b;
   end;
  write(f,y[j]:8:2);
 end;
writeln(f,'');
writeln(f,'');
{находим макс и мин для о кратных 4}
writeln('Максимальные и минимальные значения для j, кратных 4:');
writeln(f,'Maximalnye i mimimalnye znachenia dlya J, kratnyh 4:');
mx:=y[4];mn:=y[4];jmx:=4;jmn:=4;
for j:=1 to m do
if j mod 4=0 then
 begin
  if y[j]>mx then
   begin
    mx:=y[j];
    jmx:=j;
   end;
  if y[j]<mn then
   begin
    mn:=y[j];
    jmn:=j;
   end;
 end;
writeln('Максимальное=',mx:8:2,'  j=',jmx);
writeln(f,'Maximalnoe=',mx:8:2,'  j=',jmx);
writeln('Минимальное =',mn:8:2,'  j=',jmn);
writeln(f,'Minimalnoe =',mn:8:2,'  j=',jmn);
write('Fail zapisan!');
close(f);
{найдем мин. и макс.функции для графика}
mn:=y[1];mx:=y[1];
for i:=1 to m do
 begin
  if y[i]<mn then mn:=y[i];
  if y[i]>mx then mx:=y[i];
 end;
writeln;
writeln;
writeln('Max i min dlya postroenia grafika:');
writeln('Min=',mn:0:2);
writeln('Max=',mx:0:2);
write('Press Enter...');
readln;
gd:=0;
initgraph(gd,gm,'');
x0:=20;y0:=300;{начало координат}
xm:=600/28;ym:=400/((mx-mn));{коэффициенты масштабирования по осям}
xp:=getmaxX-x0;{парвая граница оси Х}
yv:=y0-round(mx*ym);{верх оси У}
yn:=y0-round(mn*ym)+30;{ее низ}
{ось Х}
line(x0,y0,xp,y0);
outtextXY(xp+5,y0-15,'j');
for j:=0 to m-1 do
 begin
  line(x0+round(j*xm),y0+3,x0+round(j*xm),y0-3);{засечки на оси}
  str(j+1,st);{в строку номера точек массива}
  outtextXY(x0+round(j*xm),y0+10,st);{выводим шкалу на ось}
 end;
{ось Y}
line(x0,yn,x0,yv);
outtextXY(x0-15,yv-15,'F(j)');
{вверх}
for j:=1 to round(mx/1000) do{значения уменьшаем в 1000 раз}
 begin
  line(x0-3,y0-round(j*ym*1000),x0+3,y0-round(j*ym*1000));{засечки}
  str(j,st);
  outtextXY(x0+5,y0-round(j*ym*1000)-3,st+'000');{подписи+000}
 end;
{вниз}
for j:=1 to abs(round(mn/1000))+1 do
 begin
  line(x0-3,y0+round(j*ym*1000),x0+3,y0+round(j*ym*1000));
  str(j,st);
  outtextXY(x0+8,y0+round(j*ym*1000)-3,'-'+st+'000');
 end;
{график}
moveto(x0,y0);
for j:=1 to m-1 do
 begin
  Setcolor(4);
  lineto(x0+round(j*xm),y0-round(y[j]*ym));{линия графика}
  Setcolor(7);
  circle(x0+round(j*xm),y0-round(y[j]*ym),2);{маркеры}
 end;
readln
end.