Using the Curve Fitting Toolbox, Matlab supports B-splines, which are a generalization of Bezier curves. A rational B-spline without internal nodes is a Bezier spline.
for example
p = spmak([0 0 0 1 1 1],[1 0;0 1]);
fnplt(p)
(0,0), (1,0), (1,1), (0,1).