% (z,信息)= BulirschStoer (dynFun t z0, tol)
%
%使用Bulirsch-Stoer方法解决了一个初值问题。这
%的方法是理想的高精度解光滑的初始值万博 尤文图斯
%的问题。
%
%计算z (t), dz / dt = dynFun (t, z),从最初的开始
% z0状态。网格点的解决方案在托尔将是准确的。
%
%如果提供的网格是不够的,这个函数会自动
%介绍中间网格点来实现所需的精度。
%
%的输入:
% dynFun =函数处理系统动力学
% dz = dynFun (t, z)
% t =标量时间
% z =[新西兰1]=国家是列向量
% dz =[新西兰1]=状态的导数是列向量
% t =(1元)=时间网格点的向量
% z0 =[新西兰1]=初始状态向量
% tol =[新西兰1]=错误宽容以及每个维度。如果托尔是一个
%标量,那么所有维度将满足错误宽容。
%
%输出:(nt = n + 1)
% z =(新西兰元)=解决初值问题
%的信息
% . error = (nx、nt) =错误估计在每个网格点
% .nFunEval =(1元)=函数评估为每个点的数量
%
%指出:
%的实现细节:
%http://web.mit.edu/ehliu/Public/Spring2006/18.304/implementation_bulirsch_stoer.pdf
%