diff options
Diffstat (limited to 'rqmc.cpp')
| -rw-r--r-- | rqmc.cpp | 37 |
1 files changed, 19 insertions, 18 deletions
@@ -5,6 +5,7 @@ double frac_part(double x){ return x - floor(x); } + double mean_rqmc(int N, double X) { sobol s(1); double sum = 0; @@ -14,24 +15,24 @@ double mean_rqmc(int N, double X) { return sum/N; } -int main() { - init_alea(0); - int I=100; - int N= 10000; - uniform U; - double m = 0; - double s = 0; - double temp; +//~ int main() { + //~ init_alea(0); + //~ int I=100; + //~ int N= 10000; + //~ uniform U; + //~ double m = 0; + //~ double s = 0; + //~ double temp; - for(int i=0;i<I;i++){ - temp = mean_rqmc(N,U()); - m+=temp; - s+=temp*temp; - } - m = m/I; - s = s/I - m*m; + //~ for(int i=0;i<I;i++){ + //~ temp = mean_rqmc(N,U()); + //~ m+=temp; + //~ s+=temp*temp; + //~ } + //~ m = m/I; + //~ s = s/I - m*m; - std::cout<<"espérance "<<m<<" taille de l'IC "<<sqrt(s)*1.96/10<<std::endl; + //~ std::cout<<"espérance "<<m<<" taille de l'IC "<<sqrt(s)*1.96/10<<std::endl; - return 0; -} + //~ return 0; +//~ } |
