Advertisement
Guest User

Untitled

a guest
Nov 21st, 2017
74
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.53 KB | None | 0 0
  1. void BoxMiller(float mi , float sigma)
  2. {
  3. std::ofstream myfile("out.txt", ios::out);
  4.  
  5. const float max = 1.0f;
  6.  
  7. for (int i = 0; i < 1000; i++)
  8. {
  9. float pi = 3.14f;
  10. float z1 = float(rand()) / RAND_MAX;
  11. float teta = 2.0*pi * float(rand()) / RAND_MAX;
  12.  
  13. float r2 = (-2 * (log(z1)));
  14. float r = sqrt(r2);
  15.  
  16. float y1 = r*(cos(teta));
  17. float y2 = r*(sin(teta));
  18.  
  19.  
  20.  
  21. float temp = (y1 * sigma) + mi;
  22.  
  23. printf("%6.5f \n", temp);
  24.  
  25. myfile << std::to_string(temp) << "\n";
  26. }
  27.  
  28. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement