Hyperchaotic system, as an important topic, has become an active research subject in nonlinear science. Over the past two decades, hyperchaotic system between nonlinear systems has been extensively studied. Although many kinds of numerical methods of the system have been announced, simple and efficient methods have always been the direction that scholars strive to pursue. Based on this problem, this paper introduces another novel numerical method to solve a class of hyperchaotic system. Barycentric Lagrange interpolation collocation method is given and illustrated with hyperchaotic system (x˙=ax+dz-yz,y˙=xz-by, 0≤t≤T,z˙=cx-z+xy,w˙=cy-w+xz,) as examples. Numerical simulations are used to verify the effectiveness of the present method.