%
%   time point
%
t0=0
%
%   x0 - initial guess for the Newton method and for the initial condition
%
%   DP1(x(t0)-x0) = 0   index 2 case
%    D(x(t0)-x0)  = 0   index 1 case
%      x(t0)-x0   = 0   index 0 case
%
x0=[100000;0;0;12000;-2;50;2.6;-0.05]
%
%   y0 - initial guess of y := R(Dx)' for the Newton method
%
y0=[0;0;0;0;0;0]
%
%   use the functions setuser_para_... to change or reset parameters of the
%   algorithm. Type   help setuser_para_general
%                     setuser_para_mainconsi
%                     setuser_para_nlsolve
%
%setuser_para_general%('s_A',1,'s_D',1,'s_B',0)
%setuser_para_mainconsi%('given_index',2)
%setuser_para_nlsolve
%
%   with additional user requirements (in 'ztraj')
%
main_consi('ftraj',y0,x0,t0,'ztraj');
%
%   standard call of main_consi (with respect of 'z')
%
%main_consi('ftraj',y0,x0,t0,'z');