Higher-order accurate finite-difference schemes for numerical solution of Burgers’ equation which often arises in mathematical modelling used to solve problems in fluid dynamics is presented. We also consider shape-preserving properties of the integro-cubic spline. The sufficient conditions for convex, monotone, or positive approximation are given. Finally, we give an application of the integro cubic spline to approximate the numerical solution of Burgers' PDE.