对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

给你两句代码，用匿名函数表示被积函数

1. 
   
2. format long
   
3. f1= quadl(@(x) erfc(x),0,0.1)
   
4. f2= quadl(@(x) erfc(x).^2,0,0.1)
   

复制代码