#include <math.h>
double free_freq=1.;// frequenza naturale dell'oscillatore 1 Hz
double f=0.5;       // frequenza forzante = 0.5 Hz;


void set_forzante(double freq) {
  f=freq;
}

double* osc_forzata(double *x, double t) { // oscillatore armonico
  double *derivata=new double[2];
  double omega =2.*M_PI*free_freq; 
  double omega2=2.*M_PI*f; ;
  derivata[0]=x[1];                            // dx/dt = v;
  derivata[1]=-omega*omega*x[0]+sin(omega2*t); // dv/dt = -omega^2 x+forzante
  return derivata;
}