Проверьте пожалуйста

#include "stdafx.h"
#include <stdlib.h>
#include <iostream>
#include <time.h>
#include <fstream>
#include <string>
#include <cmath>
using namespace std;
 
 
double func(int a, int b, int d, double x)
{
    return 1.0 * a * cos(1.0 * b*x) + pow(x, d); 
}
 
double method_of_trapezium(int a, int b, int d, double x1, double x2, double dx)
{
    double x = x1;
    double sum = (func(a, b, d, x1 + dx) + func(a, b, d, x1)) * dx / 2.0;
    while(x < x2)
    {
        sum += (func(a, b, d, x + dx) + func(a, b, d, x))*dx / 2.0;
        x += dx;
    }
 
    cout<<"result: "<<sum<<endl;
    return sum;
}
 
 
void writeToFile(double val)
{
    char * fileName = "output.txt";
    FILE * file = fopen(fileName, "w");
    if (file) // если есть доступ к файлу,
    {
        fprintf(file,"%f", val);
    }
    //Закрываем файл
    fclose(file);
 
}
 
void main(int argc, char argv[])
{
 
    clock_t t;
    t = clock();
 
    int a,b,d;
    double x1,x2;
 
    double dx = 0.001;
 
 
    //Открываем файл для чтения
    ifstream infile;
    infile.open ("file.txt");
    int i =0;
    //Считываем строки из файла до тех пор, пока не дошли до конца файла
    while(infile.eof() == false)
    {
        
        //Получаем строку из файла
        string line = "";
        getline(infile, line);
 
        if (line.find("a=") != string::npos)
        {
            line.erase(0,2);
            a = atoi(line.c_str());
        }
 
        if (line.find("b=") != string::npos)
        {
            line.erase(0,2);
            b = atoi(line.c_str());
        }
 
        if (line.find("d=") != string::npos)
        {
            line.erase(0,2);
            d = atoi(line.c_str());
        }
 
        if (line.find("x1=") != string::npos)
        {
            line.erase(0,3);
            x1 = atof(line.c_str());
        }
 
        if (line.find("x2=") != string::npos)
        {
            line.erase(0,3);
            x2 = atof(line.c_str());
        }
        if (line.find("dx=") != string::npos)
        {
            line.erase(0,3);
            dx = atof(line.c_str());
        }
    }
 
    cout<<"a="<<a<<endl;
    cout<<"b="<<b<<endl;
    cout<<"d="<<d<<endl;
    cout<<"x1="<<x1<<endl;
    cout<<"x2="<<x2<<endl;
    cout<<"dx="<<dx<<endl;
 
    double result = method_of_trapezium(a,b,d,x1,x2,dx);
 
    writeToFile(result);
 
    t = clock() - t;
 
    cout<<"File created. Time: ";
    printf("%f", (double)t/CLOCKS_PER_SEC);
    cout<<endl;
}