Курсовая работа: Численное интегрирование методом Гаусса

{$N+}

{ Вычисление интегpала десятиточечным методом Гаусса }

uses crt,graph;

var aaa,bbb,kkk: real;

{константы десятиточечного метода Гаусса}

const

g10c1=0.9739065285/6. 2012983932;

g10c2=0.8650633667/6. 2012983932;

g10c3=0.6794095683/6. 2012983932;

g10c4=0.4333953941/6. 2012983932;

g10c5=0.1488743390/6. 2012983932;

g10x1=0.0666713443/6. 2012983932;

g10x2=0.1494513492/6. 2012983932;

g10x3=0.2190863625/6. 2012983932;

g10x4=0.2692667193/6. 2012983932;

g10x5=0.2955242247/6. 2012983932;

function F(x: real): real; {интегрируемая функция}

begin

F: = kkk*(exp(-aaa*x) - exp(-bbb*x));

end;

function gauss_calc(a,b: real): real; {сам десятиточечный метод Гаусса}

var n,m,s,s1,s2,s3,s4,s5: real;

begin

m: =(b+a) /2; n: =(b-a) /2;

s1: =g10c1*(f(m+n*g10x1) +f(m-n*g10x1));

s2: =g10c2*(f(m+n*g10x2) +f(m-n*g10x2));

s3: =g10c3*(f(m+n*g10x3) +f(m-n*g10x3));

s4: =g10c4*(f(m+n*g10x4) +f(m-n*g10x4));

s5: =g10c5*(f(m+n*g10x5) +f(m-n*g10x5));