如何对代码中大量内容进行修改

burkbai

前几天发帖问过一次，现在遇到新的问题，下面是一个积分运算
clear
syms E X Nrd r t Rd Rr P0 r1 r2 N1 N2 N3 Sec G d;
d=0.7;
p=8/pi;
Rr=1;
%Rd=10^(1/4)*Rr;
Rd=2*Rr;
X=(Rd^2+r^2-Rr^2)/(2*r);
Sec=(Rd^2*acos(X/Rd)+Rr^2*acos((r-X)/Rr)-r*(Rd^2-X^2)^0.5);
P1=Sec/pi/Rr^2;
r1=d*cos(G)+((d*cos(G))^2-d^2+Rd^2)^0.5;
r2=Rd+Rr;
N1=int(p*r*P1,'r',r1,r2) %N1是G的一个大函数
在这里matlab算出来 N1 是一个G为参数的大函数，大约有几万个字符，是有无数个包含cos(G)的复杂表达式，
现在需要对G积分区间是【0 ，2*pi】，
我想用quadl积分，下面这么写，错误提示很多地方需要改为点乘，可是里面的内容是在太多，如何修改呢？谢谢！
myfun = @(G) (N1大函数的字符串);
Q = quadl(myfun,0,2*pi )

[ 本帖最后由 eight 于 2008-1-22 20:02 编辑 ]

sogooda

查找，替换？

sigma665

复杂表达式,建议还是一开始就数值积分

burkbai

多谢指点，在这里找到了方法
http://biosanger.spaces.live.com ... 9e6e2616f!276.entry