您所在位置：主题：[讨论]请教一个数值积分
ueryka
[专家分：10] 发布于 2008-04-09 01:06:00
dblquad(@NFP2,1/2,sqrt(2)/2,-pi/2,pi/2)
function out=NFP2(r,sita)
R=1;
eta=acos(1-(R^2./(2*(r.^2))));
xi=acos(r/R.*((cos(sita)).^2)-sin(abs(sita))/R*sqrt(R^2-r.^2.*(cos(sita)).^2));
delta=acos(R./(2*r));
B2=1/2*R^2.*xi-1/2.*r*R.*sin(xi);
B1=1/2*R^2.*delta-1/2.*r*R.*sin(delta)-B2;
out=1/2*r.^2.*(pi/2-abs(sita))+1/2*r.^2.*eta+B1;
这个是个面积分，但是为什么我算出来的答案是个负数呢？实在不明白，也没检查出明显错误，请教下大家
function out=NFP2(r,sita)
R=1;
eta=acos(1-(R^2./(2*(r.^2))));
xi=acos(r/R.*((cos(sita)).^2)-sin(abs(sita))/R*sqrt(R^2-r.^2.*(cos(sita)).^2));
delta=acos(R./(2*r));
B2=1/2*R^2.*xi-1/2.*r*R.*sin(xi);
B1=1/2*R^2.*delta-1/2.*r*R.*sin(delta)-B2;
out=1/2*r.^2.*(pi/2-abs(sita))+1/2*r.^2.*eta+B1;
这个是个面积分，但是为什么我算出来的答案是个负数呢？实在不明白，也没检查出明显错误，请教下大家