Trouver les extremums d'une fonctions polynome

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

int main()
{
	printf(" integration par tirage aleatoire \n");
	int D;
	float Tab[50];
	int ind;
	float P;
	float xmax;
	float xmin;
	float ymax;
	float ymin;
	float i;
	int ii;
	float yxi;
	float x;
	float ax;
	float ay;
	float ar;
	int N;
	time_t u;
	int compteur;
	float xa;
	float ya;
	float y;
	float vallapp;
	float numerateur;
	time(&u);
	srand(u);
	printf("entrer le degre du poly attention le degre est inferieur ou egal a 20 : ");
	scanf("%d", &D);
	for (ind = 0; ind <= D; ind = ind + 1)
	{
		printf(" entrer le coeff de degre %d : ", D - ind);
		scanf("%f", &Tab[ind]);

	}
	xmax = 1;
	xmin = -1;
	ymax = -10000;
	ymin = 10000;
	for (i = -1; i <= 1; i = i + 0.0000001)
	{
		x = i;
		P = 0;
		for (ind = 0; ind <= D; ind = ind + 1)
		{
			P = P*x + Tab[ind];
		}
		yxi = P;
		if (yxi <= ymin)
		{
			ymin = yxi;
		}
		if (yxi >= ymax)
		{
			ymax = yxi;
		}
	}
	printf("le ymax est %f le min %f ", ymax, ymin);
	ax = xmax - xmin;
	ay = ymax - ymin;
	ar = ax*ay;
	printf(" l'air du rectangle est %f ", ar);
	printf(" entrer le nombre de tirage souhaite");
	scanf("%d", &N);
	compteur = 0;
	for (ii = 1; ii <= N; ii = ii + 1)
	{
		xa = ((rand()*(xmax - xmin)) / RAND_MAX)  + xmin;
		ya = ((rand()*(ymax - ymin)) / RAND_MAX) + ymin;
		x = xa;
		P = 0;
		for (ind = 0; ind <= D; ind = ind + 1)
		{
			P = P*x + Tab[ind];
		}
		y = P;
		if ((ya < y) && (ya > 0))
		{
			compteur = compteur + 1;
		}
		else if ((ya > y) && (ya < 0))
		{
			compteur = compteur - 1;
		}
	}
	printf(" le compteur est %d ", compteur);
	numerateur = compteur*ar;
	float delta;
	delta = 0;
	if (ymin > 0)
		delta = ymin;
	else if (ymax < 0)
		delta = ymax;
	vallapp = numerateur / N+(xmax-xmin)*delta;
	printf(" la valeur approche est %f ", vallapp);
	scanf("%f", ymin);
}

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

typedef struct {
 float * coeff;
 int degree;
} polynome;

#define DEGREE 3

float x_pow ( float x , int degree) {
 int i;
 float pow = 1.0;
 for (i = 0 ; i < degree;i++){
  pow *= x;
 }
 return pow;
}


float polynome_x ( polynome poly , float x){
 float y = 0; 
 int i;
 
 for ( i = 0 ; i < poly.degree ; i++ ) {
  y += x_pow(x , i);
 }
 
 return y;
}

#define LIMITE 1000.0
#define PAS 0.001
int main(void)
{
 polynome my_polynome;
 int i;
 
 my_polynome.degree = DEGREE;
 my_polynome.coeff = (float * ) malloc (DEGREE * sizeof(float)) ;
 
 // On le remplis de manière simple pour l'exemple
 for (i=0;i<DEGREE;i++){
  my_polynome.coeff[i] = (float ) i;
 }
 
 
 float x_p = - LIMITE;
 float x_max;
 float max = polynome_x(my_polynome,x_p);
 float current;
 while ( x_p < LIMITE ){
  current = polynome_x(my_polynome,x_p);
  if ( current > max){
   max = current;
   x_max = x_p;
  }
  x_p += PAS;
 }
 
 printf("maximum en x = %f  avec P(x) = %f \n", x_max , max );
 
 return 0;
}

printf("entrer le degre du poly attention le degre est inferieur ou egal a 20 : ");
scanf("%d", &D);

printf(" integration par tirage aleatoire \n");
    int D;
    float Tab[500000];
    int ind;
   float P;
  int xmax;
   int xmin;
  float ymax;
   float ymin;
   float i;
  int ii;
 float yxi;
 float x;
 float ax;
 float ay;
 float ar;
 int N;
 long int u;
 int compteur;
 float xa;
 float ya;
 float y;
 float vallapp;
 float numerateur;




time(&u);
srand(u);

 printf("entrer le degre du poly attention le degre est inferieur ou egal a 20 : ");
 scanf("%d", &D);
 for (ind=0; ind<=D; ind=ind+1)
 {
 printf(" entrer le coeff de degre %d : ", D-ind);
 scanf("%f", &Tab[ind]);

 }

 P=0;



 xmax=1;
 xmin=-1;
 ymax=-10000;
 ymin=10000;
 for (i=-1; i<=1; i=i+0.0000001)
 {
x=i;

 for (ind=0; ind<=D; ind=ind+1)
 {
 P=P*x+Tab[ind];
 }


 yxi=P;
 P=0;
 if (yxi<=ymin)
 {
 ymin=yxi;
 }
 if (yxi>=ymax)
 {
 ymax=yxi;
 }
 }
 printf("le ymax est %f le min %f ", ymax, ymin);
ax=xmax-xmin;
ay=ymax-ymin;
ar=ax*ay;
printf(" l'air du rectangle est %f ", ar);
printf(" entrer le nombre de tirage souhaite");
scanf("%d", &N);

compteur=0;

for (ii=1; ii<=N; ii=ii+1)
{
    xa=((rand())/RAND_MAX)*2-1;
    ya=((rand()*(ymax-ymin))/RAND_MAX)+ymin;
    x=xa;

   for (ind=0; ind<=D; ind=ind+1)
 {
 P=P*x+Tab[ind];
 }
 y=P;
 P=0;

 if  ((ya<y) && (y>0))
 {
     compteur=compteur+1;
 }
else if ((ya>y) && (y<0))
 {
     compteur=compteur-1;
 }


}
printf(" le compteur est %d ", compteur);

numerateur=compteur*ar;



vallapp=numerateur/N;
printf(" la valeur approche est %f ", vallapp);









 scanf("%f", ymin);

if  ((ya<y) && (y>0))
 {
     compteur=compteur+1;
 }
else if ((ya>y) && (y<0))
 {
     compteur=compteur-1;
 }

if  ((ya<y) && (ya>0))
 {
     compteur=compteur+1;
 }
else if ((ya>y) && (ya<0))
 {
     compteur=compteur+1;
 }

xa=((rand())/RAND_MAX)*2-1;

xa=((rand()*(xmax-xmin))/RAND_MAX)+xmin;