Torus, how?

// http://mathworld.wolfram.com/Torus.html

MBoolean            MCmdCreateTorus::Execute        (MAttributeArray &params)
{
    m_pTheMesh = (MMesh*)MGetScene()->CreateObject ("MMesh", NULL);
    if (!m_pTheMesh)
        return NULL;
    
    MDynamicArray<MVector3D>    FacePoints;

    MFloat  fOuterRadius    = params[0].GetFloat();
    MFloat  fTubeRadius     = params[1].GetFloat();
    MInt    nVertSections   = params[2].GetInteger(); 
    MInt    nHorzSections   = params[3].GetInteger();

    if (nVertSections < 2) nVertSections = 2;
    if (nHorzSections < 2) nHorzSections = 2;

    MFloat fDU = (MMath::PI) / (MFloat)nHorzSections;
    MFloat fDV = (MMath::PI) / (MFloat)nVertSections;
    
    MFloat c = fOuterRadius-(fTubeRadius/2);
    MFloat a = fTubeRadius;

    for (MFloat fV = 0; fV < 2*MMath::PI; fV += fDV)
    {
        for (MFloat fU = 0; fU < 2*MMath::PI; fU += fDU)
        {
            FacePoints.Clear();

            FacePoints.InsertBack (MVector3D ( (c + a*cosf(fV))*cosf(fU+fDU),
                                                (a*sinf(fV)),
                                                (c + a*cosf(fV))*sinf(fU+fDU) ) );
            FacePoints.InsertBack (MVector3D ( (c + a*cosf(fV))*cosf(fU),
                                                (a*sinf(fV)),
                                                (c + a*cosf(fV))*sinf(fU) ) );
            FacePoints.InsertBack (MVector3D ( (c + a*cosf(fV+fDV))*cosf(fU),
                                                (a*sinf(fV+fDV)),
                                                (c + a*cosf(fV+fDV))*sinf(fU) ) );
            FacePoints.InsertBack (MVector3D ( (c + a*cosf(fV+fDV))*cosf(fU+fDU),
                                                (a*sinf(fV+fDV)),
                                                (c + a*cosf(fV+fDV))*sinf(fU+fDU) ) );

            m_pTheMesh->AddFace (FacePoints, TRUE, MMath::EPSILON*2);
        }
    }
    
    

    m_pTheMesh->SetName ("Torus");
    MGetSystem ()->RedrawViews ();

    return TRUE;
}