A B-spline-type basis set method for the calculation of hydrogen atom in strong magnetic fields in the frame of spheroidal coordinates has been introduced. High accurate energy levels of hydrogen in the magnetic field, with strength ranging from 0 to 1000 a.u., have been obtained. For the ground state, 1s 0 , energies with at least 11 significant digits have been obtained. For the low-lying excited state, 2 p −1 , energies with at least 9 significant digits are obtained. The method has also been applied to the calculation of hydrogen Rydberg states in laboratory magnetic fields. Energy spectra with at least 10 significant digits are presented. A comparison with other results in the literatures has been performed. Our results are comparable to the most accurate one up to date. A possible extension to the cases of parallel and crossed electric and magnetic fields have been discussed.