NSolve::ifun: Inverse functions are being used by NSolve, so some solutions may not be found; use Reduce for complete solution information.
This means what is says. however: remember that PDE may have functions as arbitrary while ODE constants as arbitrary, and moreso that some ODE you hope to find "a solution" while others you can say "are likely general" while other (by definition or by analysis) are certainly general (all solutions).
My point is before worrying whether your inverse is all solutions ask if you had any chance of finding all solutions.
Use of numbers like ".001" are a bad idea for DSolve - use arbitrary precision numbers (fractions). also - nsolve is not needed for the simple harmonic motion equation (mm can do it).
x''[t] + x[t] == c
is a linear differential equation of order two, non-homogenous
You can find an equation showing all solutions of it by knowing your theorems tell you that is the case . yp+yc, yc can be found using undetermined coefficients.
Your characteristic equation is (m^2+1) will involve imaginary number, which is parlayed into sin/cos function in a sensible way for use with wave equations. You may have a table of simple harmonic equations with the solution outcomes in your book (undamped or damped and unforced or not) so you can plug in constants to known solutions (tenenbaum).