#include "help.hpp"
  inline void help::fun(const std::vector<double>& y, 
                  const std::vector<double>& theta_dbl_, 
                  double t, std::vector<double>& ydot) const {
    using std::pow;
    #include "ode.txt"
  }
  inline void help::jacobianTheta(const std::vector<double>& y, 
                            const std::vector<double>& theta_dbl_, 
                            double t, std::vector<double>& Jtheta) const {
    // Jtheta - df1/dtheta1 ... dfn/dtheta1 ...
    // total length: N_ * M_
    // dfn/dthetam=Jtheta[m * M_ + n]
    using std::pow;
    #include "params.txt"
  }
  inline void help::jacobianY(const std::vector<double>& y, const std::vector<double>& theta_dbl_, 
                        double t, std::vector<double>& Jy) const {
    // dyn/dt = fn(y, theta, t)
    // Jy - df1/dy1 ... dfn/dy1 ...
    // total length: N_ * N_
    // dfn/dyj=Jy[j * N_ + n]
    using std::pow;
    #include "states.txt"
  }