matlab 精度不够,matlab中Romberg的精度为什么不够？帮帮忙

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

function A=romberg(s,a1,a2,eps)%参数依次为要积分的函数，下界，上届，精度。syms xs=(x^3+sin(x))/x;%('输入被积函数的表达式:')a1=0.3;%('输入积分下限:')a2=0.8;%('输入积分上限:')eps=0.00000001;%('输入精度:')h=a2-a1;A(1,1)=(h/2)*(subs(s,x,a1)+subs(s,x,a2));A(2,1)=A(1,1)/2+(h/2)*subs(s,x,a1+(h/2));i=2;j=1;n=2;while 1   m=1;   h=h/2;       while j

0.6529 0 0 0

0.6397 0.6353 0 0

0.6364 0.6353 0.6353 0

0.6355 0.6353 0.6353 0.6353