|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
%%% 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 编辑 ] |
|