Внедрение одной программы в другую

void __fastcall TForm1::Button3Click(TObject *Sender)
{
 
//изначальный код за считывание данных из файла и дальнейшего их использования( сохранение в память как я понял)
int nist, j;
if (OpenDialog2->Execute()) {
    ifstream input(OpenDialog2->FileName.c_str());
    vec FR, FF; /* для задания распределения плотности потока по радиусу мишени */
    double value1;
   input >> value1;
   while (!input.eof())  {
      FR.push_back(value1);
      input >> value1;
      FF.push_back(value1);
      input >> value1;
   }
   nist = FF.size();
 
//Построение графика, вот здесь  я пробовал менять FR на rrr[j] , FF на vsput[j]
   XYPlot1->XYPlot(FR.begin(), FF.begin(), nist, 0, clGreen);
   
 double r;
 
 
 
 
//добавлен мой код
 FILE*fp;
   if (InputOpen->Execute()) {
   fp = fopen(InputOpen->FileName.c_str(),"rt");
     if(!fp)
 
         {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);
       }
 
   rrr.resize(m);  // координаты ячеек вдоль поверхности мишени
   q.resize(m); // плотность мощности
   vsput.resize(m); // скорость распыления при бомбардировке ионов рабочего газа
 
 
 
    // Ввод данных пространственного распределения плотности мощности вдоль поверхности мишени
    // для задания распределения плотности тока по r
   hr=R/m;  //разделение подложки на м участков
   for (j=0; j<m; j++)
   {
   rrr[j]=hr/2+hr*j;
   AnsiString S_help;
   S_help=AnsiString(S)+".txt";
   ifstream input(S_help.c_str());
   }
 
// Вот здесь скорей всего неверно, поскольку я должен вводить еще один другой файл с данными в качестве значении FR FF, но //это я смогу реализовать сам
   vec f;
   f.resize(m);
   for (j=0; j<m; j++)
      f[j]=approxd(rrr[j],nist,FR.begin(),FF.begin());
   double B=0;
   for (j=0; j<m; j++)
   {
      B=B+rrr[j]*f[j]*hr;
   }
   double A=Q/(2.0*pi*B);
   for (j=0; j<m; j++)
{
     q[j]=f[j]*A;
}
FILE* file = fopen("test2.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/(ro * Na/atom);
      rrr[j]=hr/2+hr*j;
      fprintf(file, "%.3E %.3E\n", rrr[j], vsput[j]);
 
   }
fclose(file);
}
 
//здесь заканчивается мой код
 
 
//расчет некой величины, которую дальше нужно будет сверять( вот здесь как раз у меня используется FF FR
 for (int i=0; i<N; i++)
{
    r=sqrt(pow((X0-M[i].xm),2)+pow((Y0-M[i].ym),2));
    M[i].F=M[i].F+approxd(r,nist,FR.begin(),FF.begin())*nat; /* шт./м**2*с */
}
}
// Проверка
  Edit6->Text=AnsiString(FormatFloat("0.000E+00",M[1].xm));
  Edit7->Text=AnsiString(FormatFloat("0.000E+00",M[1].ym));
  Edit8->Text=AnsiString(FormatFloat("0.000E+00",M[1].F));
}