yy積分

利?MATLAB編寫?斯積分夏天四字詞語

function[m]=intgauss(n)%n代表所求節(jié)點的個數(shù)

symsx

fori=1:n

forj=1:n

y(i,j)=int(log(x)*x^(i-1)*x^(n-j),0,1);%積分的權(quán)函數(shù)不同則log(x)會改變

end

end

y%y表?權(quán)函數(shù)與正交函數(shù)與w(x)的情侶qq名字 積分值所構(gòu)成的?陣

fori=1:n

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

end

h=h'%h表?權(quán)函數(shù)與正交函數(shù)與w(x)的積分值所構(gòu)成列向量

m=y雪豹悲歌讀后感 ^(-1)*h;%m是系數(shù)向量

w=x^n;

fori=1:n

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

end

x=solve(w==0)%x表?的是節(jié)點位置

fori=1:n

forj=1:n

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

end

z(i)=y(i家風(fēng)議論文 ,n);

end

z=z';

A=l^(-紅塵來呀來 1)

A=l^(-1)全國優(yōu)秀人民警察 *z%A代表系數(shù)

fori=1:n

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

disp(X)

end