The eddy current computation in 3D steady-state problems can be derived from the solution of the vector Helmholtz equation for the magnetic vector potential. The implementation of the finite difference method to the second order vector potential (SOVP) formulation is studied. The method is developed for the numerical solution of 3D eddy current curvilinear problems. The performance is verified for two test cases described, in the cylindrical and spherical coordinate systems. The numerical results, tested against analytical solutions, indicate the technique's robustness as well as considerable savings in computational cost.