/* 25/04/95 */	
/* 28/04/95  avec calloc*/
/* ====================================================================== */	
/* resolution de l'attarcteur de Lorentz                                  */
/*                                                                        */
/*     x' = - s * x + s * y                                               */
/*     y' = - x * z + r * x - y                                           */
/*     z' =   x * y - b * z                                               */
/*                                                                        */ 
/*                                                                        */ 
/* Biblio                                                                 */ 
/*    E.A. Jackson 'Perspectives of non linear dynamics'  Cambridge 1991  */
#include <stdio.h>
#include <stdlib.h>
#include "math.h"
#include "rgraph.h"
#include "param.c"
/*#include "unix.h"*/
 void main(void)
{
int i,nbp,icoul,is,nf;
float x,y,z;
FILE *f,*ff;
ff= fopen("FF.OUT", "w");
nf=1;
nbp=1;
while(nbp>0) {
printf("nombre[%d] de points =? \n",nf);
scanf("%d",&nbp);
printf("nombre[%d] de points= %d \n",nf,nbp);
//printf("couleur \n");
//scanf("%d",&icoul);

f = fopen("F.OUT", "r");
for(i=1;i<=nbp;i++)
{
fscanf(f," %f %f %f %d \n",&x,&y,&z,&is);
if(is>0){is=nf;}
fprintf(ff,"%f %f %f %d\n",x,y,z,is);
}
nf=nf+1;
fclose(f);
}
fclose(ff);
printf("fin...\n");
ExitToShell;}