log_sin.m

来自「MATLAB的基本操作与使用」· M 代码 · 共 11 行

x1=1; x2=pi; y=1; I=0;
while abs(y)>0.00000000000001
    y1=log(x1)-sin(x1); I=I+1;
    y2=log(x2)-sin(x2);
    x=0.5*(x1+x2);  format long, disp(['x=',num2str(x)])
    y=log(x)-sin(x);
    if y*y1>0, x1=x; end
    if y*y2>0, x2=x; end
end, format long;  x, I