B-spline curve, especially cubic B-spline curve, has been widely used in CAD area. Sometimes we need to compute a B-spline curve’s control points according to points on it. The problem of computing a cubic B-spline curve’s control points is equivalent to that of solving a linear equation, which can be solved by LU decomposition or Gauss elimination. In this paper, a better algorithm, which makes use of the characteristic of the equation’s coefficient matrix, will be presented.