% function [xout, tout]=bs_lorenz(x0)
%   Solves the lorenz system for initial conditions given in x0
%   dX/dt = -aX   +aY
%   dY/dt = +rX   -Y       -XZ
%   dZ/dt =           -bZ  +XY
%   
%   a=10.;  r=28.;    b=8./3.;
%
function [xout, tout]=bs_lorenz(x0, t0, tf)

%  solve the system using the built-in ode solver ode45
[tout, xout] = ode45(@bs_lorenz_rhs, [t0, tf], x0);