HarmonicOscillator.h

Go to the documentation of this file.

#ifndef _HarmonicOscillator_h_
#define _HarmonicOscillator_h_


class HarmonicOscillator {
public:
  HarmonicOscillator(const double m, const double k)
    : fM(m), fK(k) { }

  template<typename Vector>
  void
  operator()(const double /*x*/, const Vector& y, Vector& dYdX)
    const
  {
    dYdX[0] = y[1] / fM;
    dYdX[1] = -fK * y[0];
  }

  operator unsigned int() const
  { return 2; }

private:
  const double fM;
  const double fK;
};


#endif