ncc198611
2009-05-17, 17:41
syms x;
for k=1:10
for l=1:10
ak=(k-1/2)*pi;
al=(l-1/2)*pi;
num2str(ak);
aa=strcat('(sinh((',num2str(ak),')*x).*sinh((',num2str(al),')*x))');
bb=strcat('(sin((',num2str(al),')*x).*sinh((',num2str(ak),')*x))');
pp=strcat('(',num2str(ak),').*(',num2str(al),').*((inline(aa,x))+(inline(bb,x)))')
ss(k,l)=quadl(inline(pp),0,1)
end
end
请高手帮忙看一下有什么错误,帮改一下,如果把aa,bb这项直接带到pp中能出结果,因为像aa,bb这样的项有很多,而且在strcat('( )')又不能分行,函数太长不便于打印。谢谢了。
for k=1:10
for l=1:10
ak=(k-1/2)*pi;
al=(l-1/2)*pi;
num2str(ak);
aa=strcat('(sinh((',num2str(ak),')*x).*sinh((',num2str(al),')*x))');
bb=strcat('(sin((',num2str(al),')*x).*sinh((',num2str(ak),')*x))');
pp=strcat('(',num2str(ak),').*(',num2str(al),').*((inline(aa,x))+(inline(bb,x)))')
ss(k,l)=quadl(inline(pp),0,1)
end
end
请高手帮忙看一下有什么错误,帮改一下,如果把aa,bb这项直接带到pp中能出结果,因为像aa,bb这样的项有很多,而且在strcat('( )')又不能分行,函数太长不便于打印。谢谢了。