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.