[TTT_IOU]

不太理解你的意思,看下程序可否：
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！