newton_cotes.m

为配套王能超版本数值分析简明教程编写,方便学生教师演示!

function C=Newton_Cotes(n)
syms t
for k=0:n;
    S1=1;
    S2=1;
    S3=1;
    for i=k:-1:1
       S2=S2*i;
    end
    for i=n-k:-1:1
       S3=S3*i;
    end
    for i=0:n
        if i~=k
            S1=(t-i)*S1;
        end
    end
    C(k+1)=(-1)^(n-k)/(n*S2*S3)*int(S1,t,0,n);
end