Numerical modelling of steady state continuous crystallization processes using piecewise cubic spline functions