Запись результатов в файл для набора начальных условий

int main()
{
    FILE *file;
    float D,SH,RL,AM,AN,A,RM,RLL,TCT,PB,PH,TB,NO,N,CP,R,CV,B,XB,DU;
 
    D = 0.05;
    SH = 0.05;
    RL = 0.2;
    TCT = 330;
    PB = 100000;
    PH = 500000;
    TB = 330;
    AM = 0.05;
    NO = 1000;
    N = 7200;
    CP = 1000;
    R = 287;
    CV = 712;
    AN = 1.4;
    A = 0.25;
    B = 600;
    XB = 0.8;
 
    int P = PB;
    int T = TB;
 
    float DWI = (2*3,14159)/N; /*Шаг интегрирования*/
    float WI = 3.14159;
 
    /*СЧИТАЕМ ПАРАМЕТРЫ*/
    float VH = ((3.14159*D*D*SH)/4)*(1*AM);
    float GH = VH*GH*TB;
    float U = CV*GH*TB;
    int i;
     for (i = 1; i <= N; i++) // задаем начальное значение 1, конечное N и задаем шаг цикла - 1.
    {
        float WI = WI + DWI;
        float VH1 = (3.14159*D*D*SH)/4;
        float VM = AM*VH1;
        float V = (VH1/2)*(1-cos(WI))+(RL/4)*(1-cos(2*WI))+VM;
        float VP = fabs((SH/2)*((3.14159*NO)/30)*sin(WI)+(RL/2)*sin(2*WI));
        float RO = P/(R*T); /*ПЛОТНОСТЬ*/
        float RE = (VP*D)/(RM*RO); /*КОЭф РЕЙН*/
        float ANU = A*exp(XB*log10(RE))+B; /*КОЭФ НУСИЛЬДА*/
        float ALFA = (ANU*RL)/D; /*КОЭФ ТЕПЛООТДАЧИ*/
        float FP = (3.14159*D*D)/4;
        float S = V/FP;
        float FT = (FP*2)+3.14159*D*S;
        float CN1 = (3.14159*NO)/30;
        float DTAU = DWI/CN1;
        float DQ = ALFA*FT*(TCT-T)*DTAU;
        float DU = DQ+P*VP+FP*DTAU;
        float U=U+DU;
 
        float P = ((AN-1)*U)/V;
        float T = (P*V)/(GH*R);
 
        scanf("Переменная V = ",V);
        scanf("Переменная DU = ",DU);
        scanf("Переменная P = ",P);
        scanf("Переменная T = ",T);
 
        file = fopen("fileresult.txt", "w");
        fprintf(file, "%s %u %.2f\n", V, DU, P, T);
        fclose(file);
    }
 
    }