/tmp/ccCg8Hvm.o: In function `planck': plantes.c:(.text+0x59): u

#include <stdio.h>
 #include <stdlib.h>
 #include <tgmath.h> // fonctions generiques
 #include <float.h>
 /* standard c99 pour sqrtf expf et fabsf */

 float planck(float x); /* déclaration */

 float planck(float x){
 float y;
 if(fabs(x) > sqrt(FLT_EPSILON)) {
 // appelle fabsf et sqrtf grace a tgmath.h
 y = x*x*x / (exp(x) - 1.f);
 }
 else{
 y = x * x ;
 }
 return y;
 }

 int main(void){
 int n, i;
 float x, xmin, xmax, dx ;
 float y ;
 xmin = 0.f;
 xmax = 10.f;
 dx = .5f;
 n = 1 + (int) ((xmax - xmin) / dx) ;
 for (i =1; i<= n; i++){
 x = xmin + (i -1) * dx ;
 y = planck(x) ;
 printf("%g %g\n", x, y);
 }
 exit(EXIT_SUCCESS) ;
}