weight = 5 level = 38 mult = 3 thetaBlock = {1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 4} note = "*** singularPart and weaklyHolomorphicJacobiForm are expanded to \ q^(level/4). ***" singularPart = 10 + q^(-1) + z^(-4) + 3/z^3 + 6/z^2 + 9/z + 9*z + 6*z^2 + 3*z^3 + z^4 + q*(z^(-13) + z^13) + q^4*(z^(-25) + z^25) + q^8*(3/z^35 + 3*z^35) + q^9*(9/z^37 + 9*z^37) note = "*** Humbert divisors are in the form {{D, r}, multiplicity}. ***" HumbertDRmultiplicities = {{{1, 1}, 19}, {{1, 37}, 19}, {{4, 2}, 7}, {{9, 3}, 3}, {{9, 35}, 3}, {{16, 4}, 1}, {{17, 13}, 1}, {{17, 25}, 1}, {{152, 0}, 1}} weaklyHolomorphicJacobiForm = 10 + q^(-1) + z^(-4) + 3/z^3 + 6/z^2 + 9/z + 9*z + 6*z^2 + 3*z^3 + z^4 + q*(21456 + z^(-13) + 21/z^12 + 96/z^11 + 326/z^10 + 823/z^9 + 1795/z^8 + 3363/z^7 + 5684/z^6 + 8635/z^5 + 12098/z^4 + 15558/z^3 + 18638/z^2 + 20676/z + 20676*z + 18638*z^2 + 15558*z^3 + 12098*z^4 + 8635*z^5 + 5684*z^6 + 3363*z^7 + 1795*z^8 + 823*z^9 + 326*z^10 + 96*z^11 + 21*z^12 + z^13) + q^2*(1999528 + 31/z^17 + 249/z^16 + 1196/z^15 + 4092/z^14 + 11624/z^13 + 28213/z^12 + 60806/z^11 + 118178/z^10 + 210605/z^9 + 346587/z^8 + 531596/z^7 + 763054/z^6 + 1030244/z^5 + 1311651/z^4 + 1579544/z^3 + 1801140/z^2 + 1948306/z + 1948306*z + 1801140*z^2 + 1579544*z^3 + 1311651*z^4 + 1030244*z^5 + 763054*z^6 + 531596*z^7 + 346587*z^8 + 210605*z^9 + 118178*z^10 + 60806*z^11 + 28213*z^12 + 11624*z^13 + 4092*z^14 + 1196*z^15 + 249*z^16 + 31*z^17) + q^3*(73150192 + 31/z^21 + 399/z^20 + 2408/z^19 + 10434/z^18 + 35841/z^17 + 104506/z^16 + 266809/z^15 + 611636/z^14 + 1277327/z^13 + 2460198/z^12 + 4405518/z^11 + 7386510/z^10 + 11652657/z^9 + 17371733/z^8 + 24549852/z^7 + 32980242/z^6 + 42201200/z^5 + 51526981/z^4 + 60104528/z^3 + 67050058/z^2 + 71576021/z + 71576021*z + 67050058*z^2 + 60104528*z^3 + 51526981*z^4 + 42201200*z^5 + 32980242*z^6 + 24549852*z^7 + 17371733*z^8 + 11652657*z^9 + 7386510*z^10 + 4405518*z^11 + 2460198*z^12 + 1277327*z^13 + 611636*z^14 + 266809*z^15 + 104506*z^16 + 35841*z^17 + 10434*z^18 + 2408*z^19 + 399*z^20 + 31*z^21) + q^4*(1601011450 + z^(-25) + 99/z^24 + 1196/z^23 + 7686/z^22 + 35841/z^21 + 133714/z^20 + 423862/z^19 + 1178580/z^18 + 2943040/z^17 + 6700439/z^16 + 14075227/z^15 + 27514096/z^14 + 50400159/z^13 + 86973326/z^12 + 142018460/z^11 + 220203160/z^10 + 325167524/z^9 + 458377241/z^8 + 618080481/z^7 + 798492964/z^6 + 989667842/z^5 + 1178035344/z^4 + 1347880090/z^3 + 1483329402/z^2 + 1570782629/z + 1570782629*z + 1483329402*z^2 + 1347880090*z^3 + 1178035344*z^4 + 989667842*z^5 + 798492964*z^6 + 618080481*z^7 + 458377241*z^8 + 325167524*z^9 + 220203160*z^10 + 142018460*z^11 + 86973326*z^12 + 50400159*z^13 + 27514096*z^14 + 14075227*z^15 + 6700439*z^16 + 2943040*z^17 + 1178580*z^18 + 423862*z^19 + 133714*z^20 + 35841*z^21 + 7686*z^22 + 1196*z^23 + 99*z^24 + z^25) + q^5*(24987579358 + 96/z^27 + 1504/z^26 + 11624/z^25 + 62886/z^24 + 266809/z^23 + 949990/z^22 + 2943040/z^21 + 8137839/z^20 + 20428920/z^19 + 47163076/z^18 + 101088229/z^17 + 202688687/z^16 + 382428044/z^15 + 682284922/z^14 + 1155522166/z^13 + 1863854084/z^12 + 2871057991/z^11 + 4233160262/z^10 + 5985634966/z^9 + 8129874486/z^8 + 10621183260/z^7 + 13362033244/z^6 + 16202736662/z^5 + 18951945359/z^4 + 21395865693/z^3 + 23324760138/z^2 + 24561446644/z + 24561446644*z + 23324760138*z^2 + 21395865693*z^3 + 18951945359*z^4 + 16202736662*z^5 + 13362033244*z^6 + 10621183260*z^7 + 8129874486*z^8 + 5985634966*z^9 + 4233160262*z^10 + 2871057991*z^11 + 1863854084*z^12 + 1155522166*z^13 + 682284922*z^14 + 382428044*z^15 + 202688687*z^16 + 101088229*z^17 + 47163076*z^18 + 20428920*z^19 + 8137839*z^20 + 2943040*z^21 + 949990*z^22 + 266809*z^23 + 62886*z^24 + 11624*z^25 + 1504*z^26 + 96*z^27) + q^6*(305272054792 + 22/z^30 + 823/z^29 + 8940/z^28 + 60806/z^27 + 308660/z^26 + 1277327/z^25 + 4516358/z^24 + 14075227/z^23 + 39466900/z^22 + 101088229/z^21 + 239165032/z^20 + 527240640/z^19 + 1090416424/z^18 + 2127410718/z^17 + 3933118244/z^16 + 6916314585/z^15 + 11604250970/z^14 + 18625567651/z^13 + 28662866870/z^12 + 42371584462/z^11 + 60266897194/z^10 + 82591526473/z^9 + 109183739762/z^8 + 139374612644/z^7 + 171941397230/z^6 + 205142648546/z^5 + 236843393110/z^4 + 264726335623/z^3 + 286560098312/z^2 + 300486234870/z + 300486234870*z + 286560098312*z^2 + 264726335623*z^3 + 236843393110*z^4 + 205142648546*z^5 + 171941397230*z^6 + 139374612644*z^7 + 109183739762*z^8 + 82591526473*z^9 + 60266897194*z^10 + 42371584462*z^11 + 28662866870*z^12 + 18625567651*z^13 + 11604250970*z^14 + 6916314585*z^15 + 3933118244*z^16 + 2127410718*z^17 + 1090416424*z^18 + 527240640*z^19 + 239165032*z^20 + 101088229*z^21 + 39466900*z^22 + 14075227*z^23 + 4516358*z^24 + 1277327*z^25 + 308660*z^26 + 60806*z^27 + 8940*z^28 + 823*z^29 + 22*z^30) + q^7*(3089721626669 + 168/z^32 + 3363/z^31 + 32426/z^30 + 210605/z^29 + 1058677/z^28 + 4405518/z^27 + 15825312/z^26 + 50400159/z^25 + 145045557/z^24 + 382428044/z^23 + 933574978/z^22 + 2127410718/z^21 + 4555139860/z^20 + 9212992068/z^19 + 17678686406/z^18 + 32302115317/z^17 + 56373953280/z^16 + 94215785604/z^15 + 151124448222/z^14 + 233099806460/z^13 + 346307407522/z^12 + 496263667518/z^11 + 686797800662/z^10 + 918908678156/z^9 + 1189708145696/z^8 + 1491671669497/z^7 + 1812424345192/z^6 + 2135219435344/z^5 + 2440172930241/z^4 + 2706159734889/z^3 + 2913149780818/z^2 + 3044629304356/z + 3044629304356*z + 2913149780818*z^2 + 2706159734889*z^3 + 2440172930241*z^4 + 2135219435344*z^5 + 1812424345192*z^6 + 1491671669497*z^7 + 1189708145696*z^8 + 918908678156*z^9 + 686797800662*z^10 + 496263667518*z^11 + 346307407522*z^12 + 233099806460*z^13 + 151124448222*z^14 + 94215785604*z^15 + 56373953280*z^16 + 32302115317*z^17 + 17678686406*z^18 + 9212992068*z^19 + 4555139860*z^20 + 2127410718*z^21 + 933574978*z^22 + 382428044*z^23 + 145045557*z^24 + 50400159*z^25 + 15825312*z^26 + 4405518*z^27 + 1058677*z^28 + 210605*z^29 + 32426*z^30 + 3363*z^31 + 168*z^32) + q^8*(26899893142184 + 3/z^35 + 472/z^34 + 8635/z^33 + 81098/z^32 + 531596/z^31 + 2725644/z^30 + 11652657/z^29 + 43150679/z^28 + 142018460/z^27 + 422989368/z^26 + 1155522166/z^25 + 2925013368/z^24 + 6916314585/z^23 + 15374960938/z^22 + 32302115317/z^21 + 64418904255/z^20 + 122389567120/z^19 + 222209066158/z^18 + 386551528932/z^17 + 645745030012/z^16 + 1037933272134/z^15 + 1607926632932/z^14 + 2404308131454/z^13 + 3474541878743/z^12 + 4858185765312/z^11 + 6578744521040/z^10 + 8635204085180/z^9 + 10994663980334/z^8 + 13587702299833/z^7 + 16307948192888/z^6 + 19016887149646/z^5 + 21554075550675/z^4 + 23751992839965/z^3 + 25453731070816/z^2 + 26531040204493/z + 26531040204493*z + 25453731070816*z^2 + 23751992839965*z^3 + 21554075550675*z^4 + 19016887149646*z^5 + 16307948192888*z^6 + 13587702299833*z^7 + 10994663980334*z^8 + 8635204085180*z^9 + 6578744521040*z^10 + 4858185765312*z^11 + 3474541878743*z^12 + 2404308131454*z^13 + 1607926632932*z^14 + 1037933272134*z^15 + 645745030012*z^16 + 386551528932*z^17 + 222209066158*z^18 + 122389567120*z^19 + 64418904255*z^20 + 32302115317*z^21 + 15374960938*z^22 + 6916314585*z^23 + 2925013368*z^24 + 1155522166*z^25 + 422989368*z^26 + 142018460*z^27 + 43150679*z^28 + 11652657*z^29 + 2725644*z^30 + 531596*z^31 + 81098*z^32 + 8635*z^33 + 472*z^34 + 3*z^35) + q^9*(206859933485248 + 9/z^37 + 872/z^36 + 15558/z^35 + 151288/z^34 + 1030244/z^33 + 5507447/z^32 + 24549852/z^31 + 94792200/z^30 + 325167524/z^29 + 1009079594/z^28 + 2871057991/z^27 + 7567146088/z^26 + 18625567651/z^25 + 43092660637/z^24 + 94215785604/z^23 + 195518901366/z^22 + 386551528932/z^21 + 730373043486/z^20 + 1322414540016/z^19 + 2299754663558/z^18 + 3849083373988/z^17 + 6210886191327/z^16 + 9676829388782/z^15 + 14577296867542/z^14 + 21256687396900/z^13 + 30035680013311/z^12 + 41162035446832/z^11 + 54754485301238/z^10 + 70747119832661/z^9 + 88843928464612/z^8 + 108493764357039/z^7 + 128894764818978/z^6 + 149033648127580/z^5 + 167760088828516/z^4 + 183889978433280/z^3 + 196325547890966/z^2 + 204175878226117/z + 204175878226117*z + 196325547890966*z^2 + 183889978433280*z^3 + 167760088828516*z^4 + 149033648127580*z^5 + 128894764818978*z^6 + 108493764357039*z^7 + 88843928464612*z^8 + 70747119832661*z^9 + 54754485301238*z^10 + 41162035446832*z^11 + 30035680013311*z^12 + 21256687396900*z^13 + 14577296867542*z^14 + 9676829388782*z^15 + 6210886191327*z^16 + 3849083373988*z^17 + 2299754663558*z^18 + 1322414540016*z^19 + 730373043486*z^20 + 386551528932*z^21 + 195518901366*z^22 + 94215785604*z^23 + 43092660637*z^24 + 18625567651*z^25 + 7567146088*z^26 + 2871057991*z^27 + 1009079594*z^28 + 325167524*z^29 + 94792200*z^30 + 24549852*z^31 + 5507447*z^32 + 1030244*z^33 + 151288*z^34 + 15558*z^35 + 872*z^36 + 9*z^37)