gistfile1.txt

 xydata := [[0,0],[1,9],[4,3],[6,10],[8,5],[9,3]]:
 splajn := Spline(xydata, v):
 polyn := PolynomialInterpolation(xydata, v):
 Psplajn := plot(splajn, v=0..9, color=blue):
 Ppolyn := plot(polyn, v=0..9):
 Pxydata := pointplot(xydata):
 plots[display]([Psplajn, Ppolyn, Pxydata], axes=framed);
 bsplajn := BSplineCurve(xydata, v):
 Pbsplajn := plot(bsplajn, v=0..9, color=cyan):
 plots[display]([Psplajn, Ppolyn, Pbsplajn, Pxydata], axes=framed);
 bspline1 := BSpline(1, v):
 Pbspline1 := plot(bspline1, v=0..4):
 bspline2 := BSpline(2, v):
 Pbspline2 := plot(bspline2, v=0..4, color=green):
 bspline3 := BSpline(3, v):
 Pbspline3 := plot(bspline3, v=0..4, color=blue):
 bspline4 := BSpline(4, v):
 Pbspline4 := plot(bspline4, v=0..4, color=cyan):
 plots[display]([Pbspline1, Pbspline2, Pbspline3, Pbspline4]);
 bspline3s1 := BSpline(3, v, knots=[0, 3, 5, 6]):
 bspline3s2 := BSpline(3, v, knots=[0, 1, 5, 8]):
 bspline3s3 := BSpline(3, v, knots=[0, 3, 5, 8]):
 Pbspline3s1 := plot(bspline3s1, v=0..8):
 Pbspline3s2 := plot(bspline3s2, v=0..8, color=green):
 Pbspline3s3 := plot(bspline3s3, v=0..8, color=blue):
 plots[display]([Pbspline3s1, Pbspline3s2, Pbspline3s3]);