I decided to write my problem and attach as a separate document so that everything was included and would leave room here to discuss my problem. From what I understand from the equations is that for each value of u I calculate the value of the basis function for each basis function i. Then this is multiplied by the control point. This is done for each order k. The thing that I don't understand is that how can I pull this off if u is the inner most loop? I'm not sure the proper placement of the blending function P(u). Could someone take a look at my code and logic and see if there is another way of doing this or shoot me some pointers? It would be greatly appreciated.