Catalogs of NEODyS orbital elements are available with or without covariance matrices. These catalogs can be read directly by OrbFit, or you can take the provided FORTRAN code to use in your own programs.
Plot of the NEAs affected by the main secular resonances