a:
for (int z=1;z<NN;z++)
{
	alfa[z]=F/(1+F*(2-alfa[z-1]));
	beta[z]=(teta[z]+F*beta[z-1])/(1+F*(2-alfa[z-1]));
	nn=StrToFloat(Form1->LabeledEdit14->Text);
	gamma=1-alfa[nn-1]-beta[nn-1];
}
// Вычисление уточненного значения F1:
F1=(gamma*F+K0+(1-teta[nn]))/G;
if (abs(1-F/F1)<0.001)
{
	if (F>250)
	{
		return;}
	else {
		goto b;
	}}
else {
	F=F1;
	goto a;
}
//Обратная прогонка:
b:
teta[nn]=1;
for (int x=nn-1;x==0;x--)
{
	teta[x]=alfa[x]*teta[x+1]+beta[x];
}
timex2=timex1+F*(pow(dy,2))/a;
delta1=delta1-dy;
q2=alam*(teta[1]-teta[0])*(T1-T2)/dy;
goto a;