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

double sigma=10, rho=28, beta=8.0/3.0;
double T=10;
int n=1000;

void f(double Y[3],double X[3]){
	Y[0]=sigma*(X[1]-X[0]);
	Y[1]=X[0]*(rho-X[2])-X[1];
	Y[2]=X[0]*X[1]-beta*X[2];
}

int main(){
	double Y[3]= { 1, 2, 3 };
	double F[3];
	double h=T/n;

	for(int i=0;i<=n;i++) {
 		printf("%g %g %g %g\n",
			i*h,Y[0],Y[1],Y[2]);
		f(F,Y);
		for(int j=0;j<3;j++) Y[j]+=h*F[j];
	}
	return 0;
}