为什么有时能运行有时不能

zqs310

n=1;S=[0 0 0 0 0 0];
global axb ac2 ac1 arr
for xb=2.0:0.2:3.0
    c2=xb/2;
   axb=xb;ac2=c2;
    rr=fzero('fr',[c2 3.1]);
    c1=-sqrt(rr^2-c2^2);
    ac1=c1;
    arr=rr;      
    S(n)=quadl('ff',0,xb);
    n=n+1;
   end
   
function y=fr(x)
global axb ac2
y=x.*(asin((axb-ac2)./x)+asin(ac2./x))-pi;

function y=ff(x)
global arr ac2 ac1
y=sqrt(arr^2-(x-ac2).^2)+ac1;
有时运行成功，有时提示如下："
??? Error using ==> fzero
Function values at interval endpoints must be finite and real.
Error in ==> dzz at 6
    rr=fzero('fr',[c2 3.1]);".
多运行几次就行。
并且都是在第一次循环时出错的。
我是新手，请大虾们帮忙！！！

ChaChing

无法复制楼主的情状
但猜测是否LZ函数fr中用的asin内容超过1, asin(x)已是复数

zqs310

现在一次都没有出现了那种问题了。可能第一次的时候改过了，现在好了。