A new six-parametric, axisymmetric and asymptotically flat exact solution of Einstein-Maxwell field equations having reflection symmetry is presented. It has arbitrary physical parameters of mass, angular momentum, mass-quadrupole moment, current octupole moment, electric charge and magnetic dipole, so it can represent the exterior field of a rotating, deformed, magnetized and charged object; some properties of the closed-form analytic solution such as its multipolar structure, electromagnetic fields and singularities are also presented. In the vacuum case, this analytic solution is matched to some numerical interior solutions representing neutron stars, calculated by Berti & Stergioulas [1], imposing that the multipole moments be the same. As an independent test of accuracy of the solution to describe exterior fields of neutron stars, we present an extensive comparison of the radii of innermost stable circular orbits (ISCOs)