#include  
#include
#include
#include
#include
using namespace std;

const float M_2PIf = 6.283185307179586476925286766559f;
const float M_PIf  = 3.141592653589793238462643383279f;

/* —————–RANDOM_GENERATOR———————————– */

float rand(float min, float max)
{
    return min + (max – min) * rand() / (float)RAND_MAX;
}

struct vektor {
float x,y,z;
};

/* ——————MATRIX_OSZTALY———————————— */

class Matrix {

private:
 
    float tomb[3][3];

public:
 
    void Matrix::init(vektor &a){
    tomb[0][0] = a.x * a.x;
    tomb[0][1] = a.x * a.y;
    tomb[0][2] = a.x * a.z;
    tomb[1][0] = a.y * a.x;
    tomb[1][1] = a.y * a.y;
    tomb[1][2] = a.y * a.z;
    tomb[2][0] = a.z * a.x;
    tomb[2][1] = a.z * a.y;
    tomb[2][2] = a.z * a.z;
    }

Matrix &operator +=(Matrix m){
int k,l;
for(k=0; k<3; k++)
for(l=0; l<3; l++)
  tomb[k][l] += m.tomb[k][l];
return *this;
}

Matrix &operator/=(float c){
int k,l;
for(k=0; k<3; k++)
for(l=0; l<3; l++)
  tomb[k][l] /= c;
return *this;
}

};

/* ——————————————————————– */

int main() {

vektor tar[1000];
int i;

for (i=0;i<1000;i++) {

    float phi  = rand(0.0f, M_2PIf);
    float costheta = rand(-1.0f, 1.0f);
   
    float x = 1.0f * sqrtf(1-costheta*costheta) * cosf(phi);
    float y = 1.0f * sqrtf(1-costheta*costheta) * sinf(phi);
    float z = 1.0f * costheta;

    tar[I].x = x;
    tar[I].y = y;
    tar[I].z = z;
   
    //cout << tar[1].x << " " << tar[1].y << " " << tar[1].z << endl;
   
    }[I][I][I]

//system ("PAUSE");
return 0;
}

std::ostream& operator<<(std::ostream& output, valami) {
    output << "[n";
    for (i=0,i<5,i++) {

    output << " " << valami << " ";

    output << "n";
    }
    output << "]";

    return output;
}