一样的被积函数，怎么就得到不一样的结果

sigma665

%%% example03 %%%%
z=-0.1:0.01:0.1;
syms b;
a=0.1;
kesi=(z+sqrt(z.^2-a^2))./a;
f=log(1-(cos(b)+i*sin(b))./kesi); % 被积函数
I=int(f,'b',0,pi);
II=roundn(double(I),-4);


%%% example04 %%%%
z=-0.1:0.01:0.1;
a=0.1;
kesi=(z+sqrt(z.^2-a^2))./a;
for i=1:21
    f=@(b)(log(1-(cos(b)+i*sin(b))/kesi(i))); % 被积函数
    I(i)=quad(f,0,pi);
    III(i)=roundn(double(I(i)),-4);
end

程序如上，这是为什么啊

[ 本帖最后由 sigma665 于 2007-11-13 22:00 编辑 ]

donkeyxu

第二段程序，太乱

sigma665

改过来了

donkeyxu

f=@(b)(log(1-(cos(b)+i*sin(b))/kesi(i)));
中两个i,意义不一样，不能重复用i

sigma665

谢谢，楼上的
我也刚学

再次拜谢楼上，竟然知道前一个i是虚数单位

[ 本帖最后由 sigma665 于 2007-11-14 09:45 编辑 ]