Вычислить В0 и В1 гиперболического приближения. Функция

Program Function_2;
Var x,y,bo,b1:Real;
i,N:integer;
Function sum(p,u:Real):real;
begin
p:=0;
For i:=0 to N do
begin
p:=p+1/u;
sum:=p;
end;
Function sum1(q,v:Real):real;
begin
q:=0;
For i:=0 to N do
begin
q:=q+v;
sum1:=q;
end;
Function sum2(q,v,u:Real):real;
begin
z:=0;
For i:=0 to N do
begin
z:=z+v/u;
sum2:=z;
end;
Function sumkv(c,u:Real):real;
begin
c:=0;
For i:=0 to N do
begin
c:=c+1/sqr(u);
sumkv:=c;
end;
BEGIN
Writeln('Vvedit x,y,N);
Readln(x,y,N);
b0:=sum(y,x)/sumkv(x);
Writeln('b0=',b0);
Readln(b0);
b1:=(sum1(y)-b0*N)/sum(x);
Writeln('b1=',b1);
Readln(b1);
end.

Program func;
function b0(x,y:real;m:integer):real;
Var
i:integer;
s:real;
begin
s:=0;
for i:=1 to m do
s:=s+x*y;
b0:=s;
end;
Function sumx(x:real;p:integer):real;
Var i:integer;
z:real;
begin
z:=0;
for i:=1 to p do
z:=z+x;
sumx:=z;
end;
Var x,y,bo,b1:real;
n:integer;
BEGIN
Write('Vvedite x,y,n');
Readln(x,y,n);
bo:=b0(x,y,n);
b1:=bo-bo*n*sumx(x,n);
Writeln('b0=',bo);
Writeln('b1=',b1);
END.