Здравствуйте! У меня возникла проблема. Необходимо построить графики функции и многочлена на одной координатной плоскости. Но у меня, почему-то, график функции (f(t)) строится, а многочлена (L2[c](t)) - нет. Помогите, пожалуйста.
Вот мой код:
Код:
>restart;
f := x -> 1/(1 + c*x^2);
N := 4;
a := -1.0;
b := 1.0;
h := (b - a)/N;
>c:=2;#c:=20;
x[0] := a;
y[0] := f(x[0]);
for i from 1 to N do
x[i] := x[i-1] + h;
y[i] := f(x[i]);
end do:;
>for N from 2 to 4 do
for i from 0 to N do
r:=1;for k from 0 to N do
if k <> i then r:=r*(t - x[k])/(x[i] - x[k]);
end if;
end do;
L[i] := r;
print (L[i]);
end do;
> L1[n] := 0;
for p from 0 to N do
L1[n] := L1[n] + f(x[p])*L[p];
end do;
print (L1[n]);
L2[n] := t -> L1[n];
print (L2[n]);
simplify (L1[n]);
print (____________________________________________________________________________________________________________________);
end do;
> for c from 2 to 4 do
with(plots):
Q1 := plot (f(t), t = -1..1, colour=green, thickness=2):
Q2 := plot (L2[c](t), t = -1..1, colour=red, thickness=2):
display(Q1, Q2);
end do;