Numerical solution of ordinary simultaneous differential equations of the 1st order using a method for automatic step change