A method is presented to simultaneously quantify multiple tissue parameters, including susceptibility, by using a partially RF-spoiled gradient echo sequence. The method successfully obtained whole-brain T1, T2*, PD, B1, and susceptibility maps in 1.1-mm isotropic resolution in 12 minutes at 3T. T1, T2*, PD, and B1 maps were estimated simultaneously by a least-squares fit using twelve 3D source images. The intensity function for the fit was formulated by Bloch simulations. Scan parameters were optimized for T1 and T2* values of the brain at 3T. The susceptibility map was estimated by using six multi-echo images in the source images.