квадратный корень от интеграла

#include <cmath>
#include <iostream>
using namespace std;
 
double A, T, phi0;
double u(double x);
 
int main()
{
    system("chcp 1251");
    double h = 0, eps = 0,x,S = 0; 
    cout<<"\tÂâåäèòå ïàðàìåòðû ñèíóñîèäû\n";
    cout<<"Àìïëèòóäà A = ";cin>>A;
    T = -1;
    while(T <= 0)
    {
        cout<<"Ïåðèîä    T = ";
        cin>>T;
        if(T <= 0)
            cout<<"Ïåðèîä êîëåáàíèé äîëæåí áûòü ÷èñëîì áîëüøå 0\n";
    }
    cout<<"Íà÷àëüíàÿ ôàçà phi0(â ðàäèàíàõ) = ";cin>>phi0;
    cout<<"Ââåäèòå òî÷íîñòü èíòåãðèðîâàíèÿ\n";
    cout<<"eps = ";cin>>eps;
    h = eps*(T - 0)/10.0;
    for(x = 0; x < T + h; x = x + h)
        S = S + 0.5*(pow(u(x + h),2.0) + pow(u(x),2.0))*h;
    cout<<"Èíòåãðàë $f(x)^2 ðàâåí = "<<S<<endl;
    cout<<"sqrt($f(x)^2/T)  ðàâíî = "<<sqrt(S/T)<<endl;
    system("pause");
    return 0;
}
 
double u(double x)
{
    double pi = 2*asin(1.0);
    return A*sin((2*pi/T)*x + phi0);
}

#include <cmath>
#include <fstream>
#include <iostream>
using namespace std;
 
int main()
{
    system("chcp 1251");
    int i, T;
    double S = 0;
    double *f = NULL;
    ifstream ifs("data.txt");
    if(!ifs.is_open())
        cout<<"Îøèáêà îòêðûòèÿ ôàéëà äàííûõ\n";
    else
    {
        if(!(ifs>>T))
            cout<<"Ïåðèîä íå ñ÷èòàí\n";
        else
        {
            f = new double [T];
            ifs>>f[0];
            for(i = 1; i < T; i++)
            {
                ifs>>f[i];
                S = S + 0.5*(pow(f[i],2.0) + pow(f[i - 1],2.0))*1.0;
            }
            cout<<"Èíòåãðàë $f(x)^2 ðàâåí = "<<S<<endl;
            cout<<"sqrt($f(x)^2/T)  ðàâíî = "<<sqrt(S/T)<<endl;
            ifs.close();
            delete [] f;
        }
    }
    system("pause");
    return 0;
}