vm_labs / lr_6 / LR_6

/*
  ” ©«: Lr_6.cpp


  Ђўв®а: “ᥭЄ® Ђ.‚.
  ђҐ¤. 1.0  08.10.06
*/

#include <iostream.h>
#include <math.h>
#include <conio.h>
#include <fstream.h>
#include "methods.h"

// Global variables
double	F(double);	// Function
double	F1(double);	// Derivative

double	Fi(double);	// x=Fi(x)
double	Fi1(double);	// Fi(x) derivative

double	delta=0.1;	// function calculating error

double	M=2.5,		// Max F1(x), x in [0.2;0.6]
	m=2.2;		// Min F1(x), x in [0.2;0.6]

ofstream	f_out_1("out1.txt"),	// Output data files
		f_out_2("out2.txt"),
		f_out_3("out3.txt");


double	eps=0.1,	// output data error
	eps_theor=0;	// Theoretical eps value



void main()
{
 clrscr();

 double	func_val,	// value of function
	x,		// value of root
	_F,_Fi,_Fi1;
 int	n;		// number of iterations

 // file opening

 if (f_out_1.fail() || f_out_2.fail() || f_out_3.fail())
 {
	cout<<"Impossible to open output file!"<<endl;
	exit(1);
 }



 f_out_1.setf(ios::fixed);
 f_out_1.setf(ios::showpoint);

 f_out_2.setf(ios::fixed);
 f_out_2.setf(ios::showpoint);
 f_out_2.precision(9);

 /* Calculating & table printing */
 delta=0.1;

 for (int j=1; j<=6; j++)
 {
	f_out_1.precision(j);
	f_out_1<<"Delta = "<<delta<<endl;

	eps=0.1;
	f_out_1<<"Eps\t\t"<<"X\t\t"<<"Fi(x)\t\t"<<"Fi'(x)\t\t"<<"F(x)\t\t"<<"N\t"<<"Eps_Theor\t"<<"Condit\n";
	for (int i=1; i<=6; i++)
	{
		x=0.6;
		f_out_1.precision(i);
		x=ITER(x,eps,n,f_out_2);
		_F=F(x); _Fi=Fi(x);_Fi1=Fi1(x);
		f_out_1<<eps<<"\t"; if (i<6) {f_out_1<<"\t";};
		f_out_1<<x<<"\t";   if (i<6) {f_out_1<<"\t";};
		f_out_1.precision(j);
		f_out_1<<_Fi<<"\t"; if ( (j<6) && !((j==5)&&(_Fi<0)) ) {f_out_1<<"\t";};
		f_out_1<<_Fi1<<"\t"; if ( (j<6) && !((j==5)&&(_Fi1<0)) ) {f_out_1<<"\t";};
		f_out_1<<_F<<"\t"; if ( (j<6) && !((j==5)&&(_F<0)) ) {f_out_1<<"\t";};
		f_out_1<<n<<"\t";
		eps_theor=delta/fabs(1/(1-Fi1(x)));
		f_out_1.precision(6);
		f_out_1<<eps_theor<<"\t";
		if (eps>=eps_theor) {f_out_1<<"good\n";}
		else	{f_out_1<<"bad\n";}

		eps = eps/10;
	}

	f_out_1<<endl<<endl;
	delta = delta/10;
 }


 f_out_3.setf(ios::fixed);
 f_out_3.setf(ios::showpoint);
 f_out_3<<"X\t"<<"F(x)\t\tF'(x)\t\tFi(x)\t\tFi'(x)\n";
 for (x=0;x<1;x=x+0.01)
 {
  f_out_3.precision(2);
  f_out_3<<x<<"\t";
  f_out_3.precision(9);
  f_out_3<<F(x)<<"\t"<<F1(x)<<"\t"<<Fi(x)<<"\t"<<Fi1(x)<<endl;
 }


 f_out_1.close();
 f_out_2.close();
 f_out_3.close();

 cout<<"Program finished. Press <Enter>...";
 getch();
 return;
}

// Function definition
double	F(double x)
{
	extern double delta;
	double s;

	s = asin(2*x/(1+x*x))-exp(-x*x);
	return Round(s,delta);
}

// Derivative definition
double	F1(double x)
{
	double	s;

	s=(2-4*x*x/(1+x*x))/((1+x*x)*sqrt(1-4*x*x/pow((1+x*x),2))) + 2*x*exp(-x*x);
	return Round(s,delta);
}


// x=Fi(x)
double	Fi(double x)
{
	extern	double	M,m,delta;
	double	s;

	s=x-(2/(M+m))*F(x);

	return	Round(s,delta);

}

// First derivative of Fi(x)
double	Fi1(double x)
{
	extern	double	M,m,delta;
	double	s;

	s=1-(2/(M+m))*F1(x);

	return	Round(s,delta);
}