c4_2.m

来自「这是我从一本基于MATLAB实现的数值分析算法的书中自带的光盘程序」· M 代码 · 共 20 行

%C4_2.m
%根据许多节点函数值的符号的变号来确定f(x)=0的含根区间
function C4_2
a=input('First x?');
b=input('Last x?');
x_intV=input('Increment of x?');
N=(b-a)/x_intV;
x=a+(0:N)*x_intV;
for i=1:N+1
    y(i)=func(x(i));
end
%fprintf('x ,y \n');
%fprintf('%2.0f,%5.3e \n',x,y);
for i=1:N
    if y(i)*y(i+1)<0
       fprintf('Interval that may contain a root:');
       fprintf('[ %2.0f,%2.0f ] \n',x(i),x(i+1));
    end
end