выч методы не разобраны! / Лаба№3 / Мишины / GAUSS6

{$N+,E+}
type vec=array[1..6]of extended;
var N:integer;A,B,S:extended;X,C:vec;
Function Fun(X:extended):extended;far;
Begin Fun:=sqrt(0.7*x*x+1)/(2.1*x+sqrt(0.6*x+1.7))End;
Procedure Gauss(N,K:integer;A,B:extended;var X,C:vec;var s:extended);
var i,j:integer;H,A1,B1,R,Q,S1:extended;
Begin H:=(B-A)/N;
A1:=A;
B1:=A+H;S:=0;
for i:=1to N do begin R:=(A1+B1)/2;
Q:=(B1-A1)/2;
S1:=0;
for j:=1to K do S1:=S1+C[J]*Fun(R+Q*X[j]);
S:=S+S1;
A1:=B1;B1:=A1+H
end;
S:=S*Q
End;
BEGIN X[1]:=0.932469;X[2]:=0.661209;X[3]:=0.238619;x[4]:=-X[1];X[5]:=-X[2];
X[6]:=-X[3];C[1]:=0.171324;C[2]:=0.360762;C[3]:=0.467914;C[4]:=C[1];
C[5]:=C[2];C[6]:=C[3];
repeat write('—Ёб«® Ё­вҐаў «®ў а §ЎЁҐ­Ёп Ё­вҐаў «  Ё­вҐЈаЁа®ў ­Ёп N,ЇаҐ¤Ґ«л Ё­вҐЈаЁа®ў ­Ёп A,B?');
readln(N,A,B);
Gauss(N,6,A,B,X,C,S);
writeln('€­вҐЈа «=',S)
until false
END.