3_теориЛеготкина / выч методы / Лаба №1 / LOGR

uses graph;
const n=8;
var x,y:array[0..n]of real;
i,j,t,m,b:integer;
w:real;
function lagr(x1:real):real;
var i,j:integer;
ln,q,p:real;
begin
ln:=0;
for i:=0 to n do begin
q:=1;
p:=1;
for j:=0 to n do begin
if not (i=j) then p:=p*(x1-x[j]);
if not (i=j) then q:=q*(x[i]-x[j]);
end;
ln:=ln+y[i]*p/q;
end;
lagr:=ln;
end;

begin
for i:=0 to n do x[i]:=-4+i;
y[0]:=0;
y[1]:=2;
y[2]:=0;
y[3]:=-1;
y[4]:=0;
y[5]:=1;
y[6]:=0;
y[7]:=-2;
y[8]:=0;
m:=detect;
initgraph(m,b,'m:\bgi');
line(360,0,360,340);
line(0,170,720,170);
for i:=50 to 700 do begin
t:=round(-170*(lagr(6*i/360-6))/4+170);
putpixel(i,t,3);
end;
readln;
closegraph
end.