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];

Dv = vVector * M; // This is a 1x4 mult 4x4 per X,Y,Z,W

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

uFloat = uInt * Ustep;
uVector = [uFloat³ uFloat² uFloat 1];

Vtx= Dv * uVector;//This is a 1x4 mult 1x4 per X,Y,Z,W output(Vtx);

}

}

}