#include <stdio.h> 
#include <math.h>

#define mu 0.00
#define a  1

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

  dv[0]= y;
  dv[1]= -x -mu*x*x*y+a+4*mu*y;
  
  return;
}


//Programa principal
int main(){

  int j;
  
  FILE *ptr;
  double i, v[2],t,dt,t_pre,t_max;




  //Archivo de salida
  ptr=fopen("hopf.dat","w");

  dt=0.01;
  t_max=600;


//for(j=1;j<20;j+=2){

//mu=j;
  ////Condiciones Iniciales
  
for(i=a-4;i<a+2;i+=1.9){

  v[0]=i;
  v[1]=0;
				
  t=0.;
					
  while(t<t_max){
    //Integra las ecuaciones utilizando el metodo de Runge Kutta
    rk4(ecuaciones,v,2,t,dt,mu);
    
    //Imprime la integracion
    fprintf(ptr,"%lg\t%lg\t%lg\n",t,v[0],v[1]);

    t+=dt;
  }
        fprintf(ptr,"\n");
  }
   
  fclose(ptr);
  return(0);
}