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

if x1<0.5 then x1: =-0.5;

y2: =f(ln(bbb/aaa) /(bbb-aaa)) *1.2;

y1: =-y2;

{Линия y=0}

setcolor(15);

y: =0;

j: =trunc(480*(y-y2) /(y1-y2));

line(0,j,639,j);

{Линии x=a,x=b}

setcolor(5);

j: =trunc(480*(-y2) /(y1-y2));

i: =trunc(640*(a-x1) /(x2-x1));

j1: =trunc(480*(f(a) - y2) /(y1-y2));

line(i,j, i,j1);

i: =trunc(640*(b-x1) /(x2-x1));

j: =trunc(480*(-y2) /(y1-y2));

j1: =trunc(480*(f(b) - y2) /(y1-y2));

line(i,j, i,j1);

{Сам график}

setcolor(14);

setlinestyle(0,0,3);

t: =(b-a) /n;

k: =0;

j1: =trunc(480*(-y2) /(y1-y2));

for i: =0 to 640 do begin

x: =(x2-x1) *i/640+x1;

y: =f(x);

j: =trunc(480*(y-y2) /(y1-y2));

if j>479 then j: =479;