Advertisement
Guest User

Untitled

a guest
Jan 28th, 2020
104
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.63 KB | None | 0 0
  1. #include <iostream>
  2. #include <bits/stdc++.h>
  3. using namespace std;
  4.  
  5. double FofX(double input) {
  6. return 1/input - cos(input);
  7. }
  8.  
  9. double root(double ak, double bk) {
  10. return (ak * FofX(bk) - bk * FofX(ak)) / (FofX(bk) - FofX(ak));
  11. }
  12.  
  13. int main() {
  14. int k = 0;
  15. double a = 4.0;
  16. double b = 6.0;
  17. double x = 0.0;
  18. cout << "k\t\tak\t\tbk\t\txk\t\t" << endl;
  19.  
  20. for (k = 0; k <= 5; k++) {
  21. x = root(a, b);
  22.  
  23. cout << k << "\t\t" << a << "\t\t" << b << "\t\t" << x << endl;
  24.  
  25. if (FofX(x) * FofX(a) < 0) {
  26. b = x;
  27. } else {
  28. a = x;
  29. }
  30. }
  31.  
  32. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement