Решение уравнения методом ломаных

    
        int i,L=16;
    float x[100],p[100],d[100],E[100],x1[100],x2[100],P[100],e;
    printf("\n");
    cout << "Начало отрезка а = "; cin >> a;
    cout << "Конец отрезка b = "; cin >> b;
    cout << "точность e = "; cin >> e;
 
    x[0]=(function(a)-function(b)+L*(a+b))/2/L;
    //printf("%7.2f",x[0]);
 
    p[0]=(function(a)+function(b)+L*(a-b))/2;
    //printf("%7.2f",p[0]);
    printf("|   n   |   x*   |    p*    |  delta  |   E   |   x1'  |   x2'  |    P   |\n");
    for (i=0;i<60;i++)
        {d[i]=(function(x[i])-p[i])/2/L;
    //printf ("%7.2f",d[i]);
    
    E[i]=2*L*d[i];
    //printf ("%7.2f",E[i]);
    
    x1[i]=x[i]-d[i];
    x2[i]=x[i]+d[i];
    //printf ("%7.2f %7.2f",x1[i],x2[i]);
    
    P[i]=(function(x[i])+p[i])/2;
    //printf ("%7.2f",P[i]);
    printf("\n");
    x[2*i+1]=x1[i];
    x[2*i+2]=x2[i];
    p[2*i+1]=P[i];
    p[2*i+2]=P[i];
    //printf ("%7.2f %7.2f %7.2f %7.2f",x[2*i+1],x[2*i+2],p[2*i+1],p[2*i+2]);
 
printf("|   %d   | %6.2f |  %6.2f  | %7.2f | %5.2f | %6.2f | %6.2f | %6.2f |\n",i,x[i],p[i],d[i],E[i],x1[i],x2[i],P[i]);
 
 
    }
    
 
    printf("\n");
    system("pause");
 return 0;

float x[100],p[100],d[100],E[100],x1[100],x2[100],P[100],e;

 x[2*i+1]=x1[i];
    x[2*i+2]=x2[i];
    p[2*i+1]=P[i];
    p[2*i+2]=P[i];