Funcion de tercer grado

Este programa nos permite encontrar los valor de f(x) de una función de tercer grado o cúbica para los valores de x que el usuario acota y su da el valor de incremento.

1. //funcion de grado 3
2. #include<stdio.h>
3. #include<math.h>
4. int main(){
5.    printf("Funcion de grado 3\n"); 
6.    double a[4],x,xmin,xmax,xincr,y,n;
7.    int i,condicion1;   
8.    do{
9.       for(i=0;i<4;i++){
10.       printf("ingrese el coeficiente %d\n",i);   
11.       scanf("%lf",&a[i]); 
12.       }  
13.      do{
14.         printf("ingrese el x minimo\n");
15.         scanf("%lf",&xmin);
16.         printf("ingrese el x maximo\n");
17.         scanf("%lf",&xmax);
18.         printf("ingrese el incremento de x\n");
19.         scanf("%lf",&xincr);     
20.         for(x=xmin;x<=xmax;x=x+xincr){    
21.           n=3;
22.           y=0;                         
23.           for(i=0;i<4;i++){
24.              y=a[i]*pow(x,n)+y;
25.              n--;          
26.           } 
27.           printf("f(%lf.8)=%lf\n",x,y);         
28.         }        
29.         printf("hallar para otros valores de x\n presione\n");
30.         printf(" 1 continuar\n 2 finalizar\n");
31.         scanf("%d",&condicion1);
32.      }while(condicion1==1);
33.      printf("hallar para otros valores de los coeficientes\n presione\n");
34.      printf(" 1 continuar\n 2 finalizar\n");
35.      scanf("%d",&condicion1);             
36.      }while(condicion1==1);
37.      return 0;
38. }