ps = {bz -> 0, b3 -> 311/600, e3 -> 1/9, e4 -> 2/7, e5 -> 1/8}
 
sysfr1 = {0.07473730671270793866375082859642714519`25., 
     0.12132667010819138083437791781614754242`25., 
     0.14831084933379935525301466082679992356`25., 
     0.21910953776674254588761519618300949815`25., 
     0.36960180159744308337387816515628191061`25., 
     0.39989178365943975658115403249073630133`25., 
     0.05933307677072320593653577363255553156`25., 
     0.36960180159744308337387816513433703604`25., 
     0.37396216756964492111029420792414061357`25., 
     0.33348172166627248449634897435419197919`25., 
     0.25366684792814554466838797143730385817`25., 
     0.37894821875228624425712403331353420247`25., 
     0.31385717811607559535051925750028179046`25., 
     0.33348172166627248449634897435406433252`25., 
     0.31875788856435619357499275590348247923`25., 
     0.18194193368393766381176448610404912399`25., 
     0.18935328168577393859847892838540856445`25., 
     0.28022669075701447864017032340550506526`25., 
     0.18935328168577393859847892838028370723`25., 
     0.28257638511627096502682163121158853094`25., 
     0.19225407329693063866234713092985781288`25.}