// - 

#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <windows.h>
#include <string.h>
#include <iostream.h>
#include <time.h>

#define ID_BUTTON_LOSUJ 2001

#define R_FSNEDECORA	0
#define R_CHIKWADRAT	1
#define R_TSTUDENTA		2
#define R_BETA			3

#define SCALE_FSNEDECORA 200
#define SCALE_CHIKWADRAT 50
#define SCALEX_CHIKWADRAT 500
#define IL_KLAS 20

#define MAX_DO_LOSOWANIA 1000

LPCSTR lpszAppName = "Generator liczb losowych o zadanym rozkładzie";
double wartf[1000], wMAX;
double dystr[1000];
double wylosowane[MAX_DO_LOSOWANIA];
unsigned int klasy[IL_KLAS];
unsigned int pfx, pfy;    // wsp punktu na ekranie okreslajacego poczatek ukl wsp dla rozkladu
unsigned int IleLiczb;    // ile liczb ma zostac wylosowanych
double step, scalef;
unsigned int pa, pb;			  // parametry rozkładu
unsigned char cfr, cfg, cfb; // składowe koloru wykresu gestosci
unsigned char chr, chg, chb; // składowe koloru wykresu histogramu

char buf[100];
static char *endptr;

static HWND EditIleLiczb; // ile liczb ma zostac wylosowanych
static HWND ButtonLosuj;  // losuj
static HWND ListBoxWylosowane; // list box z wylosowanymi liczbami
static HWND ComboBoxTyp; // rodzaj rozkladu
static HWND Edit1, Edit2; // parametry rozkładu
static HWND EditRf, EditGf, EditBf; // składowe koloru gestosci
static HWND EditRh, EditGh, EditBh; // składowe koloru histogramu

// zmienne wykorzystywane przez generator
static double uu[97];
static ip=97;
static jp=33;
static double cc=362436.0/16777216.0;
static double cd=7654321.0/16777216.0;
static double cm=16777213.0/16777216.0;

void rozklad_f_snedecora(double m, double n);
void rozklad_chi_kwadrat(double m);
double fun_podc_beta(double m, double n, double x);
double fun_podc_gamma(double m, double x);
double SimpBeta(double m, double n);
double SimpGamma(double m);
void Dystrybuanta(void);
void Losuj(void);
void Sortuj(void);
void GEN_Restart(int i, int j, int k, int l);
double GEN_uni(void);

// Deklarcja procedury okna
LRESULT CALLBACK WndProc(HWND hWnd, UINT message, WPARAM wParam, LPARAM lParam);

// **************************************************************************************
int PASCAL WinMain(HINSTANCE hInstance, HINSTANCE hPrevInstance, LPSTR lpCmdLine, int nCmdShow)
{
    MSG             msg;            // Struktura komunikatu
    WNDCLASS        wc;             // Struktura klasy okna
    HWND            hWnd;           // Uchwyt głównego okna

    // Określenie klasy okna
    wc.style            = CS_HREDRAW | CS_VREDRAW;
    wc.lpfnWndProc      = (WNDPROC) WndProc;
    wc.cbClsExtra       = 0;
    wc.cbWndExtra       = 0;
    wc.hInstance        = hInstance;
    wc.hIcon            = NULL;
    wc.hCursor          = LoadCursor(NULL, IDC_ARROW);

    wc.hbrBackground    = NULL;

    wc.lpszMenuName     = NULL;
    wc.lpszClassName    = lpszAppName;

    // Zarejestrowanie klasy okna
    if(RegisterClass(&wc) == 0) return FALSE;


    // Stworzenie głównego okna aplikacji
    hWnd = CreateWindow(
                lpszAppName,
                lpszAppName,
                WS_OVERLAPPEDWINDOW | WS_CLIPCHILDREN | WS_CLIPSIBLINGS | WS_SYSMENU | WS_MINIMIZEBOX,
				100, 100,
                550, 550,
                NULL,
                NULL,
                hInstance,
                NULL);

    // Wyjście, jeśli nie powiodło się utworzenie okna
    if(hWnd == NULL) return FALSE;

    // Wyświetlenie okna
    ShowWindow(hWnd, SW_MAXIMIZE);
    UpdateWindow(hWnd);

    // Przetwarzanie komunikatów aż do momentu zakończenia działania aplikacji
    while( GetMessage(&msg, NULL, 0, 0))
    {
        TranslateMessage(&msg);
        DispatchMessage(&msg);
    }

    return msg.wParam;
}


                                                                     
// Procedura okna, obsługuje wszystkie komunikaty skierowane do tego okna
LRESULT CALLBACK WndProc(HWND hWnd, UINT message, WPARAM wParam, LPARAM lParam)
{
	HINSTANCE zad;
	HRGN hrgn;
	HBRUSH hBrush;
	HPEN hPen;
	RECT rc, rcHist;
	static HDC hDC;
	double step = 0.05, st=0;
	unsigned int licz=0, i, j;

	switch (message)
    {
        // Utworzenie okna
		case WM_CREATE:
			hDC = GetDC(hWnd);
			zad = ((LPCREATESTRUCT)lParam)->hInstance;

			ComboBoxTyp = CreateWindowEx(0, "COMBOBOX", "F Snedecora", CBS_DROPDOWN | WS_CHILD | WS_VISIBLE, 
									  87, 18, 120, 200, hWnd, NULL, zad, NULL);
			SendMessage(ComboBoxTyp, CB_ADDSTRING, 0, (LPARAM) "F-Snedecora");
			SendMessage(ComboBoxTyp, CB_ADDSTRING, 0, (LPARAM) "chi-kwadrat");
			SendMessage(ComboBoxTyp, CB_ADDSTRING, 0, (LPARAM) "t-Studenta");
			SendMessage(ComboBoxTyp, CB_ADDSTRING, 0, (LPARAM) "beta");
			SendMessage(ComboBoxTyp, CB_SETCURSEL, 0, (LPARAM) 0);

			pa = 8;
			Edit1 = CreateWindowEx(0, "EDIT", "8", WS_CHILD | WS_VISIBLE | WS_BORDER, 
								   306, 60, 20, 20, hWnd, NULL, zad, NULL);

			pb = 8;
			Edit2 = CreateWindowEx(0, "EDIT", "8", WS_CHILD | WS_VISIBLE | WS_BORDER, 
								   372, 60, 20, 20, hWnd, NULL, zad, NULL);

			ButtonLosuj = CreateWindow("BUTTON","Losuj",WS_CHILD|WS_VISIBLE,
			  300, 18, 50, 23, hWnd,(HMENU)ID_BUTTON_LOSUJ, zad, NULL);

			cfr = 200;
			EditRf = CreateWindowEx(0, "EDIT", "200", WS_CHILD | WS_VISIBLE | WS_BORDER, 
								   370, 20, 30, 20, hWnd, NULL, zad, NULL);
			cfg = 0;
			EditGf = CreateWindowEx(0, "EDIT", "0", WS_CHILD | WS_VISIBLE | WS_BORDER, 
								   420, 20, 30, 20, hWnd, NULL, zad, NULL);
			cfb = 0;
			EditBf = CreateWindowEx(0, "EDIT", "0", WS_CHILD | WS_VISIBLE | WS_BORDER, 
								   470, 20, 30, 20, hWnd, NULL, zad, NULL);

			chr = 0;
			EditRh = CreateWindowEx(0, "EDIT", "0", WS_CHILD | WS_VISIBLE | WS_BORDER, 
								   370, 45, 30, 20, hWnd, NULL, zad, NULL);
			chg = 0;
			EditGh = CreateWindowEx(0, "EDIT", "0", WS_CHILD | WS_VISIBLE | WS_BORDER, 
								   420, 45, 30, 20, hWnd, NULL, zad, NULL);
			chb = 200;
			EditBh = CreateWindowEx(0, "EDIT", "200", WS_CHILD | WS_VISIBLE | WS_BORDER, 
								   470, 45, 30, 20, hWnd, NULL, zad, NULL);

			IleLiczb = 100;
			EditIleLiczb = CreateWindowEx(0, "EDIT", "100", WS_CHILD | WS_VISIBLE | WS_BORDER, 
								   306, 60, 30, 20, hWnd, NULL, zad, NULL);

			ListBoxWylosowane = CreateWindowEx(0, "LISTBOX", "100", WS_CHILD | WS_VISIBLE | WS_BORDER | WS_HSCROLL | WS_VSCROLL | LBS_MULTICOLUMN,
							 	   700, 360, 300, 300, hWnd, NULL, zad, NULL);
			
			
			SendMessage(ListBoxWylosowane, LB_SETCOLUMNWIDTH, 50, (LPARAM) NULL);

			GEN_Restart(12,34,56,78);
	
			rozklad_f_snedecora(pa, pb);
//			rozklad_chi_kwadrat(pa);
			Dystrybuanta();
			Losuj();
			Sortuj();	
			for(i=0;i<IleLiczb;i++)
			{
				sprintf(buf, "%.3f", wylosowane[i]);
				SendMessage(ListBoxWylosowane, LB_ADDSTRING, i, (LPARAM) buf);
			}

		break;

        // Okno jest niszczone, więc przeprowadzamy porządki
	    case WM_DESTROY:
		    PostQuitMessage(0);
	    break;

        // Zmianiają się rozmiary okna
		case WM_SIZE:
		break;

        // Funkcja rysująca. 
		case WM_PAINT:
			GetClientRect(hWnd, &rc);
			
			MoveWindow(ComboBoxTyp, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-60, 60, 120, 200, true);
			MoveWindow(Edit1, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-40, 160, 20, 20, true);
			MoveWindow(Edit2, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6+40, 160, 20, 20, true);
			MoveWindow(ListBoxWylosowane, rc.left+30, rc.bottom-290, 630, 270, true);
			MoveWindow(EditRf, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-30, 268, 30, 20, true);
			MoveWindow(EditGf, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6+10, 268, 30, 20, true);
			MoveWindow(EditBf, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6+50, 268, 30, 20, true);
			MoveWindow(EditRh, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-30, 293, 30, 20, true);
			MoveWindow(EditGh, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6+10, 293, 30, 20, true);
			MoveWindow(EditBh, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6+50, 293, 30, 20, true);
			MoveWindow(ButtonLosuj, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6+50, 364, 50, 22, true);
			MoveWindow(EditIleLiczb, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-25, 365, 30, 20, true);
			hBrush = CreateSolidBrush(RGB(200, 200, 200));
			FillRect(hDC, &rc, hBrush);
			DeleteObject(hBrush);

			// narysowanie ramek
			hPen = CreatePen(PS_SOLID, 2, RGB(0, 0, 125));
			SelectObject(hDC, hPen);
			RoundRect(hDC, rc.left+10, rc.top+10, rc.right-(rc.right-rc.left)/3-5, rc.top+370, 10, 10);
			RoundRect(hDC, rc.right-(rc.right-rc.left)/3+5, rc.top+10, rc.right-10, rc.top+100, 10, 10);
			RoundRect(hDC, rc.right-(rc.right-rc.left)/3+5, rc.top+110, rc.right-10, rc.top+200, 10, 10);
			RoundRect(hDC, rc.right-(rc.right-rc.left)/3+5, rc.top+210, rc.right-10, rc.top+340, 10, 10);
			RoundRect(hDC, rc.left+10, rc.top+380, rc.right-(rc.right-rc.left)/3-5, rc.bottom-10, 10, 10);
			RoundRect(hDC, rc.right-(rc.right-rc.left)/3+5, rc.top+350, rc.right-10, rc.top+400, 10, 10);
			RoundRect(hDC, rc.right-(rc.right-rc.left)/3+5, rc.top+410, rc.right-10, rc.bottom-10, 10, 10);
			DeleteObject(hPen);

			TextOut(hDC, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-20, 230, "Kolory:", 7);
			SetTextColor(hDC, RGB(255, 0, 0));
			TextOut(hDC, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-22, 250, "R", 1);
			SetTextColor(hDC, RGB(0, 255, 0));
			TextOut(hDC, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6+20, 250, "G", 1);
			SetTextColor(hDC, RGB(0, 0, 255));
			TextOut(hDC, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6+60, 250, "B", 1);

			SetTextColor(hDC, RGB(cfr, cfg, cfb));
			TextOut(hDC, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-100, 270, "Gęstość:", 8);
			SetTextColor(hDC, RGB(chr, chg, chb));
			TextOut(hDC, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-114, 295, "Histogram:", 10);

			SetTextColor(hDC, RGB(0, 0, 0));
			TextOut(hDC, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-115, 367, "Ilość liczb:",12);
			TextOut(hDC, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-50, 30, "Typ rozkładu: ", 14);
			TextOut(hDC, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-60, 130, "Parametry rozkładu: ", 20);
			TextOut(hDC, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6-70, 160, "m = ", 4);
			TextOut(hDC, ((rc.right-rc.left)/3*2)+(rc.right-rc.left)/6+15, 160, "n = ", 4);
			TextOut(hDC, rc.left+100, rc.top+400, "Wygenerowane liczby (po posortowaniu):", 38);
			TextOut(hDC, rc.right-(rc.right-rc.left)/3+25, rc.top+440, "Program generuje liczby o rozkładach:",36);
			TextOut(hDC, rc.right-(rc.right-rc.left)/3+25, rc.top+460, "   - F-Snedecora",16);
			TextOut(hDC, rc.right-(rc.right-rc.left)/3+25, rc.top+480, "   - chi kwadrat",16);
			TextOut(hDC, rc.right-(rc.right-rc.left)/3+25, rc.top+500, "   - t Studenta",15);
			TextOut(hDC, rc.right-(rc.right-rc.left)/3+25, rc.top+520, "   - beta",9);
			TextOut(hDC, rc.right-(rc.right-rc.left)/3+25, rc.top+560, "Zalecana rozdzielczość 1024x768.",32);
			TextOut(hDC, rc.right-(rc.right-rc.left)/3+25, rc.top+600, "wersja 1.0, autor: Jacek Kuć, WM UŁ",35);
			TextOut(hDC, rc.right-(rc.right-rc.left)/3+25, rc.top+660, "-------------------------- Łódź 2004 --------------------------",63);

			// narysowanie ukladow wspolrzednych
			hPen = CreatePen(PS_SOLID, 2, RGB(0, 0, 0));
			SelectObject(hDC, hPen);
			// gorny uklad wspolrzednych
			MoveToEx(hDC, rc.left+20, rc.bottom-(rc.bottom-rc.top)/2-10, NULL);
			LineTo(hDC, rc.right-(rc.right-rc.left)/3-15, rc.bottom-(rc.bottom-rc.top)/2-10);
			MoveToEx(hDC, rc.left+25, rc.bottom-(rc.bottom-rc.top)/2-5, NULL);
			LineTo(hDC, rc.left+25, rc.top+20);
			// strzalka na osi pionowej
			MoveToEx(hDC, rc.left+25, rc.top+20, NULL);
			LineTo(hDC, rc.left+25-2, rc.top+30);
			MoveToEx(hDC, rc.left+25, rc.top+20, NULL);
			LineTo(hDC, rc.left+25+2, rc.top+30);
			// strzalka na osi poziomej
			MoveToEx(hDC, rc.right-(rc.right-rc.left)/3-15, rc.bottom-(rc.bottom-rc.top)/2-10, NULL);
			LineTo(hDC, rc.right-(rc.right-rc.left)/3-25, rc.bottom-(rc.bottom-rc.top)/2-8);
			MoveToEx(hDC, rc.right-(rc.right-rc.left)/3-15, rc.bottom-(rc.bottom-rc.top)/2-10, NULL);
			LineTo(hDC, rc.right-(rc.right-rc.left)/3-25, rc.bottom-(rc.bottom-rc.top)/2-12);

			pfx = rc.left+25;
			pfy = rc.top+(rc.bottom-rc.left)/2-10;

			st = 0.00;
			step = 5.0/IL_KLAS;

			for(i=0; i<IL_KLAS; i++)
			{
				licz=0;
				j=0;
				for(j=0;j<IleLiczb;j++)
				{
					if(wylosowane[j]>=st && wylosowane[j]<st+step) licz++;
				}
				st+=step;
				klasy[i] = licz;
			}

			// narysowanie histogramu
 			hBrush = CreateSolidBrush(RGB(chr, chg, chb));
			step = (double)(rc.right-rc.left-(rc.right-rc.left)/3-40)/IL_KLAS;
			for(i=0; i<IL_KLAS; i++)
			{
				rcHist.left = pfx+i*step+2;
				rcHist.right = pfx+i*step+step;
				rcHist.top = pfy-klasy[i];
				rcHist.bottom = pfy;
				FillRect(hDC, &rcHist, hBrush);
			}
			DeleteObject(hBrush);

			hPen = CreatePen(PS_SOLID, 2, RGB(cfr, cfg, cfb));
			SelectObject(hDC, hPen);
			step = (double)(rc.right-rc.left-(rc.right-rc.left)/3-40)/1000;
			// znalezienie max w tablicy wartf[]
			wMAX = 0; for(i=0;i<999;i++) if(wartf[i]>wMAX) wMAX=wartf[i];
			scalef = 300;//((rc.bottom-rc.top)/2-50)/wMAX;
			for(i=0; i<999; i++)
			{
				// rozkład
				MoveToEx(hDC, rc.left+25+step*i, rc.bottom-(rc.bottom-rc.top)/2-10-wartf[i]*scalef, NULL);
				LineTo(hDC, rc.left+25+step*(i+1), rc.bottom-(rc.bottom-rc.top)/2-10-wartf[i+1]*scalef);
			}
			DeleteObject(hPen);

			ValidateRect(hWnd, NULL);
		break;

		case WM_COMMAND:
			switch(wParam)
			{
				case ID_BUTTON_LOSUJ:
					SendMessage(EditIleLiczb, WM_GETTEXT, 5, (LPARAM)buf);
    				IleLiczb=strtod(buf, &endptr);
					SendMessage(Edit1, WM_GETTEXT, 5, (LPARAM)buf);
    				pa=strtod(buf, &endptr);
      				SendMessage(Edit2, WM_GETTEXT, 5, (LPARAM)buf);
					pb=strtod(buf, &endptr);
					rozklad_f_snedecora(pa, pb);
//					rozklad_chi_kwadrat(pa);
					Dystrybuanta();
					Losuj();
					Sortuj();	
					SendMessage(ListBoxWylosowane, LB_RESETCONTENT, 0, 0);
					for(i=0;i<IleLiczb;i++)
					{
						sprintf(buf, "%.3f", wylosowane[i]);
						SendMessage(ListBoxWylosowane, LB_ADDSTRING, i, (LPARAM) buf);
					}
      				SendMessage(EditRf, WM_GETTEXT, 5, (LPARAM)buf);
					cfr=strtod(buf, &endptr);
      				SendMessage(EditGf, WM_GETTEXT, 5, (LPARAM)buf);
					cfg=strtod(buf, &endptr);
      				SendMessage(EditBf, WM_GETTEXT, 5, (LPARAM)buf);
					cfb=strtod(buf, &endptr);
      				SendMessage(EditRh, WM_GETTEXT, 5, (LPARAM)buf);
					chr=strtod(buf, &endptr);
      				SendMessage(EditGh, WM_GETTEXT, 5, (LPARAM)buf);
					chg=strtod(buf, &endptr);
      				SendMessage(EditBh, WM_GETTEXT, 5, (LPARAM)buf);
					chb=strtod(buf, &endptr);
	  			    InvalidateRect(hWnd, NULL, TRUE);
				break;
			}
		break;			

	// Domyślna obsługa wszystkich nieprzetwarzanych komunikatów        
		default:   
		return (DefWindowProc(hWnd, message, wParam, lParam));
    }
    return (0L);
}

void rozklad_f_snedecora(double m, double n)
{
	double xx;
	double beta;
	beta = SimpBeta(m/2, n/2);
	for(unsigned int x=0; x<1000; x++)
	{
		xx = (double)x/SCALE_FSNEDECORA;
		wartf[x] = ((pow(m/n, m/2) * pow(xx, (m/2)-1)) / beta)*pow(1+(m/n)*(xx), -(m+n)/2);
	}
}

void rozklad_chi_kwadrat(double m)
{
	double xx;
	double gamma;

	gamma = SimpGamma(m/2);
	for(unsigned int x=0; x<1000; x++)
	{
		xx = (double)(x+1)/SCALE_CHIKWADRAT;
		wartf[x] = (1/(pow(2,m/2)*gamma))*pow(xx,(m/2)-1)*exp(-xx/2);
	}
}


double fun_podc_beta(double m, double n, double x)
{
	return( pow(x, m-1) / pow(1+x, m+n) );
}

double fun_podc_gamma(double m, double x)
{
	return( pow(x, m-1)*exp(-x) );
}

double SimpBeta(double m, double n)
{
	double x=0;
	double s=0, s1, s2, s3;
	double delta=0.05;
	double a, b, c, d;

	d = delta;

	do
	{
		s1 = fun_podc_beta(m, n, x);
		x+=delta*0.5;
		s3 = fun_podc_beta(m, n, x);
		x+=delta*0.5;
		s2 = fun_podc_beta(m, n, x);

		a = 2*(s1+s2-2*s3)/(d*d);
		b = (s2-s1)/d-a*(2*x+d);
		c = s1-a*x*x-b*x;
		
		s+=a*(3*x*x*d+3*x*d*d+d*d*d)/3+b*(2*d*x+d*d)*0.5+c*d;
	} while(x<1000);

	return( s );
}

// funkcja Gamma metoda Simpsona
double SimpGamma(double m)
{
	double x=0;
	double s=0, s1, s2, s3;
	double delta=0.05;
	double a, b, c, d;

	d = delta;

	if(m==0.5) return( 1.7724538509 );
	else
	{
		do
		{
			s1 = fun_podc_gamma(m, x);
			x+=delta*0.5;
			s3 = fun_podc_gamma(m, x);
			x+=delta*0.5;
			s2 = fun_podc_gamma(m, x);

			a = 2*(s1+s2-2*s3)/(d*d);
			b = (s2-s1)/d-a*(2*x+d);
			c = s1-a*x*x-b*x;
		
			s+=a*(3*x*x*d+3*x*d*d+d*d*d)/3+b*(2*d*x+d*d)*0.5+c*d;
		} while(x<1000);

		return( s );
	}
}

void Dystrybuanta()
{
	double s, s1, s2, s3, a, b, c, d=0.00001, x;
	double dziel;

//	switch(typ)
//	{
//		case R_FSNEDECORA: dziel = SCALE_FSNEDECORA; break;
//		case R_CHIKWADRAT: dziel = 10; break;
//		case R_TSTUDENTA:  dziel = 30; break;
//		case R_BETA:       dziel = 300; break;
//	}
	dziel = SCALE_FSNEDECORA;

	FILE *plik;
	plik = fopen("dystr.txt","wt");
	for(unsigned int i=0; i<998; i++)
	{
		s = 0;
		for(unsigned int j=0; j<i; j++)
		{
			s1 = wartf[j];
			s3 = wartf[j+1];
			s2 = wartf[j+2];

			x = (double)j/dziel;

			a = 2*(s1+s2-2*s3)/(d*d);
			b = (s2-s1)/d-a*(2*x+d);
			c = s1-a*x*x-b*x;

			s+=a*(3*x*x*d+3*x*d*d+d*d*d)/3+b*(2*d*x+d*d)*0.5+c*d;
		}
		dystr[i] = (s/dziel)/d;
		fprintf(plik,"%d %.3f\n",i,dystr[i]);
	}
	fclose(plik);
}

void Losuj(void)
{
	double l;
	unsigned int j;	

//	srand( (unsigned)time(NULL) );

	FILE *plik;
	plik = fopen("test.txt","wt");
	for(unsigned int i=0; i<IleLiczb; i++)
	{
		//l = ((double)(rand()%1000))/1000;
		l = GEN_uni();
		
		j = 0;
		while(dystr[j]<l)
		{
			j++;
		}
		if(j!=0)
		{
			if(dystr[j]-l>l-dystr[j-1])
			{
				j--;
			}
		}
		wylosowane[i] = ((double)j)/SCALE_FSNEDECORA;

		fprintf(plik, "%d %.3f %.3f %.3f\n", i, l, dystr[j], wylosowane[i]);
	}
	fclose(plik);
}

void Sortuj(void)
{
	unsigned int i, j;
	double t;

	for(i=0; i<IleLiczb; i++)                                      
	{                                                       
		for(j=IleLiczb-1; j>i; j--)                            
		{                                               
			if(wylosowane[j]<wylosowane[j-1])             
			{                                       
				t = wylosowane[j-1];               
				wylosowane[j-1] = wylosowane[j];      
				wylosowane[j] = t;                 
			}                                       
		}                                               
	}                                                       
}

void GEN_Restart(int i, int j, int k, int l)
{
	int ii, jj, m, wi, wj, wk, wl;
	double s, t;
	wi=i; wj=j; wk=k; wl=l;
	for(ii=0; ii<97; ii++)
	{
		s=0;t=0.5;
		for(jj=1;jj<=24;jj++)
		{
			m=(((wi*wj)%179)*wk)%179;
			wi=wj; wj=wk; wk=m;
			wl=(53*wl+1)%169;
			if((wl*m)%64>=32) s+=t;
			t*=0.6;
		}
		uu[ii]=s;
	}
}

double GEN_uni(void)
{
	double pom;
	pom=uu[ip-1]-uu[jp-1];
	if(pom<0.0) pom+=1;
	uu[ip-1]=pom;
	ip--;
	if(ip==0) ip=97;
	jp--;
	if(jp==0) jp=97;
	cc-=cd;
	if(cc<0.0) cc+=cm;
	pom-=cc;
	if(pom<0.0) pom+=1;
	return(pom);

}