为什么有时能运行有时不能
n=1;S=;global axb ac2 ac1 arr
for xb=2.0:0.2:3.0
c2=xb/2;
axb=xb;ac2=c2;
rr=fzero('fr',);
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',);".
多运行几次就行。
并且都是在第一次循环时出错的。
我是新手，请大虾们帮忙！！！ 无法复制楼主的情状
但猜测是否LZ函数fr中用的asin内容超过1, asin(x)已是复数
感谢
现在一次都没有出现了那种问题了。可能第一次的时候改过了，现在好了。
页:
[1]