paraxial pupil calculation function added