// In header: <boost/numeric/odeint/stepper/rosenbrock4.hpp>
template<typename Value,
         typename Coefficients = default_rosenbrock_coefficients< Value >,
         typename Resizer = initially_resizer>
class rosenbrock4 {
public:
  // types
  typedef Value                                               value_type;
  typedef boost::numeric::ublas::vector< value_type >         state_type;
  typedef state_type                                          deriv_type;
  typedef value_type                                          time_type;
  typedef boost::numeric::ublas::matrix< value_type >         matrix_type;
  typedef boost::numeric::ublas::permutation_matrix< size_t > pmatrix_type;
  typedef Resizer                                             resizer_type;
  typedef Coefficients                                        rosenbrock_coefficients;
  typedef stepper_tag                                         stepper_category;
  typedef unsigned short                                      order_type;
  typedef state_wrapper< state_type >                         wrapped_state_type;
  typedef state_wrapper< deriv_type >                         wrapped_deriv_type;
  typedef state_wrapper< matrix_type >                        wrapped_matrix_type;
  typedef state_wrapper< pmatrix_type >                       wrapped_pmatrix_type;
  typedef rosenbrock4< Value, Coefficients, Resizer >         stepper_type;
  // construct/copy/destruct
  rosenbrock4(void);
  // public member functions
  order_type order() const;
  template<typename System>
    void do_step(System, const state_type &, time_type, state_type &,
                 time_type, state_type &);
  template<typename System>
    void do_step(System, state_type &, time_type, time_type, state_type &);
  template<typename System>
    void do_step(System, const state_type &, time_type, state_type &,
                 time_type);
  template<typename System>
    void do_step(System, state_type &, time_type, time_type);
  void prepare_dense_output();
  void calc_state(time_type, state_type &, const state_type &, time_type,
                  const state_type &, time_type);
  template<typename StateType> void adjust_size(const StateType &);
  // protected member functions
  template<typename StateIn> bool resize_impl(const StateIn &);
  template<typename StateIn> bool resize_x_err(const StateIn &);
  // public data members
  static const order_type stepper_order;
  static const order_type error_order;
};