Ошибка (71 строка). Ожидалась функция, встречена процедура Invoke

program Gauss;
uses crt;
type vek1=array[1..2]of real; reafun1=function(x1:real);reafun2=function(x2:real);reafun3=function(x3:real);reafun4=function(x4:real);
var h,a,b,x,integral:real; i,n,m:byte;
function f1 (x1:real):real;
begin 
f1:=(x1*x1+1)/((x1*x1*x1+3*x1+1)*(x1*x1*x1+3*x1+1));
end;
function f2 (x2:real):real;
begin 
f2:=(arcsin(x2))/(sqrt(1+x2));
end;
function f3 (x3:real):real;
begin 
f3:=1/(x3+1);
end;
function f4 (x4:real):real;
begin 
f4:=sin(x4)/(5+3*sin(x4));
end;
 
procedure gaussint1(n1:byte;niz1,ver1:real;f1:reafun1;var int1:real);
var a1,x1:array[1..20]of real; j1:byte;
begin
if n1=2 then
begin x1[1]:=0.577350269189625764509148780502; a1[1]:=1.0; end;
if n1=4 then
begin
x1[1]:=0.861136311594052575223946488893;
x1[2]:=0.339981043584856264802665759103;
a1[1]:=0.347854845137453857373063949222;
a1[2]:=0.652145154862546142626936050778;
end;
if n1=6 then
begin
x1[1]:=0.932469514203152027812301554494;
x1[2]:=0.661209386466264513661399595020;
x1[3]:=0.238619186083196908630501721681;
a1[1]:=0.171324492379170345040296142173;
a1[2]:=0.360761573048138607569833513838;
a1[3]:=0.467913934572691047389870343990;
end;
if n1=8 then
begin
x1[1]:=0.960289856497536231683560868569;
x1[2]:=0.796666477413626739591553936476;
x1[3]:=0.525532409916328985817739049189;
x1[4]:=0.183434642495649804939476142360;
a1[1]:=0.101228536290376259152531354310;
a1[2]:=0.222381034453374470544355994426;
a1[3]:=0.313706645877887287337962201987;
a1[4]:=0.362683783378361982965150449277;
end;
if n1=10 then
begin
x1[1]:=0.973906528517171720077964012084;
x1[2]:=0.865063366688984510732096688423;
x1[3]:=0.679409568299024406234327365115;
x1[4]:=0.433395394129247190799265943166;
x1[5]:=0.148874338981631210884826001130;
a1[1]:=0.0666713443086881375935688098933;
a1[2]:=0.149451349150580593145776339658;
a1[3]:=0.219086362515982043995534934228;
a1[4]:=0.269266719309996355091226921569;
a1[5]:=0.295524224714752870173892994651;
end;
for j1:=1 to n1 do
begin x1[n1-(j1-1)]:=-x1[j1]; a1[n1-(j1-1)]:=a1[j1] end;
for j1:=1 to n1 do x1[j1]:=(ver1+niz1)/2+(ver1-niz1)/2*x1[j1];
int1:=0;
for j1:=1 to n1 do int1:=int1+a1[j1]*f1(x1[j1]);
int1:=(ver1-niz1)/2*int1;
begin
a1:=0;
b1:=sqrt(2)/2;
n1:=10;
m1:=100;
h1:=(b1-a1)/m1;
x1:=0;
for i1:=1 to m1 do begin gaussint(n1,a1,a1+h1,f1,integral1);
a1:=a1+h1; x1:=x1+integral1; end;
textcolor(11);
writeln('a1= ',a1,' b1= ',b1,' m1= ',m1,' x1= ',x1:0:6);
end;
end;
 
procedure gaussint2(n2:byte;niz2,ver2:real;f2:reafun2;var int2:real);
var a2,x2:array[1..20]of real; j2:byte;
begin
if n2=2 then
begin x2[1]:=0.577350269189625764509148780502; a2[1]:=1.0; end;
if n2=4 then
begin;
x2[1]:=0.861136311594052575223946488893;
x2[2]:=0.339981043584856264802665759103;
a2[1]:=0.347854845137453857373063949222;
a2[2]:=0.652145154862546142626936050778;
end;
if n2=6 then
begin
x2[1]:=0.932469514203152027812301554494;
x2[2]:=0.661209386466264513661399595020;
x2[3]:=0.238619186083196908630501721681;
a2[1]:=0.171324492379170345040296142173;
a2[2]:=0.360761573048138607569833513838;
a2[3]:=0.467913934572691047389870343990;
end;
if n2=8 then
begin
x2[1]:=0.960289856497536231683560868569;
x2[2]:=0.796666477413626739591553936476;
x2[3]:=0.525532409916328985817739049189;
x2[4]:=0.183434642495649804939476142360;
a2[1]:=0.101228536290376259152531354310;
a2[2]:=0.222381034453374470544355994426;
a2[3]:=0.313706645877887287337962201987;
a2[4]:=0.362683783378361982965150449277;
end;
if n2=10 then
begin
x2[1]:=0.973906528517171720077964012084;
x2[2]:=0.865063366688984510732096688423;
x2[3]:=0.679409568299024406234327365115;
x2[4]:=0.433395394129247190799265943166;
x2[5]:=0.148874338981631210884826001130;
a2[1]:=0.0666713443086881375935688098933;
a2[2]:=0.149451349150580593145776339658;
a2[3]:=0.219086362515982043995534934228;
a2[4]:=0.269266719309996355091226921569;
a2[5]:=0.295524224714752870173892994651;
end;
for j2:=1 to n2 do
begin x2[n2-(j2-1)]:=-x2[j2]; a2[n2-(j2-1)]:=a2[j2] end;
for j2:=1 to n2 do x2[j2]:=(ver2+niz2)/2+(ver2-niz2)/2*x2[j2];
int2:=0;
for j2:=1 to n2 do int2:=int2+a2[j2]*f2(x2[j2]);
int2:=(ver2-niz2)/2*int2;
begin
a2:=0;
b2:=sqrt(2)/2;
n2:=10;
m2:=100;
h2:=(b2-a2)/m2;
x2:=0;
for i2:=1 to m2 do begin gaussint2(n2,a2,a2+h2,f2,integral2);
a2:=a2+h2; x2:=x2+integral2; end;
textcolor(15);
writeln('a2= ',a2,' b2= ',b2,' m2= ',m2,' x2= ',x2:0:6);
end;
end;
 
procedure gaussint3(n3:byte;niz3,ver3:real;f3:reafun3;var int3:real);
var a3,x3:array[1..20]of real; j3:byte;
begin
if n3=2 then
begin x3[1]:=0.577350269189625764509148780502; a3[1]:=1.0; end;
if n3=4 then
begin
x3[1]:=0.861136311594052575223946488893;
x3[2]:=0.339981043584856264802665759103;
a3[1]:=0.347854845137453857373063949222;
a3[2]:=0.652145154862546142626936050778;
end;
if n3=6 then
begin
x3[1]:=0.932469514203152027812301554494;
x3[2]:=0.661209386466264513661399595020;
x3[3]:=0.238619186083196908630501721681;
a3[1]:=0.171324492379170345040296142173;
a3[2]:=0.360761573048138607569833513838;
a3[3]:=0.467913934572691047389870343990;
end;
if n3=8 then
begin
x3[1]:=0.960289856497536231683560868569;
x3[2]:=0.796666477413626739591553936476;
x3[3]:=0.525532409916328985817739049189;
x3[4]:=0.183434642495649804939476142360;
a3[1]:=0.101228536290376259152531354310;
a3[2]:=0.222381034453374470544355994426;
a3[3]:=0.313706645877887287337962201987;
a3[4]:=0.362683783378361982965150449277;
end;
if n3=10 then
begin
x3[1]:=0.973906528517171720077964012084;
x3[2]:=0.865063366688984510732096688423;
x3[3]:=0.679409568299024406234327365115;
x3[4]:=0.433395394129247190799265943166;
x3[5]:=0.148874338981631210884826001130;
a3[1]:=0.0666713443086881375935688098933;
a3[2]:=0.149451349150580593145776339658;
a3[3]:=0.219086362515982043995534934228;
a3[4]:=0.269266719309996355091226921569;
a3[5]:=0.295524224714752870173892994651;
end;
for j3:=1 to n3 do
begin x3[n3-(j3-1)]:=-x3[j3]; a3[n3-(j3-1)]:=a3[j3] end;
for j3:=1 to n3 do x3[j3]:=(ver3+niz3)/2+(ver3-niz3)/2*x3[j3];
int3:=0;
for j3:=1 to n3 do int3:=int3+a3[j3]*f3(x3[j3]);
int3:=(ver3-niz3)/2*int3;
begin
a3:=0;
b3:=sqrt(2)/2;
n3:=10;
m3:=100;
h3:=(b3-a3)/m3;
x3:=0;
for i3:=1 to m3 do begin gaussint3(n3,a3,a3+h3,podint3,integral3);
a3:=a3+h3; x3:=x3+integral3; end;
textcolor(13);
writeln('a3= ',a3,' b3= ',b3,' m3= ',m3,' x3= ',x3:0:6);
end;
end;
 
procedure gaussint4(n4:byte;niz4,ver4:real;f4:reafun4;var int4:real);
var a4,x4:array[1..20]of real; j4:byte;
begin
if n4=2 then
begin x4[1]:=0.577350269189625764509148780502; a4[1]:=1.0; end;
if n4=4 then
begin
x4[1]:=0.861136311594052575223946488893;
x4[2]:=0.339981043584856264802665759103;
a4[1]:=0.347854845137453857373063949222;
a4[2]:=0.652145154862546142626936050778;
end;
if n4=6 then
begin
x4[1]:=0.932469514203152027812301554494;
x4[2]:=0.661209386466264513661399595020;
x4[3]:=0.238619186083196908630501721681;
a4[1]:=0.171324492379170345040296142173;
a4[2]:=0.360761573048138607569833513838;
a4[3]:=0.467913934572691047389870343990;
end;
if n4=8 then
begin
x4[1]:=0.960289856497536231683560868569;
x4[2]:=0.796666477413626739591553936476;
x4[3]:=0.525532409916328985817739049189;
x4[4]:=0.183434642495649804939476142360;
a4[1]:=0.101228536290376259152531354310;
a4[2]:=0.222381034453374470544355994426;
a4[3]:=0.313706645877887287337962201987;
a4[4]:=0.362683783378361982965150449277;
end;
if n4=10 then
begin
x4[1]:=0.973906528517171720077964012084;
x4[2]:=0.865063366688984510732096688423;
x4[3]:=0.679409568299024406234327365115;
x4[4]:=0.433395394129247190799265943166;
x4[5]:=0.148874338981631210884826001130;
a4[1]:=0.0666713443086881375935688098933;
a4[2]:=0.149451349150580593145776339658;
a4[3]:=0.219086362515982043995534934228;
a4[4]:=0.269266719309996355091226921569;
a4[5]:=0.295524224714752870173892994651;
end;
for j4:=1 to n4 do
begin x4[n4-(j4-1)]:=-x4[j4]; a4[n4-(j4-1)]:=a4[j4] end;
for j4:=1 to n4 do x4[j4]:=(ver4+niz4)/2+(ver4-niz4)/2*x4[j4];
int4:=0;
for j4:=1 to n4 do int4:=int4+a4[j4]*f4(x4[j4]);
int4:=(ver4-niz4)/2*int4;
begin
a4:=0;
b4:=sqrt(2)/2;
n4:=10;
m4:=100;
h4:=(b4-a4)/m4;
x4:=0;
for i4:=1 to m4 do begin gaussint4(n4,a4,a4+h4,f4,integral4);
a4:=a4+h4; x4:=x4+integral4; end;
textcolor(12);
writeln('a4= ',a4,' b4= ',b4,' m4= ',m4,' x4= ',x4:0:6);
end;
end;
 
begin
gaussint1(1,1,1,1,1,1,1);
gaussint2(1,1,1,1,1,1,2);
gaussint3(1,1,1,1,1,1,3);
gaussint4(1,1,1,1,1,1,4);
end.