#ODE for BS22
## note that Der[z] = (d/dz) operator

9216*z^3*(-26625 + 1720480*z^2 + 51739038*z^4 - 859601508*z^6 - 
   11749680884*z^8 + 45406216016*z^10 + 491812757184*z^12 - 
   7942039456976*z^14 - 13310875733504*z^16 - 117988763318592*z^18 - 
   111055585700864*z^20 + 363673147745280*z^22 + 46309759475712*z^24 + 
   1049345483702272*z^26 - 805406283661312*z^28 - 2848015276572672*z^30 - 
   618200865701888*z^32 - 47188470136832*z^34 + 23342610382848*z^36 + 
   277025390592*z^38) + 36*(-1000 - 125385*z^2 - 71990140*z^4 - 
   120604500*z^6 + 78114649052*z^8 + 796182580508*z^10 - 
   60846368251776*z^12 + 186066871948464*z^14 + 2928324026033856*z^16 - 
   17266157745076224*z^18 - 18585855123550976*z^20 - 
   582535702683573248*z^22 - 413925863971700736*z^24 + 
   520066966066135040*z^26 - 1894078532639260672*z^28 + 
   4416192752646094848*z^30 + 5685559397242109952*z^32 - 
   10352230728883765248*z^34 - 2741925169218977792*z^36 - 
   197523937477787648*z^38 + 148477500459057152*z^40 + 2056636499755008*z^42)*
  Der[z] + 36*z*(1000 + 147385*z^2 - 161071620*z^4 - 4776616524*z^6 + 
   56880629988*z^8 + 5844799651108*z^10 - 115041489722240*z^12 + 
   58558457928016*z^14 + 7288935096052800*z^16 - 23726291558683136*z^18 + 
   92302370144040704*z^20 - 1234161083190774784*z^22 - 
   665746032772186112*z^24 + 1274368110092058624*z^26 - 
   8030313621102264320*z^28 + 2666456123139162112*z^30 + 
   31450776193434910720*z^32 - 15182972897200177152*z^34 - 
   5643877030652018688*z^36 - 660757471979110400*z^38 + 
   394557892201545728*z^40 + 6181729249263616*z^42)*Der[z]^2 + 
 3*(200 + 145077*z^2 + 12838748*z^4 - 342578876*z^6 - 38395064748*z^8 - 
   539777568268*z^10 + 67503472139968*z^12 - 706943411694192*z^14 - 
   2909655044410816*z^16 + 58194283431431680*z^18 - 183445221587972352*z^20 + 
   1704026816198503424*z^22 - 8761277686740795392*z^24 - 
   4496063395371909120*z^26 + 18072860021693087744*z^28 - 
   73367104702260445184*z^30 - 71449186120765538304*z^32 + 
   374533068794347126784*z^34 - 53750265260850282496*z^36 - 
   39680713257979478016*z^38 - 9169141202795102208*z^40 + 
   3538955195369127936*z^42 + 62266442992582656*z^44)*Der[z]^3 + 
 3*z*(-200 + 108923*z^2 + 10645472*z^4 + 94094764*z^6 - 11520805796*z^8 - 
   473247703172*z^10 + 20667672522288*z^12 - 102673585289360*z^14 - 
   1565960668522112*z^16 + 14215057789896192*z^18 - 70151659507958528*z^20 + 
   627334116421614592*z^22 - 1877204332784566272*z^24 - 
   1442100417067384832*z^26 + 7212603053436043264*z^28 - 
   14573898330097582080*z^30 - 49431946553785843712*z^32 + 
   119505298694455951360*z^34 + 3629758095250423808*z^36 - 
   8032065939224657920*z^38 - 3898380272269262848*z^40 + 
   978218177924169728*z^42 + 19360750497693696*z^44)*Der[z]^4 + 
 12*z^2*(-1 + 2*z)*(1 + 2*z)*(100 + 2541*z^2 - 402368*z^4 - 16320304*z^6 + 
   403999410*z^8 + 26872452988*z^10 - 469901061032*z^12 - 
   1984678852064*z^14 + 47110956136928*z^16 - 189867612164096*z^18 + 
   2116125071942144*z^20 - 10724592327776256*z^22 - 6361659161698304*z^24 + 
   33726897275142144*z^26 - 100150342418432000*z^28 - 
   329847444032454656*z^30 + 894763323596210176*z^32 + 
   119871444699054080*z^34 - 39727066119143424*z^36 - 
   37106498802810880*z^38 + 6624162420359168*z^40 + 147746874982400*z^42)*
  Der[z]^5 + z^3*(-1 + 2*z)^2*(1 + 2*z)^2*(-1 + 4*z)*(1 + 4*z)*
  (1 + 12*z^2 - 384*z^4 + 64*z^6)*(200 + 17223*z^2 + 359016*z^4 - 
   25768772*z^6 - 160112268*z^8 + 2074374156*z^10 + 3094012848*z^12 + 
   161340317424*z^14 + 95008052544*z^16 - 154402840576*z^18 + 
   2150587987968*z^20 + 4014853226496*z^22 - 16673457176576*z^24 - 
   6351589736448*z^26 - 559700508672*z^28 + 786180341760*z^30 + 
   17314086912*z^32)*Der[z]^6