方程求根

function xy=findrootinterval ( fun,xmin,xmax,n )

%findrootinterval 在区间[xmin,xmax]上求函数 f(x) 的有限根区间

%输入参数：fun： 字符串，可以是表达函数 f(x) 的.m文件名，字符表达式或符号变量表达式

%                  max：初始搜索区间的左、右端点

%                      n：（可选）等分区间的份数 缺省：n=20

%output：xy：两列元素的矩阵，保存有根区间的左、右端点。若没找到有根区间，则xy=[ ].

%这个算法的主要目的是确定有根区间，用它来求根当然速度不快

==========================================================================

===========================================================================

二分法求根：

function x = bisection(fun,a,b,xtol,ftol,n,action)

%bisection 用二分法求f (x) = 0 的根

%输入参数： fun ：字符串，可以使表示函数 f(x)  的.m文件名，字符表达式或符号变量表达式

%                     a，b：初始有根区间的左，右端点

%                      xtol：（可选）根的容许误差。缺省：xtol=5*eps

%                   ftol：（可选）函数f(x) 的容许误差，缺省：ftol=5*eps

%                      n：（可选）最大的二分次数，缺省：n=50

%                    action：（可选)决定是否输出中间结果，缺省：action=0，不输出

%输出参数：x：用二分法求得的 f(x) = 0 的近似根

%二分法的优点是简单可靠，总能收敛。

%二分法的缺点一是首先确定有根区间，并非易事，而是收敛速度太慢

========================================================================