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

main()
{
  float L=4.7;//longueur
  float H=4.1;//largeur
  float W=3.1;//hauteur
  int p,q,r,i;//liste des indices
  float f;
  i=0;
  FILE *dat=fopen("data.dat","w");
  for (p=0;p<=100;p++)
  {
    for (q=0;q<=100;q++)
    {
      for (r=0;r<=100;r++)
      {
	if (p!=0 || q!=0 || r!=0)
	{
	  f=(340/2)*sqrt((p*p)/(L*L)+(q*q)/(W*W)+(r*r)/(H*H));
	  fprintf (dat,"%f\n",f);
	  i++;
	}
      }
    }
  }

  printf ("imax=%d",i);
}//fin du main