Курсовая, нужно в виде графика вывести смену электродвижущей силы, все условия и формулы с вводом данных задал, при комп

#define _CRT_SECURE_NO_WARNINGS
#include <stdio.h>
#include <iostream>
#include <cmath>
#define k = 1,380649;
double kc;
int kk;
struct data
{
    double T, t, delt, Tk0, m, m0, n, R0, h, delF;
};
 
double getT(double T, double Tk0, double t, double m);
double getM(double m0, double n, double t double T);
double getR(double h, double T, double t, double R0);
 
void main()
{
    struct data inp[3];
 
    FILE* inp_file = fopen_s("input.txt", "r");
    FILE* outp_file = fopen_s("output.txt", "w");
 
    if (inp_file == 0)
    {
        printf("File not found");
        return 0;
    }
    double E;
    for (int i = 0; i < 3; i++) 
    {
        fscanf(inp_file, "%lf%lf%lf%lf%lf%lf%lf%lf", &inp[i].T, &inp[i].delt, &inp[i].Tk0, &inp[i].m0, &inp[i].n, &inp[i].R0, &inp[i].h, &inp[i].delF);
 
        printf("%d:\n", i + 1);
        for (double t = 0; t <= inp[i].T; t += inp[i].delT) 
        {
            E = sqrt(4 * k * getT(inp[i].T, inp[i].Tk0, t, inp[i].m) * getR(inp[i].h, inp[i].T, t, inp[i].R0)*inp[i].delF);
            fprintf(outp_file, "%.2lf,%lf\n", t, Tk)
            printf("t = %.2lf Tk = %lf\n", t, Tk);
        }
        puts("");
        fputs("\n", outp_file);
    }
 
    return 0;
}
double getT(double T, double Tk0, double t, double m)
{
    if (t >= 0 && t <= T / 4)
    {
        return Tk0 * (1 - exp((-1*getM(double m0, double n, double t double T) * t));
    }
    else
    {
        return Tk0 * (1 - exp((-m) * (T / 4));
    }
}
double getM(double m0, double n, double t double T);
{
    if (t <= T / 4)
    {
        return m0 * (1 + n * t);
    }
    else if (t > T / 4 && t <= ((3 / 4) * T))
    {
        return m0 * (1 + (n * (T / 4)));
    }
    else
    {
        return m0 * (1 + (n * T / 4) - n * (t - 3 * T / 4));
    }
}
double getR(double h, double T, double t, double R0);
{
    if (t >= 0 && <= T / 3)
    {
        return R0 * (1 + h * t);
    }
    else if (t > T / 3 && t <= ((2 / 3) * T)
    {
        return R0(1+h*T/3);
    }
    else
    {
        return R0 * (1 + h * T / 3) * exp(-h * (t - (T * (2 / 3)));
    }
}

double getT(double T, double Tk0, double t, double m)
{
if (t >= 0 && t <= T / 4)
{
return Tk0 * (1 - exp((-1*getM(double m0, double n, double t double T) * t));
}
else
{
return Tk0 * (1 - exp((-m) * (T / 4));
}
}
double getM(double m0, double n, double t double T);
{
if (t <= T / 4)
{
return m0 * (1 + n * t);
}
else if (t > T / 4 && t <= ((3 / 4) * T))
{
return m0 * (1 + (n * (T / 4)));
}
else
{
return m0 * (1 + (n * T / 4) - n * (t - 3 * T / 4));
}
}