sekwencyjne pi

#include <stdio.h>
#include <stdlib.h>
#include <time.h>

#define PI 3.14159265358979311600

double losuj()
{
    return ( (double)rand() / (double)RAND_MAX );
}


int main(int argc, char *argv[])
{

    double x,y;
    double pi;
    double blad,precyzja;

    int koniec = 0;
    long int kolo = 0;
    long int kwadrat = 0;
    int czas1 = 0;
    int czas2 = 0;

    srand(time(NULL));

    printf("Okresl precyzje obliczen: ");
    scanf("%lf",&precyzja);
    printf("Licze....\n");

    czas1 = time(NULL);

    while(!koniec)
    {
        x = ( (double)rand() / (double)RAND_MAX );
        y = ( (double)rand() / (double)RAND_MAX );

        if( x * x + y * y < 1.0)
            kolo++;
        else
            kwadrat++;

        pi = ( (double)kolo * 4.0 ) / ( (double)kwadrat + (double)kolo );

        blad = fabs( pi - PI );
        if( blad < precyzja )
            koniec = 1;

    }

    czas2 = time(NULL);

    printf("Czas potrzebny na obliczenie: %d\n", czas2 - czas1);
    printf("Dana liczba pi: %1.20f\n", PI);
    printf("obliczona wartość pi: %1.20f dla %ld/%ld trafień do koła z wszystkich losowań\n", pi, kolo, kolo + kwadrat);



    return 0;
}