realspin

    int n=0;
    float1 h=6.626e-34;
    float3 Bv=float3(0,0,0);
    for(int y=0;y<900;y++)
    for(int x=0;x<3600;x++)
    {
        float1 a1=(float1)x*M_PI*2/3600.0;
        float1 a2=(float1)y*M_PI*2/3600.0;
        Bv=Bv+float3(cos(a1)*cos(a2),sin(a1)*cos(a2),sin(a2))*h;
        n++;
    }
    Bv=Bv/n;
    printf("%Lf  \n",atanl(h/Bv.z)*180.0/M_PI);
    printf("%Lf  \n",((h/Bv.z)/(M_PI/2)));