7楼: Originally posted by 309065816 at 2013-07-11 09:22:05
=CombineTraprl('sqrt(4-x^2)',0,2,1.0e-4)
这是第一个程序里面最后的命令，而且按照编的程序可以看出，默认精度就是小数点后4位的，为什么出来的结果是这样的：
q =

   3.138365760415841


s =

    51...

9楼: Originally posted by 309065816 at 2013-07-11 09:23:43
第二个程序是要怎么写才能控制精度呢？
>> trapm(10000,0,2)

ans =

   3.141591477611321...

12楼: Originally posted by 月只蓝 at 2013-07-11 09:33:18
answer=trapm(10000,0,2)；
vpa(answer,5)    %显示5位有效数字...

11楼: Originally posted by 月只蓝 at 2013-07-11 09:30:31
“默认精度就是小数点后4位的”，这指的是计算精度，或者说实际上能达到的精度；
而q =3.138365760415841，这是显示精度，这很可能是format long命令引起的。
用下面一个例子来解释，计算精度和显示精度的区别： ...

14楼: Originally posted by 309065816 at 2013-07-11 09:40:05
嗯嗯，明白了。对了，2个程序运行出来的结果有点不一样，这个没多大影响吧？
answer=trapm(10000,0,2);
vpa(answer,5)    %显示5位有效数字

ans =

3.1416

>> =CombineTraprl('sqrt(4-x^2)',0,2 ...

15楼: Originally posted by 月只蓝 at 2013-07-11 09:59:55
两个程序结果不会一样的，本身两种程序涉及的算法就不一样。
被求函数存在原函数的，你可以把积分端点代入原函数，得出精确解，然后用两个程序的结果来比较。...