Least squares methods for solving partial differential equations by using Bézier control points

We discuss the least squares methods for solving partial differential equations on arbitrary polygon domain by using Bézier control points. It uses triangular Bézier patches of degree n   with CkCk continuity to approximate the exact solution of partial differential equations. Numerical examples are provided to illustrate the proposed method is flexible.