Ошибка в программе, скорей всего где-то скобки не поставил или пробелов не хватает, но это не точно

__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
    pi=3.141592654;
    electron=1.6e-19; /* заряд электрона */
    NAvogadro=6.023e23; /* число Авогадро, шт./моль */
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::Button1Click(TObject *Sender)
{
 
//тут идет считывание данных с тхт документа
 int j;
 FILE*fp;
   if (InputOpen->Execute()) {
      if(!(fp = fopen(InputOpen->FileName.c_str(),"rt")))
         {Application->MessageBoxA("Файл не открылся", "Ошибка", MB_OK);}
      else
         {int ii=fscanf(fp," m= %d "
           " ro= %le atom= %le "
           " l1= %le "
           " R= %le "
           " Q= %le U= %le kcur= %le kenergy= %le "
           " Ysput= %le Yself_sput= %le Mion= %le "
           " Esput= %le Eself_sput= %le "
           " S= %8s "
           ,&m,&ro,&atom,&l1,&R,&Q,&U,&kcur,&kenergy,&Ysput,&Yself_sput,&Mion,&Esput,&Eself_sput,S);
          if (ii!=15) Application->MessageBoxA("Не все данные введены",
                 "Ошибка", MB_OK);
       fclose(fp);
       }
   }
 
 
   n0=ro*NAvogadro/atom;   /* ядерная плотность */
 
   Edit6->Text=AnsiString(FormatFloat("0.000E+00",R));
   Edit7->Text=AnsiString(FormatFloat("0.000E+00",R));
   Edit8->Text=AnsiString(FormatFloat("0.000E+00",R));
 
   r.resize(m);  /* координаты ячеек вдоль поверхности мишени */
   q.resize(m); /* плотность мощности */
   vsput.resize(m); /* скорость распыления при бомбардировке ионов рабочего газа */
 
   // Ввод данных пространственного распределения плотности мощности вдоль поверхности мишени
   vec FR, FF; /* для задания распределения плотности тока по r */
   AnsiString S_help;
   S_help=AnsiString(S)+".txt";
   ifstream input(S_help.c_str());
   double value1;
   input >> value1;
   while (!input.eof())
   {
      FR.push_back(value1);
      input >> value1;
      FF.push_back(value1);
      input >> value1;
   }
   nist = FF.size();
 
 
// Расчет одной величины и вбитие этой величины в цифрку которая будет выводить это в файл 
   vec f;
   f.resize(m);
   for (j=0; j<m; j++)
      f[j]=approxd(r[j],nist,FR.begin(),FF.begin());
   double B=0;
   for (j=0; j<m; j++)
   hr=R/m;
   r[j]=hr/2+hr*j;
      B=B+r[j]*f[j]*hr;
        double A=Q/(2.0*pi*B);
        q[j]=f[j]*A;
 
//Расчет апроксимаций мощности
double P_approx=0.0;
for (j=0; j<m; j++)
  P_approx=P_approx+r[j]*q[j];
P_approx=P_approx*hr*2.0*pi;
Edit2->Text=AnsiString(FormatFloat("0.00E+00",P_approx));
 
 
//Вывод результатов в файл(ВОТ ТУТ ВЫВОДИТ НУЛИ в Файл)
FILE* file = fopen("result.txt", "wt");
if (file == NULL)
{
Application->MessageBoxA("Файл не создался", "Ошибка", MB_OK);
return;
}
 
 
for (j = 0; j < m; j++)
{
  vsput[j]=(q[j]/(U*electron))*kcur*kenergy*Ysput/n0;
  fprintf(file, "%.1f %.1f\n", vsput[j], r[j]);
}
fclose(file);
}

//Вывод результатов в файл(ВОТ ТУТ ВЫВОДИТ НУЛИ в Файл)
FILE* file = fopen("result.txt", "wt");
if (file == NULL)
{
Application->MessageBoxA("Файл не создался", "Ошибка", MB_OK);
return;
}
 
 
for (j = 0; j < m; j++)
{
  vsput[j]=(q[j]/(U*electron))*kcur*kenergy*Ysput/n0;
  fprintf(file, "%.1f %.1f\n", vsput[j], r[j]);
}
fclose(file);
}

   for (j=0; j<m; j++)
   hr=R/m;
   r[j]=hr/2+hr*j;
      B=B+r[j]*f[j]*hr;
        double A=Q/(2.0*pi*B);
        q[j]=f[j]*A;

  
for (j=0; j<m; j++)
{
   hr=R/m;
   r[j]=hr/2+hr*j;
   B=B+r[j]*f[j]*hr;
   double A=Q/(2.0*pi*B);
   q[j]=f[j]*A;
}

for(int j = 0; j < m; j++)