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