Программа виснет при выполнении

#include <math.h>
#include <vcl.h>
#pragma hdrstop
 
#include "Main.h"
#include "Result.h"
#include "Error.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
float W;
float A;
float F;
float T;
float tay;
float n;
float K;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
vA1->Clear();
vA2->Clear();
vA3->Clear();
vA4->Clear();
vA5->Clear();
 
vF1->Clear();
vF2->Clear();
vF3->Clear();
vF4->Clear();
vF5->Clear();
 
vW1->Clear();
vW2->Clear();
vW3->Clear();
vW4->Clear();
vW5->Clear();
 
vtay->Clear();
vZ->Clear();
vn->Clear();
 
Edit1->Clear();
Edit2->Clear();
Edit3->Clear();
Edit4->Clear();
 
 
}
//---------------------------------------------------------------------------
 
float Amod();
float Fmod();
float Tr();
float Kr();
 
void __fastcall TForm1::Button1Click(TObject *Sender)
{
float W1=0.0129; //StrToFloat(vW1->Text);
float W2=0.2719; //StrToFloat(vW2->Text);
float W3=0.5309; //StrToFloat(vW3->Text);
float W4=0.7899; //StrToFloat(vW4->Text);
float W5=1.0489; //StrToFloat(vW5->Text);
float A1=1.9992; //StrToFloat(vA1->Text);
float A2=1.7131; //StrToFloat(vA2->Text);
float A3=1.2207; //StrToFloat(vA3->Text);
float A4=0.8288; //StrToFloat(vA4->Text);
float A5=0.5728; //StrToFloat(vA5->Text);
float F11=-1.2697; //StrToFloat(vF1->Text);
float F21=-50.5946; //StrToFloat(vF2->Text);
float F31=-89.1217; //StrToFloat(vF3->Text);
float F41=-117.427; //StrToFloat(vF4->Text);
float F51=-138.649; //StrToFloat(vF5->Text);
float n=2; //StrToFloat(vn->Text);
float tay=StrToFloat(vtay->Text);
float Z=StrToFloat(vZ->Text);
 
 
float F1=F11*3.14/180;
float F2=F21*3.14/180;
float F3=F31*3.14/180;
float F4=F41*3.14/180;
float F5=F51*3.14/180;
 
W=W1;
A=A1;
F=F1;
 
while (abs (tay*W>=(1/2)*F))
        tay=tay/2;
 
T=Tr();
K=Kr();
 
 
W=W2;
A=A2;
F=F2;
 
float Am=Amod();
float Fm=Fmod();
 
if (abs ((A-Am)/A)<=Z)
        if (abs ((F-Fm)/F)>=Z)
                return;
        else
                {
                n=2*(log(K/A)/log(pow(pow(T,2)*pow(W,2)+1,n/2)));
                tay=(-F-n*atan(T*W))/W;
 
                W=W3;
                A=A3;
                F=F3;
 
                Am=Amod();
                Fm=Fmod();
 
                if (abs ((A-Am)/A)<=Z)
                        if (abs ((F-Fm)/F)>=Z)
                                return;
                        else
                                {
                                T=Tr();
                                K=Kr();
 
                                W=W4;
                                A=A4;
                                F=F4;
 
                                Am=Amod();
                                Fm=Fmod();
 
                                if (abs ((A-Am)/A)<=Z)
                                        if (abs ((F-Fm)/F)>=Z)
                                                return;
                                        else
                                                {
                                                n=2*(log(K/A)/log(pow(pow(T,2)*pow(W,2)+1,n/2)));
                                                tay=(-F-n*atan(T*W))/W;
 
                                                W=W5;
                                                A=A5;
                                                F=F5;
 
                                                Am=Amod();
                                                Fm=Fmod();
 
                                                 if (abs ((A-Am)/A)<=Z)
                                                        if (abs ((F-Fm)/F)>=Z)
                                                                return;
                                                        else
                                                        {                                                     
                                                        Edit5->Text="сменить структуру модели";
                                                        }
                                                }
                                }
                }
 
 
Edit1->Text=n;
Edit2->Text=T;
Edit3->Text=K;
Edit4->Text=n;
}
 
/*вычисление T*/
float Tr()
        {
        double T=tan(-F-tay*W)/n/W;
        return T;
        }
 
/*вычисление K*/
float Kr()
        {
        double K=A*pow(pow(T,2)*pow(W,2)+1,n/2);
        return K;
        }
 
/*вычисление модуля модели*/
float Amod()
        {
        double Am;
         Am=K/pow(pow(T,2)*pow(W,2)+1,n/2);
         return Am;
        }
/*вычисление фазы модели*/
float Fmod()
        {
        double Fm;
        Fm=-n*(atan(T*W)-tay*W);
        return Fm;
        }