Calculo numérico - erro metodo da bissecao projeto em c++

#include <iostream>
#include <stdlib.h>
#include <stdio.h>
#include <math.h> using namespace std;

void Inicio (){ cout << "\n \t\tJULIANA, LEONARDO, MARIA CARLA, ROMULO, EVERTON \n\n"; cout << "\t\t *********************************************\n\t\t| PROJETO COMPUTACIONAL - BISSECAO |\n\t\t *********************************************\n\n";
}
int main()
{
double l, c, X[1000]; double Ri[1000]; double Si[1000]; Ri[0]=0; double CritPara1, CritPara2, Volume; int i=0; Inicio(); cout << "\nPara o projeto da caixa, siga as instrucoes abaixo:\n\n(Valores dados pelo exercicio: largura 40 e comprimento 52)"; cout <<"\n\n\nInsira a largura: "; cin>>l;//LARGURA cout << "\nInsira o comprimento: " ; cin>>c; //COMPRIMENTO Si[0]=l;

for(i=0; i<13; i++){ X[i] = (Ri[i] + Si[i]) /2; //DERIVADA DA FUNÇÃO VOLUME CritPara1 = (((12 * Ri[i]*Ri[i]) - (4*c*Ri[i]) - (4*l*Ri[i]) + l*c) * ((12 * X[i]*X[i]) - (4*c*X[i]) - (4*l*X[i]) + l*c)); if(CritPara1 < 0){ Ri[i+1] = Ri[i]; Si[i+1] = X[i]; if(i != 0){ CritPara2 = (X[i] - X[i-1]) / (X[i]); if (CritPara2 < 0){ CritPara2 = CritPara2*(-1); } } }else{ Ri[i+1] = X[i]; Si[i+1] = Si[i]; if(i != 0){ CritPara2 = (X[i] - X[i-1]) / (X[i]); if (CritPara2 < 0){ CritPara2 = CritPara2*(-1); } } } }

system ("cls"); cout << "*******************************************************************************\n\t\t\t | RESULTADOS |\n*******************************************************************************"; cout<<