Как в с++ записать данные в Excel?

#include <iostream>
#include <stdlib.h>
#include <time.h>
#include <cmath>
 
using namespace std;
 
int main()
{
    long Q=5000,Z=1000;
    int k=0, x1=0, x2=0, N=100;
    double m=0, g[1000],s=0.,w=0.,otvet,summa,znach;
    double p=0.;
    srand(time(NULL));
    for(int h=0;h<Z;h++)
    {
        for(int j=0;j<Q;j++)
        {
                x1=0;
                x2=0;
                for(int i = 0; i < N; i++)
                {
                        m = 0.001 * (rand() % 1001);
                        if (0<=m && m<0.125){
                            x2=x2+1; //x
                        }
                        if (0.125<=m && m<0.25){
                            x2=x2-1; //-x
                        }
                        if (0.25<=m && m<0.375){
                            x1=x1+1; //y
                        }
                        if (0.375<=m && m<0.5){
                            x1=x1-1; //-y
                        }
                        if (0.5<=m && m<0.625){
                            x2=x2+1; //x
                            x1=x1+1; //y
                        }
                        if (0.625<=m && m<0.75){
                            x2=x2-1; //-x
                            x1=x1-1; //-y
                        }
                        if (0.75<=m && m<0.875){
                            x2=x2+1; //x
                            x1=x1-1; //-y
                        }
                        if (0.875<=m && m<1){
                            x2=x2-1; //-x
                            x1=x1+1; //y
                        }
                        if (x1==0 && x2==0){
                                k=k+1;
                                break;
                        }
                }
 
        }
        p=(double)((double)k / (double) Q)*(double)100*(double)1.;
        cout<<"Vivod p="<<p<<endl<<endl;
        k=0;
        g[h]=p;
        s=s+p;
    }
    znach=0.;
    summa=0.;
    otvet=0.;
    w= (double)s/ (double)Z*1.; //среднее p
    for (int i=0;i<Z;i++){
        znach=(double) pow((double)((double)w-g[i]),2)*1.;
        summa=(summa+znach)*1.;
    }
    otvet= (double) pow(((double)((double)summa/(double)Z)),0.5)*1.;
    cout<<endl<<endl;
    cout<<"Vivod srednee (p) w="<<w<<endl<<endl;
    cout<<"Vivod sigma="<<otvet<<endl<<endl;
    cout<<"Doveritelnyi interval s doveritelnoy veroyatnostyo 95% ot "<<w-otvet*1.96<<" do "<< w+otvet*1.96 <<endl<<endl;
    cout<<endl<<endl;
    system("pause");
}

hr=OLEMethod(DISPATCH_PROPERTYPUT, NULL, pWApp, L"Visible", 1, x);