Решение системы уравнений методом прогонки

int main(int argc, char* argv[]){
    
 ifstream inp("input.txt"); // открыли файл для чтения
 ofstream out("output.txt");
 if(!inp){
    cout<<"Open file error\n";
    system("pause");
   return 1;
};
 
 int arraySize = 0;
 inp>>arraySize;
 double mass[105][105];
 //Эти два вложенных цикла используются для ввода элементов массива
for (int i=0; i<arraySize; i++)
{
 
for (int j=0; j<arraySize; j++)
{ inp>>mass[i][j];}
 
}
double R[105];
for (int k=0; k<arraySize;k++){
    inp>>R[k];}
 
for (int i=0; i<arraySize; i++)
{
 
for (int j=0; j<arraySize; j++)
{ out<<mass[i][j]<<" ";}
out<<endl;
}
out<<endl;
out<<"vector R   ";double delta[105],lambda[105];
for (int k=0; k<arraySize;k++){
    out<<R[k]<<" ";}
out<<endl;
//нахождение коэффициентов лямбда и дельта в разложении х(i-1)=delta(i)*x(i)+lambda(i)
for (int i=1; i<arraySize; i++) 
{ delta[0]=-mass[0][1]/mass[0][0];lambda[0]=R[0]/mass[0][0];
delta[arraySize]=0;mass[i][arraySize]=0;
    delta[i]=-(mass[i][i+1])/(mass[i][i]+(delta[i-1]*mass[i][i-1]));
   lambda[i]=(R[i]-mass[i][i-1]*lambda[i-1])/(mass[i][i]-mass[i][i-1]*delta[i-1]);
   out<<"delta"<<i<<"="<<delta[i]<<endl;
   out<<"lambda"<<i<<"="<<lambda[i]<<endl;
out<<"delta0"<<"="<<delta[0]<<endl;
   out<<"lambda0"<<"="<<lambda[0]<<endl;
 
 
   for (int i=arraySize-1;i>=0;i--){
   double x[105];
   x[arraySize-1]=lambda[arraySize-1];
   x[i]=(delta[i]*x[i+1])+lambda[i];
   out<<"x"<<i<<"="<<x[i]<<endl;}}