An approximate function approach has been developed using the subsystem synthesis method for real-time multibody vehicle dynamics models. In this approach, instead of solving loop closure constraint equations of the suspension linkage, approximate functions are used. The approximate function represents the functional relationship between dependent coordinates and independent coordinates of the suspension subsystem. This kinematic relationship is also included in the suspension subsystem equations of motion. Different order of polynomial functions are tried to find out the best candidate functions. The proposed method is also compared with the conventional subsystem synthesis method to verify its efficiency and accuracy.