yy积分

利⽤MATLAB编写⾼斯积分夏天四字词语

function[m]=intgauss(n)%n代表所求节点的个数

symsx

fori=1:n

forj=1:n

y(i,j)=int(log(x)*x^(i-1)*x^(n-j),0,1);%积分的权函数不同则log(x)会改变

end

end

y%y表⽰权函数与正交函数与w(x)的情侣qq名字 积分值所构成的⽅阵

fori=1:n

h(i)=-int(log(x)*x^(i-1)*x^n,0,1);

end

h=h'%h表⽰权函数与正交函数与w(x)的积分值所构成列向量

m=y雪豹悲歌读后感 ^(-1)*h;%m是系数向量

w=x^n;

fori=1:n

w=w+m(i)*x^(n-i);

end

x=solve(w==0)%x表⽰的是节点位置

fori=1:n

forj=1:n

l(i,j我的法兰西岁月 )=x(j)^(i-1);

end

z(i)=y(i家风议论文 ,n);

end

z=z';

A=l^(-红尘来呀来 1)

A=l^(-1)全国优秀人民警察 *z%A代表系数

fori=1:n

X=sprintf('%f*f(%f)',A(i),x(i));

disp(X)

end