#include
#include
#include
#include
#include
using namespace std;
# define M_PI           3.14159265358979323846  /* pi */
int main (){
    int n;
    
    /*cin >> n;
    double *r;
    r=new double[n];*/
    double r[50][2];//because n=100
    double xa=0,xb=0;
    double pi=M_PI;
    double ra=0 ,rb=0;

    for(int i=0;i<50;i++){
        
        xa=(rand()%10001)/10000.0;
    xb=(rand()%10001)/10000.0;
        ra=sqrt(-2*log(xa))*cos(2*M_PI*xb);
        rb=sqrt(-2*log(xa))*sin(2*M_PI*xb);

        r[i][0]=ra;
        r[i][1]=rb;

        cout<<"r"<         
        cout << endl;

    }

    return 0;
}

#include
#include
#include
#include
#include
using namespace std;
# define M_PI           3.14159265358979323846  /* pi */
int main (){
    int n;
    
    
double xa;
double xb;
    
    int k[7]={1,2,3,4,5,6,7};
    double fx[7]={0};

    for(int k=0;k<7;k++){

        xa=0.125*(k+1);
    fx[k]=(1/(sqrt(2*M_PI))) * exp((-1/2) * pow(xa,2));
    cout<< fx[k]<

    }

    return 0;
}