A program to calculate the eigenfunctions of the random phase approximation for two electron systems