|
Здравствуйте!Требуется приближенно посчитать интеграл методом Гаусса для 3 точек,делением отрезка интегрирования на r частей. Ниженаписанное не работает.
f(x)=sin(x^2);
b:=-1:
r:=3:
> MetGauss:=0:
> for i from 0 to r while b<=1 do
> b:=-1+2/r*(i+1):
> a:=-1+2/r*i:
> eqs1:={c1+c2+c3=b-a,c1*x1+c2*x2+c3*x3=1/2*(b^2-a^2),c1*x1^2+c2*x2^2+c3*x3^2=1/3*(b^3-a^3),c1*x1^3+c2*x2^3+c3*x3^3=1/4*(b^4-a^4),c1*x1^4+c2*x2^4+c3*x3^4=1/5*(b^5-a^5),c1*x1^5+c2*x2^5+c3*x3^5=1/6*(b^6-a^6)}:
> l:=solve(eqs1,[x1,x2,x3,c1,c2,c3]);
> MetGauss:=evalf[13](MetGauss+l[1][4]*f(l[1][1])+l[1][5]*f(l[1][2])+l[1][6]*f(l[1][3])):
> od;
Решая систему eqs1 ,получаем несколько различных решений,которые образуют список списков вида [[x1=bla,x2=bla..],[x1=bla1...]];как изменить этот список списков,чтобы в каждом его списке у каждого элемента отсутствовали его первые два символа?
символы х1,х2 и прочие мешают,значение функции от элементов списка списков уже не найдешь..
|
