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

// #define mu -1


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

  dv[0]= mu*x-4*x*x*x;
/*  dv[1]= x + a*y*(x*x + y*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("ej2.dat","w");

  dt=0.01;
  t_max=2;

double mu;

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

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

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

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