这个表达式：
wd=ld^(-rr);
其中的rr是什么？
另外，用int求积分，很多情况下是得不到解析结果的。你说其他积分方法精度低，请问你需要什么精度？
出错在E1赋值时用了double转换，但你定义的a1就一符号，你让matlab怎么帮你把它转成数值?
