lab6 / main

яю#include "pch.h"

#include <iostream>

#include <stdio.h>

#include <math.h>

#include <stdlib.h>

#include "METHODS.H"



double delta, c, d;



double Round(double X, double Delta){

	if (delta <= 1E-9) { puts("No okruglenue--\n"); exit(1); }

	if (X > 0.0) return (delta*(long((X / delta) + 0.5)));

	else    return (delta*(long((X / delta) - 0.5)));

}



double F(double x) { // f(x ) = (1+cos(x))/(3 - sin(x))

	double s;

	long int S;

	s = (1 + cos(x)) / (3 - sin(x)) ;





	if (s / delta < 0)

		S = s / delta - 0.5;

	else

		S = s / delta + 0.5;

	s = S * delta;



	s = Round(s, delta);

	return s;

}





double F1(double x) {

	double s;

	long int S;

	s = (cos(x)*cos(x) + cos(x) + 3 * sin(x) - 9) / (3 - sin(x)) / (3 - sin(x));



	if (s / delta < 0)

		S = s / delta - 0.5;

	else

		S = s / delta + 0.5;

	s = S * delta;



	s = Round(s, delta);

	return s;

}



double ITER(double X0, double Eps, int &N){

	if (Eps <= 0.0) { puts("525@=>5 7040=85 B>G=>AB8\n"); exit(1); }

	double X1 = F(X0);

	double X2 = F(X1);

	N = 2;

	while ((X1 - X2)*(X1 - X2) > fabs((2 * X1 - X0 - X2)*Eps)){

		X0 = X1;

		X1 = X2;

		X2 = F(X1);

		N++;

	}

	return(X2);

}



int main() {

	int k;

	long int s;

	double a, b, eps, x;

	double F(double);



	a = -1.0;

	b = 2.0;



	printf("delta:");

	std::cin >> delta;



	eps = 0.1;

	x = ITER(b, eps, k);

	printf("x=%lf k=%d\n", x, k);

	eps = 0.01;

	x = ITER(b, eps, k);

	printf("x=%lf k=%d\n", x, k);

	eps = 0.001;

	x = ITER(b, eps, k);

	printf("x=%lf k=%d\n", x, k);

	eps = 0.0001;

	x = ITER(b, eps, k);

	printf("x=%lf k=%d\n", x, k);

	eps = 0.00001;

	x = ITER(b, eps, k);

	printf("x=%lf k=%d\n", x, k);

	eps = 0.000001;

	x = ITER(b, eps, k);

	printf("x=%lf k=%d\n", x, k);



	return 0;

}