See

• The implementation ode.cpp
• The description ODE

<project-file type=“source”/>

<content> #ifndef ODE_H #define ODE_H

#include <vector> #include <networkstructure.h> #include <equations.h>

class ODEResults { public:

  ODEResults(int max_household_size, int max_external_degree);
  void addPoint(double t, double y[]);
  float t(int i) const { return T[i]; }
  float It(int i) const;
  float IH(int i) const;
  int points() const { return T.size(); }

private:

  int MaxHouseholdSize;
  int MaxExternalDegree;
  std::vector<float> T;
  std::vector<std::vector<float> > H;
  std::vector<std::vector<float> > S;
  std::vector<std::vector<float> > I;

};

ODEResults *ode_simulation(ODEParameters &p, double t0, const NetworkStructure &net);

#endif ODE_H </content> <use name=“equations.h”/> <use name=“networkstructure.h”/>