关于积分quad2d函数的用法问题

江山美人

编了个程序如下：积分之后求和的问题，会弹出错误使用 quad2d (line 114)C 必须为有限标量浮点常量或函数句柄。请各位大神指点一下，怎么修改，谢谢！

1. syms x y k
   
2. a=1;
   
3. b=2;
   
4. c=3;
   
5. d=2;
   
6. br=0.1;
   
7. hr=0.2;
   
8. h=0.01;
   
9. Nr=4;
   
10. t1=-k./(Nr+1)-br./2;
    
11. t2=-k./(Nr+1)+br./2;
    
12. fr=quad2d(@(x,y) (1+y).^2.*x.^(a+c-1).*y.^(b+d),1+h/2,1+h/2+hr,t1,t2);
    
13. z=symsum(fr,k,1,Nr-1)
    

复制代码

chybeyond

把t1和t2改为函数句柄形式t1=@(k) -k./(Nr+1)-br./2;
t2=@(k) -k./(Nr+1)+br./2;

江山美人

chybeyond 发表于 2014-6-19 10:26
把t1和t2改为函数句柄形式t1=@(k) -k./(Nr+1)-br./2;
t2=@(k) -k./(Nr+1)+br./2;

OK了，谢谢啊