#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;
}