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