A numerical solver 0f 3D Poisson and Nernst-Planck equations, application to ion channel conduction