To obtain the dσ/dε you can indeed use commercial software as Matlab and Origin. However, if you are aiming to determine the strain hardening rate in metal alloys, you need to extract the derivative of the stress-strain curve in the upturn which is after the yield point and before the necking. Also consider that it is more accurate to work in terms of true stress - true strain instead the traditional engineering stress - strain.