newton_cotes.asv

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

function C=Newton_Cotes(n)
syms t
S1=1;
S2=1;
S3=1;
for k=0:n;
    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
    
end