#include <stdio.h> 
#include <math.h>
#define N_NAME1_AUX 50
#define a -1

// define parametros para usar en el sistem en todo el codigo
struct Par{
  double  gamma, omega;
} aa;


// Ecuaciones del sistema
void ecuaciones(int n, double v[], double dv[], double t){
  double x,y;
  x=v[0];/*y=v[1];*/
  y=v[1];

  dv[0]= y;
  dv[1]= - 2*aa.gamma*y - aa.omega*aa.omega*x;
/*  dv[1]= x + a*y*(x*x + y*y);*/
  
  return;
}

//----------------funcion que genera muchos nombres segun parametros--//
void nombrar1_par1(char *name_space1,int par1){
  sprintf(name_space1,"oscilador.%d.dat",par1);
return;
}

//Programa principal
int main(){

  int i,j;
  FILE *ptr, *ptr1;
  double v[2],t,dt,t_pre,t_max;
  char nombre_archivo1[N_NAME1_AUX]; // define una variable del tipo string
  

//Archivo de salida del comando para graficar todo junto en gnuplot
  ptr=fopen("valeria_dibujando.plt","w");

  dt=0.01;
  t_max=4;

  //Condiciones Iniciales
  
  v[0]=1;
  v[1]=0;
  aa.gamma=0;
  aa.omega=1;

//----para hacer muchos archivos--//

// loop sobre el valor de omega--//
for(j=1;j<15;j++){aa.omega=j;
 // genera muchos archivos segun el parametro de iteracion j//				
  nombrar1_par1(nombre_archivo1,j);
  ptr1=fopen(nombre_archivo1,"w");

  t=0;				
  while(t<t_max){
    //Integra las ecuaciones utilizando el metodo de Runge Kutta
    rk4(ecuaciones,v,2,t,dt);
    
    //Imprime la integracion
   
    fprintf(ptr1,"%lg\t%lg\n",t,v[0]);
    
  
    t+=dt;
   }
//-- Crea un archivo donde tiene una instruccion para graficar de a uno y vas pasando con enter
 fprintf(ptr, "p [][]'%s' u 1:2 w l\n pause -1\n", nombre_archivo1);
  }    
  fclose(ptr);
  return(0);
}