LSODA
can also automatically detect stiff systems of ODEs and adjust its integration scheme so as to handle them
efficiently.