悠若谷 发表于 2009-6-11 17:22

对erfc函数的平方怎么求定积分?

我知道对erfc(x)函数作定积分可以这样写:
ierfcy=int('erfc(x)',0,0.1) ; interfcx=subs(ierfcy);
>> interfcx
interfcx =
   0.09436748855146

那么对erfc(x)^2求定积分呢?
ierfcy=int('erfc(x)^2',0,0.1) ; interfcy=subs(ierfcy);
Warning: Explicit integral could not be found.
> In sym.int at 58
In char.int at 9
>> interfcy

interfcy =
int(erfc(x)^2,x = 0 .. 1/10)

这个应该怎么改呢?哪位同学帮忙看看O(∩_∩)O~

[ 本帖最后由 ChaChing 于 2009-6-21 15:40 编辑 ]

rocwoods 发表于 2009-6-11 22:36

给你两句代码,用匿名函数表示被积函数
format long
f1= quadl(@(x) erfc(x),0,0.1)
f2= quadl(@(x) erfc(x).^2,0,0.1)
页: [1]
查看完整版本: 对erfc函数的平方怎么求定积分?