Скачиваний:
17
Добавлен:
06.02.2016
Размер:
833 б
Скачать
Program L;
const
n=4;
var
xx,f:array[0..n] of real;
j:integer;
ff:real;
function lag (x:real):real;
var k,i:integer;
p1,p2,sum:real;
begin
sum:=0;
for k:=0 to n do
begin
p1:=1;
p2:=1;
for i:=0 to n do
begin
if i<>k then
begin
p1:=p1*(x-xx[i]);
p2:=p2*(xx[k]-xx[i]);
end;
end;
sum:=sum+p1/p2*f[k] ;
end;
lag:=sum;
end;
begin
for j:=0 to n do
begin
// xx[j]:=2+4/n*j; // равномерная сетка - 0 и 2 отрезок
xx[j]:=3+cos((2*j+1)*Pi/(2*(n+1))); //чебышевская сетка - коэфициенты определяются исходя из промежутка
f[j]:=sqrt(Abs(xx[j]-Pi)); // своя формула
end;
for j:=0 to 100 do
begin
ff:=lag(2+4/100*j);
// writeln(0+2/100*j); // значения
writeln(ff); // ответы
end;
end.

Соседние файлы в папке Новая папка