为了用高斯公式和误差值计算积分，我在matlab中编写了这个程序：
%calculate integral with Gauss 3-point
clc
syms x ;
f=input('f(x) : ');
a=input('a = : ');
b=input('b = : ');
x=.5*((b-a)*(-1*(3/5)^0.5)+(a+b));
w0f0=(5/9)*eval(f);
x=.5*((b-a)*0+(a+b));
w1f1=(8/9)*eval(f);
x=.5*((b-a)*((.6)^.5)+(a+b));
w2f2=(5/9)*eval(f);
antegral=w0f0+w1f1+w2f2;
antegral=0.5*(b-a)*antegral;
disp(antegral)
%cal Error value = maximum of f^(6)/1570
syms x
m6=diff(f,6);
m6
m6=-1*m6;
[x,fval]=fminbnd(fun,x0,xn);
En=fval/15750;
disp (En);
此代码从命令行获取F（x）并计算积分。
我有2个问题：
如果输入f（x）= sin(x)，a = 0，b = pi在命令行中进行计算，可能会看到2.97708877614822e-009。我想成为数字格式
为某些E(n)计算f(x)，例如(1+x)^3 a = 0，b = 2有错误，必须计算0但计算{ {1}}