function JH_data(N)
   t0 := Cputime();
   ans := [* *];
   for d in Divisors((N-1) div 2) do 
      if d eq (N-1) div 2 then
         continue; 
      end if;
      H := (N-1) div d;
      t := Cputime();
      M := ModularSymbolsH(N, H,2,0);
      A := CuspidalSubspace(M);
      dim := Dimension(A) div 2;
      L := LRatio(A,1);
      T := TorsionBound(A,60);
      Hbar := H div 2;   // order of quotient of H by {+-1}.
      Phi := Hbar/GCD(Hbar,6);
      Append(~ans,[* N, d, dim, L, T, T^2*L, Phi, Cputime(t) *]);
   end for;
   return ans;
end function;

ans := JH_data(arg);
ans;
