weight = 8
level = 17
mult = 15
thetaBlock = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2}
note = "*** singularPart and weaklyHolomorphicJacobiForm are expanded to \
q^(level/4). ***"
singularPart = 16 + q^(-1) + 6/z^2 + 10/z + 10*z + 6*z^2 + 
     q^2*(z^(-12) + z^12)
note = "*** Humbert divisors are  in the form {{D, r}, multiplicity}. ***"
HumbertDRmultiplicities = {{{1, 1}, 16}, {{4, 2}, 6}, {{8, 12}, 1}, 
     {{17, 17}, 1}, {{68, 0}, 1}}
weaklyHolomorphicJacobiForm = 16 + q^(-1) + 6/z^2 + 10/z + 10*z + 6*z^2 + 
     q*(32018 + 30/z^8 + 290/z^7 + 1265/z^6 + 3800/z^5 + 8541/z^4 + 
       15540/z^3 + 23317/z^2 + 29650/z + 29650*z + 23317*z^2 + 15540*z^3 + 
       8541*z^4 + 3800*z^5 + 1265*z^6 + 290*z^7 + 30*z^8) + 
     q^2*(2989838 + z^(-12) + 158/z^11 + 1901/z^10 + 11028/z^9 + 43611/z^8 + 
       131132/z^7 + 319886/z^6 + 653982/z^5 + 1148437/z^4 + 1756542/z^3 + 
       2365197/z^2 + 2820086/z + 2820086*z + 2365197*z^2 + 1756542*z^3 + 
       1148437*z^4 + 653982*z^5 + 319886*z^6 + 131132*z^7 + 43611*z^8 + 
       11028*z^9 + 1901*z^10 + 158*z^11 + z^12) + 
     q^3*(109364106 + 60/z^14 + 1738/z^13 + 16821/z^12 + 98756/z^11 + 
       416467/z^10 + 1379996/z^9 + 3767927/z^8 + 8750350/z^7 + 17641897/z^6 + 
       31341534/z^5 + 49569280/z^4 + 70329228/z^3 + 89977656/z^2 + 
       104176222/z + 104176222*z + 89977656*z^2 + 70329228*z^3 + 
       49569280*z^4 + 31341534*z^5 + 17641897*z^6 + 8750350*z^7 + 
       3767927*z^8 + 1379996*z^9 + 416467*z^10 + 98756*z^11 + 16821*z^12 + 
       1738*z^13 + 60*z^14) + q^4*(2393664106 + 195/z^16 + 5442/z^15 + 
       59041/z^14 + 389896/z^13 + 1865223/z^12 + 7014342/z^11 + 
       21782998/z^10 + 57625424/z^9 + 132698880/z^8 + 270087732/z^7 + 
       491418136/z^6 + 806019886/z^5 + 1199341857/z^4 + 1626563098/z^3 + 
       2017477969/z^2 + 2293745956/z + 2293745956*z + 2017477969*z^2 + 
       1626563098*z^3 + 1199341857*z^4 + 806019886*z^5 + 491418136*z^6 + 
       270087732*z^7 + 132698880*z^8 + 57625424*z^9 + 21782998*z^10 + 
       7014342*z^11 + 1865223*z^12 + 389896*z^13 + 59041*z^14 + 5442*z^15 + 
       195*z^16)