Вычислительный эксперимент и методы вычислений. Ловецкий К.П - 13 стр.

UptoLike

ВУЗ: 

РУДН | Москва

Составители: 

Рубрика: 

Пример 2.
13
x
Рассмотрим функцию
sin 4y
=
на отрезке
[
]
0, / 2
π
и вы-
числим полиномы Лагранжа ().
n
Lx
Приведем текст подпрограмм-функций на паскале для вы-
числения значений многочлена Лагранжа для произвольно-
го значения и погрешности интерполяции на заданном
количестве точек.
n
Процедура вычисления полиномов Лагранжа ():
n
Lx
type
OneArgFunction = function (x:real) : real;
function L_n(x1:real; x : array of real; n:integer; f :
OneArgFunction):real;
var
s, sum : real;
i, j : integer;
begin
sum:=0;
for i:=1 to n do
begin
s:=1;
for j:=1 to n do
begin
if j<>i then
s:=s*(x1-X[j])/(X[i]-X[j])
end;
sum:=sum+s*f(X[i]);
end;
result := sum;
end;
Процедура оценки погрешности на отрезке
[
]
0, / 2
π
по
nPoints значениям:
function maxDeviation(xLeft, xRight : real; nPoints : integer) :
real;

Страницы