Rigorous high-precision solution of Fuchsian ODEs and applications