Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include <iostream>
- #include <math.h>
- using namespace std;
- double f1(double x){
- double res = 1/pow(3*atan(x), 1.0/3);
- return res;
- }
- double f2(double x){
- double res = atan(x) - 1/ (3*pow(x,3));
- return res;
- }
- int main()
- {
- setlocale(LC_ALL, "rus");
- cout << "1 метод" << endl;
- double x = 0.76;
- double eps = 1e-3;
- double prev_res = f1(x);
- double curr_res;
- while(abs(prev_res - f1(prev_res)) > eps ){
- cout<<prev_res<<endl;
- prev_res = f1(prev_res);
- curr_res = f1(prev_res);
- }
- cout<<curr_res<<' ';
- cout << "\n\n" << "2 метод" << endl;
- double a = 0, b = 1;
- double prev_x = (a + b) / 2;
- if((f2(a) * f2(prev_x)) < 0){
- b = prev_x;
- }
- else{
- a = prev_x;
- }
- double curr_x = (a + b) / 2;
- while(abs(prev_x - curr_x) > eps){
- prev_x = curr_x;
- if(f2(a) * f2(prev_x) < 0){
- b = prev_x;
- }
- else{
- a = prev_x;
- }
- curr_x = (a + b) / 2;
- cout<<prev_x<<' ';
- cout<<f2(a) * f2(prev_x) << ' ';
- cout << f2(b) * f2(prev_x) << endl;
- }
- cout<<curr_x<<' ';
- cout<<f2(a) * f2(curr_x) << ' ';
- cout << f2(b) * f2(curr_x) << endl;
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement