Блок-схема вычисления функции

Uses GraphABC;
type mas= array [1..30,1..30] of real;
     mas1=array [1..30] of real;
var ii:text;
a: mas;
h,x1,y1,x,y,b: mas1;
vv,vz,k,i,j,m,n,l,r,max,mx,my,x0,y0,t: integer;
f,g,s,q1,z,p: real;
zt:array[-10..4000]of integer;
ft:array[-10..4000]of real;
v:text;
kl:file of integer;
 
procedure p1;
var f:real;
begin
max:=0; f:=a[k,k];
for j:=k+1 to n do
if f <abs(a[j,k]) then
begin  f:=abs(a[j,k]); max:=j;
end;
end;
 
procedure p2;
var l:integer;
begin
for l:=k to n do
begin  q1:=a[k,l]; a[k,l]:=a[max,l]; a[max,l]:=q1;
end;
q1:=b[k]; b[k]:=b[max]; b[max]:=q1;
end;
 
procedure p3;
var q1:real;
begin
q1:=a[k,k];
for j:=k to n do
a[k,j]:=a[k,j]/q1;
b[k]:=b[k]/q1;
for i:=k+1 to n do
begin   q1:=a[i,k];
for j:=k to n do a[i,j]:=a[i,j]-a[k,j]*q1;
b[i]:=b[i]-b[k]*q1;
end;
end;
 
procedure p4;
begin
for i:=n-1 downto 1 do
for j:=n downto i+1 do
b[i]:=b[i]-a[i,j]*b[j];
end;
 
 
procedure Aproksimatsiy;
begin
n:=4;
m:=7;
for i:=1 to n do
begin
for j:=1 to m do
begin
l:=5-i;
p:=1;
for t:=l to n-1 do
p:=p*x1[j];
y[i]:=y[i]+y1[j]*p;
end;
b[i]:=y[i];
end;
x[1]:=m;
for i:=2 to m do
begin
for j:=1 to m do
begin
l:=m+1-i;
p:=1;
for t:=l to m-1 do
p:=p*x1[j];
x[i]:=x[i]+p;
end;
end;
for i:=1 to n do
begin
for j:=1 to n do
begin
l:=i-1;
a[i,j]:=x[j+l];
end;
end;
 
{gaus}
for k:=1 to n do
begin P1;
if max<>0 then P2;
P3;
end;
P4;
f:=0;
n:=n-1;
{Рисование графика}
SetPenColor(rgb(random(200),random(200),random(200)));
zt[round(x1[1])-2]:=round(x1[1]-2);
for r:=round(x1[1]-1) to round(x1[m]+1) do
begin
f:=0;
zt[r]:=zt[r-1]+1;
for i:=1 to 4 do
begin
k:=i-1;
p:=1;
for j:=1 to k do
begin
p:=p*zt[r];
end;
f:=f+b[i]*p;
end;
ft[r]:=f;
end;
for r:=round(x1[1]+1) to round(x1[m]) do
begin
line(x0+(round(zt[r-1]))*mx,y0-(round(ft[r-1]))*mx,x0+(round(zt[r]))*mx,y0-(round(ft[r]))*my);
end;
writeln('Апроксимационный многочлен:');
{Вывод функции}
Write('y = ');
for i:=1 to n+1 do
begin
Write('(x^',n-i+1,')*',b[i]:1:5);
end;
end;
 
procedure Interpolytsiy;
begin
n:=7;
SetPenColor(rgb(random(255),random(255),random(255)));
for i:=1 to n-1 do
begin
h[i]:=(y1[i+1]-y1[i])/(x1[i+1]-x1[i]);
b[i]:=y1[i]-h[i]*x1[i];
end;
moveto(x0+round(x1[1])*mx,y0-round(y1[1])*my);
for i:=1 to n-1 do
begin
z:=x1[i];
while z<= x1[i+1] do
begin
f:=h[i]*z+b[i];
z:=z+0.01;
lineto(x0+round(z)*mx,y0-round(f)*my);
end;
end;
writeln;
writeln('Таблица значений:');
writeln;
writeln('   x        y');
assign(v,'c:\funktsiy.txt');
rewrite(v);
writeln(v,'Интерполяционный многочлен:');
for i:=1 to n-1 do
begin
writeln(v,'от x',i,'до x',i+1,'  y = ',h[i],'*x+(',b[i],')');
writeln(round(x1[i]):4,'   ',round(h[i]*(x1[i])+b[i]):7);
writeln(round(x1[i]+(x1[i+1]-x1[i])/2):4,'   ',round(h[i]*(x1[i]+(x1[i]-x1[i+1])/2)+b[i]):7);
end;
close(v);
end;
 
begin
 
 
 
{Масштабирование}
mx:=1;
my:=1;
Writeln;
 
 
{Центр осей}
x0:=120;
y0:=550*my;
 
 
{Длина оси ОХ}
t:=x0+530;
 
 
{Оси}
Line(0,y0,t,y0);
Line(x0,0,x0,(y0+200)*my);
line(t,y0,t-25,y0+5);
line(t,y0,(t-25),y0-5);
Line(x0,0,x0-5,25);
Line(x0,0,x0+5,25);
 
 
{Подписи осей}
textOut(x0+10,5,'Y');
textOut(t-10,y0+10,'X');
textOut(x0-10,y0+6,'0');
 
 
{Деления}
l:=-3;
while (l<11) do
begin
r:=l*50;
Line(r*mx+x0,y0+5,r*mx+x0,y0-5);
l:=l+1;
end;
l:=-20;
while (l<101) do
begin
r:=l*5;
Line(r*mx+x0,y0+3,r*mx+x0,y0-3);
l:=l+1;
end;
l:=-4;
while (l<11) do
begin
r:=l*50;
Line(x0-5,y0-r*my,x0+5,y0-r*my);
l:=l+1;
end;
l:=-31;
while (l<101) do
begin
r:=l*5;
Line(x0-3,y0-r*my,x0+3,y0-r*my);
l:=l+1;
end;
 
{подписи делений}
textOut(x0+50*mx-3,y0+6,'50');
textOut(x0+100*mx-3,y0+6,'100');
textOut(x0+150*mx-3,y0+6,'150');
textOut(x0+200*mx-3,y0+6,'200');
textOut(x0+250*mx-3,y0+6,'250');
textOut(x0+300*mx-3,y0+6,'300');
textOut(x0+350*mx-3,y0+6,'350');
textOut(x0+400*mx-3,y0+6,'400');
textOut(x0+450*mx-3,y0+6,'450');
textOut(x0+500*mx-3,y0+6,'500');
textOut(x0-50*mx-8,y0+6,'-50');
textOut(x0-100*mx-8,y0+6,'-100');
textOut(x0-20,y0-50*my-2,'50');
textOut(x0-26,y0-100*my-2,'100');
textOut(x0-26,y0-150*my-2,'150');
textOut(x0-26,y0-200*my-2,'200');
textOut(x0-26,y0-250*my-2,'250');
textOut(x0-26,y0-300*my-2,'300');
textOut(x0-26,y0-350*my-2,'350');
textOut(x0-26,y0-400*my-2,'400');
textOut(x0-26,y0-450*my-2,'450');
textOut(x0-26,y0-500*my-2,'500');
textOut(x0-26,y0+50*my-2,'-50');
textOut(x0-31,y0+100*my-2,'-100');
textOut(x0-31,y0+150*my-2,'-150');
 
assign(kl,'c:\znach.txt');
reset(kl);
for i:=1 to 7 do begin
read(kl,l);
x1[i]:=l;
read(kl,l);
y1[i]:=l;
end;
close(kl);
 
Interpolytsiy;
 
read(vz);
if (vz>=x1[1]) and (vz<=x1[7]) then begin
for i:=1 to 7 do
if vz>=x1[i]
then vv:=i;
{Рисование графика}
SetPenColor(rgb(random(255),random(255),random(255)));
circle(x0+(round(vz))*mx,y0-(round(h[vv]*vz+b[vv]))*mx,5);
circle(x0+(round(vz))*mx,y0-(round(h[vv]*vz+b[vv]))*mx,3);
SetPenColor(rgb(random(255),random(255),random(255)));
moveto(x0+round(x1[vv])*mx,y0-round(y1[vv])*my);
z:=x1[vv];
while z<= x1[vv+1] do
begin
f:=h[vv]*z+b[vv];
z:=z+0.01;
lineto(x0+round(z)*mx,y0-round(f)*my);
end;
end
else
Writeln('Точка не входит в заданный интервал');
 
Aproksimatsiy;
 
Writeln;
Writeln;
{Вывод  введённых точек}
SetPenColor(rgb(random(255),random(255),random(255)));
for i:=1 to 7 do
begin
circle(x0+round(x1[i]),y0-round(y1[i]),3 );
end;
end.