E1 := 6 * hypergeom([1/3,2/3],[2],27*t*(2-3*t)/(1-4*t)^3)/(1-4*t)/(1-64*t); g2 := (1-4*t) * (1-60*t+120*t^2-64*t^3); J := g2^3 / (1728*(1-t)^2*t^3*(2-3*t)^3*(1-64*t)); H := hypergeom([1/12, 5/12], [1], 1/J) / g2^(1/4); E2 := (t-1) / 2 / (1+6*t) * (4*H + (4*t-1) * diff(H, t)); E := E1 - E2; series(E, t, 20); gfun:-holexprtodiffeq(E, y(t));