float x1[3];
	
x1[0] = x[0]/3.0;
x1[1] = -x[1]/3.0+5;
x1[2] = x[2]/3.0;

	float xx = x1[0]; 
    float y = x1[1]; 
    float z = x1[2]; 

    // head 
	float head1Cent[] = {0, 4.5, 0}; 
    float head = hfEllipsoid(x1, head1Cent, 6, 5.5, 6); 

    // nose 
	float hana1Cent[] = {0, 6.5, 5.9}; 
    float hana = hfEllipsoid(x1, hana1Cent, 0.75, 0.75, 0.75); 

    // face 
	float xface1Cent[] = {0, 4.5, 0.3}; 
    float xface1 = hfEllipsoid(x1, xface1Cent, 6, 5.5, 6); 
    float xface2Cent[] = {0, 3.7, 0}; 
    float xface2 = hfEllCylZ(x1, xface2Cent, 4.7, 4.5); 
    float atama1 = R_int(xface1, xface2); 

    // mouth
    float xkuti1Cent[] = {0, 4.5, 6}; 
    float xkuti1 = hfEllipsoid(x1, xkuti1Cent, 4, 4, 4); 
    float xkuti2 = 4.5-y; 
    float xkuti = R_int(xkuti1, xkuti2); 

    // eyes 
    float eye1Cent[] = {0, 4.5, 0.4}; 
    float eye1 = hfEllipsoid(x1, eye1Cent, 6, 5.5, 6); 
    float eye2Cent[] = {1.5, 7.5, 0}; 
	//float eye2Cent[] = {1.55, 7.5, 0}; 
    float eye2 = hfEllCylZ(x1, eye2Cent, 1.5, 1.4); 
    float eye3 = R_int(eye1, eye2); 

    float eye4Cent[] = {-1.5, 7.5, 0}; 
	//float eye4Cent[] = {-1.55, 7.5, 0}; 
    float eye4 = hfEllCylZ(x1, eye4Cent, 1.5, 1.4); 
    float eye5= R_int(eye1, eye4); 

    float eye6Cent[] = {0.3, 7.5, 5.5}; 
	//float eye6Cent[] = {0.35, 7.5, 5.5}; 
    float eye6 = hfEllipsoid(x1, eye6Cent, 0.3, 0.3, 0.3); 
    float eye7 = R_int(eye3, -eye6); 

    float eye8Cent[] = {-0.3, 7.5, 5.5}; 
	//float eye8Cent[] = {-0.35, 7.5, 5.5}; 
    float eye8 = hfEllipsoid(x1, eye8Cent, 0.3, 0.3, 0.3); 
    float eye9 = R_int(eye5, -eye8); 

    //head final 
    float atama = R_int(R_uni(eye9, R_uni(eye7, R_uni(atama1, R_uni(head, hana)))), -xkuti); 
	
	return atama;