Advertisement
Guest User

Untitled

a guest
May 24th, 2015
182
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.87 KB | None | 0 0
  1. #include <stdio.h>
  2. #include <cmath>
  3.  
  4.  
  5. #define NINTERVALOS 100000
  6.  
  7. struct Intervalo
  8. {
  9.  
  10. double a, b;
  11.  
  12. };
  13.  
  14. int main(int argc, char** argv)
  15. {
  16.  
  17. struct Intervalo intervalo [NINTERVALOS];
  18. int n, m, nintervalos;
  19.  
  20. intervalo [0].a=0.0;
  21. intervalo [0].b=1.0;
  22. nintervalos=1;
  23.  
  24.  
  25. do {
  26.  
  27. for(n=nintervalos-1;
  28.  
  29. n>=0 && nintervalos<NINTERVALOS; n--)
  30.  
  31. {
  32. m=nintervalos;
  33.  
  34. intervalo [m].a=intervalo[n].a+2*(
  35. intervalo [n].b
  36. -intervalo[n].a)/3;
  37. intervalo [m].b=intervalo[n].b;
  38.  
  39.  
  40. intervalo [n].a=intervalo[n].a;
  41. intervalo [n].b=intervalo[n].a+(
  42. intervalo [n].b
  43. -intervalo[n].a)/3;
  44.  
  45. nintervalos++;
  46.  
  47.  
  48. } while(nintervalos<NINTERVALOS);
  49.  
  50.  
  51. for(n=0; n<nintervalos; n++)
  52. {
  53.  
  54. printf (" %g\t%g\n", intervalo [n].a,
  55. intervalo [n].b);
  56.  
  57. }
  58.  
  59.  
  60. return 0;
  61. }
  62.  
  63. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement