A class of multistage and multistep integration methods which can obtain r new values simultaneously at each integration step was developed. Their stability regions were derived and sketched by MATLAB, and their regions are either A-stable or A(a ) -stabl
