AbstractExtraction of crystallization kinetics, the fundamental information which governs the performance of a crystallization process, is characterized by experimental difficulties in MSMPR and mathematical difficulties in batch crystallization. Here, a rigorous approach is taken to estimate kinetic parameters from a batch crystallizer. The two step LaxWendorff technique is adapted for the solution of the population balance in a batch crystallizer and an iterative self‐correcting least squares algorithm is implemented for the estimation of the kinetic parameters. The need for multi‐response estimation as opposed to single‐response from terminal CSD is demonstrated. The kinetics extracted are average ones, representing primary and secondary nucleation kinetics. It is found that the kinetic parameters estimated by multi‐response technique have a superior predictive capability as opposed to those obtained using the terminal CSD only. An advantage of the proposed algorithm is that the measurement of course of supersaturation, which is difficult to perform, is not r