The Cramér-Rao Lower Bound (CRLB) is a metric for optimizing quantitative sequences that requires an analytical expression for the signal. The CRLB for Magnetic Resonance Fingerprinting (MRF) has a complex formulation that makes it difficult to account for system imperfections or relevant signal contributions such as diffusion. We apply automatic differentiation to Bloch simulations and choose flip angles and repetition times that optimize the CRLB of the MRF sequence without deriving an explicit analytical expression for the signal. This method is computationally efficient and can easily be extended to include B0 and B1 inhomogeneities. Results are validated with in-vivo measurements.