// phys.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"
#include <iostream>
#include "math.h"
#include "windows.h"

using namespace std;
#define _QUO_ 5
#include "locale.h"

int _tmain(int argc, _TCHAR* argv[])
{
	double sum(0), mas[5]={0}, mas2[5]={0}, Sigma(0), mid(0), stud(0), xs(0), xp(0), midwidth(0), delta(0);	
	
    HANDLE hStdOut;
	hStdOut=GetStdHandle(STD_OUTPUT_HANDLE);

	setlocale(LC_ALL,"russian");
    for(int i(0); i<_QUO_; i++)
	{
		cout<<"Введите значение "<<i<<" : ";
		cin>>mas[i];				
		sum+=mas[i];
	}
	cout<<endl<<"---------------------------"<<endl;
	cout.setf(ios::left);
	cout.fill(' ');
	cout<<"Сумма значений                      : ";
	cout.width(8);
	cout.precision(3);
	cout<<sum<<endl;
	cout<<"Средняя величина                    : ";
	mid=sum/_QUO_;
	cout.width(8);
	cout.precision(3);	
	cout<<mid<<endl;
    
	/*cout<<"Количество измерений : "<<_QUO_<<endl;
	cout<<"Вероятность          : 0.9"<<endl;*/	
	sum=0;
    for(int i(0); i<_QUO_;i++)
	{
        sum+=pow((mas[i]-mid),2);
	}

	Sigma=sqrt(sum/(_QUO_*(_QUO_-1)));
	cout<<"Сигма                               : "<<Sigma<<endl;
    cout<<endl<<"---------------------------"<<endl;
	cout<<"3*Сигма                             : "<<3*Sigma<<endl;
	 cout<<endl<<"---------------------------"<<endl;
	 
	SetConsoleTextAttribute(hStdOut,FOREGROUND_RED|FOREGROUND_INTENSITY); 
	DWORD nLen=80;
	LPDWORD lpNLen= new DWORD;
    COORD crd;
	crd.X=0;
	PCONSOLE_SCREEN_BUFFER_INFO lpConsoleScreenBufferInfo;
	lpConsoleScreenBufferInfo=new CONSOLE_SCREEN_BUFFER_INFO;	
	for(int i(0); i<_QUO_; i++)
	{	
		if(abs(mas[i]-mid)>3*Sigma){
			crd.Y=i;
		    FillConsoleOutputAttribute(hStdOut, 
				FOREGROUND_RED|FOREGROUND_INTENSITY,
							   nLen,
							   crd,
							   lpNLen
                );	
			GetConsoleScreenBufferInfo(hStdOut,lpConsoleScreenBufferInfo);
			crd.X=_tcslen(_T("Введите значение "))+10;
			SetConsoleCursorPosition(hStdOut,crd);    
			nLen=_tcslen(_T("Значение вне допустимого интервала погрешности!!!"));
			TCHAR *lpBuffer= new TCHAR[nLen+1];
			memset(lpBuffer,0,nLen+1);			
			_tcscpy(lpBuffer,_T("Значение вне допустимого интервала погрешности!!!"));
			WriteConsole(hStdOut,lpBuffer,nLen,lpNLen,NULL);
			crd=lpConsoleScreenBufferInfo->dwCursorPosition;
            SetConsoleCursorPosition(hStdOut,crd);    
			cout<<"Значение вне допустимого интервала погрешности : ";
		    cout.width(9);
	        cout.precision(4);
			cout<<mas[i];			
			cout<<" mas["<<i<<"]"<<endl;
			delta=abs(mas[i]-mid);			
			cout<<"Delta                                          : ";
			cout.width(20);
	        cout.precision(5);
			cout<<delta<<endl;
		}
	}
	delete lpNLen;	
	delete lpConsoleScreenBufferInfo;    
	SetConsoleTextAttribute(hStdOut,FOREGROUND_RED|FOREGROUND_BLUE|FOREGROUND_GREEN);
	
	 cout<<endl<<"---------------------------"<<endl;
	cout.precision(3);
	stud=2.132;
	cout<<"Коэффициент Стьюдента               : "<<stud<<endl;	
	xs=Sigma*stud;
	cout<<"X-случайное                         : "<<xs<<endl;
	cout<<"Введите приборное   Х               : ";
	cin>>xp;
	getchar();
	pow(xp,2);
	pow(xs,2);
    midwidth=sqrt(xp+xs);
	cout<<endl<<"---------------------------"<<endl;
	
	cout<<"Полуширина доверительного интервала : "<<midwidth<<endl;
	cout<<endl<<"---------------------------"<<endl;	
	cout<<"Продолжить вычисления по  М5?(y/n).."<<endl;
	if(getchar()=='y'){
		cout<<endl;
		double D=0, dD(0.0001), h=0,dh(0.0005), g(9.807), dg(0.001), m(0.155), dm(0.0005);		
		cout<<"Введите диаметр блока(метры) : ";
		cin>>D;
		cout<<"Введите высоту               : ";
		cin>>h;
		cout<<endl<<"---------------------------"<<endl;
		double Z=(m*pow(D,2)/4)*(g*pow(mid,2)/(2*h) - 1);
		SetConsoleTextAttribute(hStdOut,FOREGROUND_GREEN);
		cout<<"Момент инерции(I=mD^2(gt^2-2h)/8h) : "<<Z<<endl;
		SetConsoleTextAttribute(hStdOut,FOREGROUND_RED|FOREGROUND_BLUE|FOREGROUND_GREEN);
		double div=sqrt(4*pow(dD/D,2)+pow(dm/m,2)
			           +(4*pow(midwidth/mid,2)+pow(dg/g,2)+pow(dh/h,2))
			           /pow(1-(2*h/(g*pow(mid,2))),2));
        cout<<"dI/I                               : "<<div<<endl;
		SetConsoleTextAttribute(hStdOut,FOREGROUND_GREEN);
		cout<<"=> dI                              : "<<div*Z<<endl;
		SetConsoleTextAttribute(hStdOut,FOREGROUND_RED|FOREGROUND_BLUE|FOREGROUND_GREEN);
        cout<<endl<<"---------------------------"<<endl;
		Z=4*h/(pow(mid,2)*D);
		SetConsoleTextAttribute(hStdOut,FOREGROUND_GREEN);
		cout<<"Угловое ускорение Е(4h/[D*t^2])    : "<<Z<<endl;
		SetConsoleTextAttribute(hStdOut,FOREGROUND_RED|FOREGROUND_BLUE|FOREGROUND_GREEN);
		div=2*midwidth/mid;
		cout<<"dE/E (2dt/midt)                    : "<<div<<endl;
		SetConsoleTextAttribute(hStdOut,FOREGROUND_GREEN);
		cout<<"=> dE                              : "<<div*Z<<endl;		
		SetConsoleTextAttribute(hStdOut,FOREGROUND_RED|FOREGROUND_BLUE|FOREGROUND_GREEN);
		getchar();
		getchar();
	}	
	CloseHandle(hStdOut);
	return 0;
}