In the previous papers, 25, we have obtained analytical results for global parameters, vital for the description of the chemical reaction: Δ𝑁, Δ𝐸 and Δ𝜇. They were calculated within the alchemical approach, up to the second order for both open and closed systems, using the operating approximations to the local softness (eq. 3) and its derivative over N (eq. 5). This has motivated the computational results of this present paper. With the use of 𝑠(𝐫) and