The aim of this study is to analyze the numerical solution of magnetohydrodynamic Jeffery fluid past over the upper horizontal parabolic surface with the help of Adam-Milne Predictor Corrector method along with the RK method. Adams predictor-corrector technique is very significant because it improve accuracy of results as compared to using either method alone. The predictor step gives an initial approximation and the corrector step refines this approximation based on the implicit equation. The assumption based on the boundary layer and stagnation point flow of magnetohydrodynamic Jeffery fluid which is past on the melting upper horizontal parabolic surface and the physical aspects are examined with the variable fluid properties. The velocity slip effect on the surface of paraboloid is used to determine its influence on the movement of fluid. The thermal and solutal transfer rates has crucial role in the chemical reactions, climate changes, electronic devices, distillation and separation processes, water and air pollution. Therefore we considered both the thermal and solutal transfer rates with the effects of Joule heating, viscous dissipation, heat source/sink, activation energy and Soret effect. The implementation of all the assumption on the basic conservation laws gives us the governing equation in the form of PDE’s and then the similarity variables are translated these equations into the form of ODE’s. The numerical technique named as ‘Adams-Milne Predictor-Corrector method’ is adopted to solve the numerical solutions. The results are examined in the numerical and graphical forms. The graphical behavior of numerous parameters on the velocity, concentration and temperature regions are analyzed. The numerical findings of skin friction and Nusselt number are also placed here and compared the results with the Bvp5c and Adams-Milne (Predictor-Corrector) method. Graphical Abstract: The slip parameter, ratio of relaxation to retardation parameter, viscosity parameter, Deborah number and Hartmann number drops the velocity for both Newtonian and non-Newtonian cases whereas the velocity increases due to the stretching ratio parameter and melting surface coefficient. The heat source/sink parameter, Eckert number, viscosity parameter, thermal conduction coefficient and Hartmann number. The amplification in concentration region is examined by the consideration of Soret number, thermal diffusion and activation energy, while the reaction rate coefficient drops the concentration.