Anonymous

//Paragwgh tyxaiwn arithmwn kanonikhs katanomhs me Monte-Carlo methodo.//To aoristo oloklhrwma ths f(x):Gauss den yparxei kai efarmozw monte-carlo.//Me anakateu8unsh mporeis na pareis to arxeio pou 8es twra. Ta mhnumata pane sto stderr, ta data sto stdout.#include<iostream>#include<cstdlib>#include<cmath>#include<time.h>using namespace std;int main(){     const int a=0,b=1000;  double x,y,f,m,s;  int i=0,j;  cerr<<"\n";  cerr<<" PARGWGH TYXAIWN ARITHMWN KANONIKHS KATANOMHS "<<"\n";  cerr<<"\n";  cerr<<"   Dwse thesh koryfhs:"<<"\n";  cin>>m;  cerr<<"   Dwse diaspora:"<<"\n";  cin>>s;  cerr<<"   Plhthos tyxaiwn arithmwn:"<<"\n";  cin>>j;  cerr<<"\n";     time_t seconds;     time(&seconds);     srand((unsigned int) seconds);  while(i<j){  //tha paraksw tyxaious omoiomorfhs katanomhs sto (m-4s,m+4s)     x=(float)(rand()%(b-a)+a)/1000.0;//x :tyxaio omoiomorfhs katanomhs sto (0,1)    x=(x-0.5)*8*s+m;//to (0,1) metaferetai sto (-0.5,0.5) kai diasteletai kata                     // 8 diaspores ,meta metaferetai gyro apo thn mesh timh m.  //Grafw thn f(x):Gauss synarthsei twn m,s      f=(pow(2.718281828,-(pow((x-m)/s,2))/2)) / (s*(sqrt(2*3.141592654)));  //Zari... y tyxaio omoiomorfhs katanomhs sto (o , f max)    y=(float)(rand()%(b-a)+a)/1000.0;  //y sto (0,1)    y=y/(s*(sqrt(2*3.1415926)));       //y sto (0 , f max)      if(f>=y) {cout<<x<<endl;i++;}  }return 0;}