The electric potential near a wall for a multi-species plasma with volume produced negative ions in presence of axially varying magnetic field is studied following an analytical-numerical approach. A constant negative ion source is assumed throughout the plasma volume, along with finite temperature positive ions and Boltzmann electrons. The particles are assumed to be guided by an open magnetic field that has its maximum at the centre, and field strength decreasing towards the walls. The one dimensional (1D) Poisson equation is derived using an analytical approach, and then solved numerically to study the potential profiles. Effect of (a) negative ion production rate, (b) magnetic field profile, and (c) negative ion temperature on the potential profile has been investigated. A potential peak appears near the wall when the negative ion temperature and density are sufficiently high. Also, the presence of negative ions further decreases the potential in the plasma region for a finite Debye Length (λD).