anchor.diff:      call get_parameters(param)
      
       coeff_a              = (Omega**2 * Omega_0**2 * Sigma * Sigma_0)/2._x_precision
  -    coeff_b              = (-1._x_precision/3._x_precision) * cst_rad*T**4 * T_0**4 /rs
  +    coeff_b              = (-1._x_precision/3._x_precision) * cst_rad*T**4 * T_0**4 / rs
       coeff_c              = (-1._x_precision * param%RTM * T  *  Sigma * Sigma_0)/(2._x_precision * rs**2)
   
       call quadratic(coeff_a , coeff_b , coeff_c , H)
   
  -  !  H                    = H / rs
       rho                  = Sigma / H
       P_rad                = T**4
       P_gaz                = rho * T
  
positive.code:     integer                                                :: i     = 0
      integer                                                :: j     = 0
      real(kind = x_precision)                               :: temp  = 0.0d0
      real(kind = x_precision), parameter                    :: t_min = 9.715d-6
      real(kind = x_precision), parameter                    :: t_max = 9.715d-4
      integer,                  parameter                    :: nb_it = 50
      real(kind = x_precision)                               :: sigma = 0.0d0
      real(kind = x_precision)                               :: Smin  = 0.0d0
      real(kind = x_precision)                               :: Smax  = 0.0d0
      real(kind = x_precision)                               :: eps   = 1.0d-15
      real(kind = x_precision)                               :: omega = 0.0d0
      real(kind = x_precision)                               :: r     = 0.0d0
      real(kind = x_precision)                               :: H     = 0.0d0    
  
negative.code:                   tau_eff, P_rad, P_gaz,E_ff,Fz,f_min,Sigma_0, Omega_0,rs,T_0, rho_0)
            endif
            
            dichotomy       = (Smin + Smax) * 1.0 / 2._x_precision
  
            call variables(T, dichotomy, Omega, H, rho, cs, nu, Q_plus, Q_minus, K_ff,&
                 K_e, tau_eff, P_rad, P_gaz,E_ff,Fz,f_dichotomy,Sigma_0, Omega_0,rs,T_0, rho_0)
  