The electron repulsion integrals over Slater−type orbitals with non−integer principal quantum numbers are investigated. These integrals are useful in both non−relativistic and relativistic calculations of many−electron systems. They involve hyper−geometric functions that are practically difficult to compute. Relationships free from hyper−geometric functions for expectation values of Coulomb potential (r21−1) are derived. These relationships are new and show that the complication coming from two−range nature of Laplace expansion for the Coulomb potential is removed. This is achieved by utilizing auxiliary functions represented in finite power series. They serve as essential components in deriving straightforward recurrence relationships for electron repulsion integrals. In the context of computing the expectation values of potentials with arbitrary power, the methodology presented here for evaluation of these integrals forms the initial condition. It is also adapted to multi−center integrals. Program summaryProgram Title: HyperRAFCPC Library link to program files:https://doi.org/10.17632/6pbv2y7s42.1Developer's repository link:https://github.com/abagciphys/HyperRAF.gitLicensing provisions: MITProgramming language:Julia Programming Language [1]Supplementary material: An exploratory variant of the software program written in the Mathematica Programming Language [2].External routines/libraries:Nemo, a computer algebra package for the Julia programming language [3], JRAF, a Julia package for computation of relativistic molecular auxiliary functions [4].Nature of problem: Definite integrals involving higher transcendental functions given by,fmn1(a,b,x)=xm−1e−bxΓ[n,ax], andfmn2(a,b,x)=xm−1e−bxγ[n,ax] are frequently encountered in atomic physics, with the electron repulsion integral being a notable illustration.For exclusive solutions to these integrals, one can refer to Erdélyi's [5] or Gradshteyn and Ryzhik's [6] books. The solutions are obtained by using the series representation of incomplete gamma functions. The result is hyper−geometric functions of the form F12[1,b,c;z], where b=m+n, c=m+1, z=b/(a+b) for fmn1(a,b,x) and c=n+1, z=a/(a+b) for fmn2(a,b,x), respectively. Due to the non−trivial structure of infinite series that are used to define them, the computation for hyper−geometric poses challenges. Convergence of their series strictly depends on the values of parameters. Computational issues such as cancellation or round−off error emerge.Solution method: This research introduces novel bi−directional hyper−radial functions that are used to establish fresh recurrence relationships for electron repulsion integrals, eliminating the dependence on hyper−geometric functions. Dual functionality is inherent in the hyper−radial functions as they offer alternative solutions for definite integrals involving higher transcendental functions. Additionally, they transform the representation of the hyper−geometric functions into finite power series.