knots = {0, 0, 0, 0, 1, 1, 1, 1};
GraphicsGrid[
Table[Plot3D[
BSplineBasis[{3, knots}, i, x] BSplineBasis[{3, knots}, j, y], {x,
0, 1}, {y, 0, 1}, PlotRange -> All, Mesh -> 3, Axes -> False,
Boxed -> True,
PlotLabel ->
TraditionalForm[
BSplineBasis[{3, knots}, i, u] BSplineBasis[{3, knots}, j,
v]]], {i, 0, 3}, {j, 0, 3}], ImageSize -> {500, 500}]