Plots of zeroj0.y and zeroj1.y should display steps with
lower-right corners at ( t , y ) = ( xi,Modelica ,
xi,GSL ) and such values should be equal at these
corners.
As not all special functions' zeros are provided via GSL
routines, the displayed use of Modelica's zero-crossing logical
blocks could be useful for finding such information for these other
functions.