Vstep = 1 / (NumberOfVVertices-1);
Ustep = 1 / (NumberOfUVertices-1);

M = B * P * B^T;

for(unsigned int vInt=0; vInt < NumberOfVVertices; vInt++)

{

vFloat = vInt * Vstep;
vVector = [vFloat³ vFloat² vFloat 1];

D = (vVector * M) * [Ustep³ Ustep² Ustep 1];

Vtx = ???;
duVtx = ???;
dduVtx = ???;
ddduVtx = ???;

for(unsigned int uInt=0; uInt < NumberOfUVertices; uInt++)
{

output(Vtx);

Vtx += duVtx;
duVtx += dduVtx;
dduVtx += ddduVtx;

}

}

}