qujian.m

来自「应用模拟退火方法求解非线形方程」· M 代码 · 共 20 行

function [X]=qujian(L,R,h)
fprintf('根的搜索区间为：[%g,%g]\n',L,R)
p=L;q=R;k=1;
X=zeros(1,2);
while L+h<=R
    if f(L)*f(L+h)<=0
        fprintf('有根区间%d：[%g,%g]\n',k,L,L+h);
        X(k,:)=[L,L+h];
        k=k+1;
        if f(L)==0
            L=L+h/2;
        end
    end
    L=L+h;
end
t=p:0.1:q;
for i=1:length(t)
    y(i)=f(t(i));y1(i)=0;
end
plot(t,y,t,y1)