A direct method for obtaining discrete relaxation spectra from creep data is proposed. The conversion of creep data to relaxation data is avoided and the discrete relaxation times are freely adjustable. The nonnegative least square method is used to generate nonnegative discrete relaxation strength.