#ifndef harmonic_h
#define harmonic_h
#include "math.h"
// calculate harmonic oscillator base functions, for 1d and 3d
// calculate also the kinetic energy part of the wave function,
// using H=T+V hence T=E-V

class harmonic {
	public:
		harmonic(double m=1, double c=1) {mass=m; string=c;};
		~harmonic() {};
		double hermite(int k, double x);
		double wav1d(int k, double x);
		double wav3d(int k, int l, double x);
		void setstring(double par) {string=par;};
		void setmass(double par) {mass=par;};
		double secdiv(int k, double x);
		double E(int k) { return (k+0.5)*sqrt(string/mass);};
	private:
		double mass;
		double string;
};
#endif