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



void ecuaciones(int n, double v[], double dv[], double t,double A, double w )
{
  double x,y,z;
  x=v[0];y=v[1];z=v[2];

  dv[0]= 10*(y-x);
  dv[1]= 28*x-y-x*z;
  dv[2]=x*y-8*z/3.; 

  return;
}


int main(){

	FILE *ptr;
	double v[3],t,dt,t_pre,t_max;

	ptr=fopen("salida.dat","w");


	v[0]=10.;
	v[1]=10.;
	v[2]=40.;
        
	
	dt=0.001;
		
	t=0.;
	t_pre=0;
	t_max=50.;
			
			
	while(t<t_max){
		rk4(ecuaciones,v,3,t,dt);
		if(t>t_pre) fprintf(ptr,"%lg\t%lg\t%lg\n",v[0],v[1],v[2]);
		t+=dt;
	}	

	fclose(ptr);
	return(0);

}