Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#4 |
普通会员
注册日期: 2009-02-22
帖子: 91
声望力: 18 ![]() |
![]()
不太理解你的意思,看下程序可否:
syms x format long e ss=zeros(10,10); for k=1:3 for l=1:3 ak=(k-1/2)*pi; al=(l-1/2)*pi; aa=sinh(ak*x)*sinh(al*x); bb=sin(al*x)*sinh(ak*x); pp=ak*al*(aa+bb); ss(k,l)=quadl(inline(pp),0,1); end end 不过,上界1应小点,不然会出现warning! |
![]() |
![]() |