see local file "/usr/lib/python2.7/dist-packages/sage/quadratic_forms/quadratic_formmassSiegel_densities.py"

Q = DiagonalQuadraticForm(ZZ, [1,1]);
for p in [3,5,7,11,13,17,19] :
    print([p,Q.Pall_mass_density_at_odd_prime(p)])
[(0, Quadratic form in 2 variables over Integer Ring with coefficients: 
[ 1 0 ]
[ * 1 ])] [(0, 2, 4)] [2/3] 2/3
[3, 2/3]
[(0, Quadratic form in 2 variables over Integer Ring with coefficients: 
[ 1 0 ]
[ * 1 ])] [(0, 2, 4)] [6/5] 6/5
[5, 6/5]
[(0, Quadratic form in 2 variables over Integer Ring with coefficients: 
[ 1 0 ]
[ * 1 ])] [(0, 2, 4)] [6/7] 6/7
[7, 6/7]
[(0, Quadratic form in 2 variables over Integer Ring with coefficients: 
[ 1 0 ]
[ * 1 ])] [(0, 2, 4)] [10/11] 10/11
[11, 10/11]
[(0, Quadratic form in 2 variables over Integer Ring with coefficients: 
[ 1 0 ]
[ * 1 ])] [(0, 2, 4)] [14/13] 14/13
[13, 14/13]
[(0, Quadratic form in 2 variables over Integer Ring with coefficients: 
[ 1 0 ]
[ * 1 ])] [(0, 2, 4)] [18/17] 18/17
[17, 18/17]
[(0, Quadratic form in 2 variables over Integer Ring with coefficients: 
[ 1 0 ]
[ * 1 ])] [(0, 2, 4)] [18/19] 18/19
[19, 18/19]
Q=QuadraticForm(ZZ, 4, [1,1,1,1,1,0,1,2,2,2])
Q.Gram_matrix_rational()
[  1 1/2 1/2 1/2]
[1/2   1   0 1/2]
[1/2   0   2   1]
[1/2 1/2   1   2]
Q = DiagonalQuadraticForm(ZZ, [1,3])
Q.Pall_mass_density_at_odd_prime(3)
[(0, Quadratic form in 1 variables over Integer Ring with coefficients: 
[ 1 ]), (1, Quadratic form in 1 variables over Integer Ring with coefficients: 
[ 1 ])] [(0, 1, 2), (1, 1, 2)] [1, 1] 2
2
Q.Pall_mass_density_at_odd_prime(5)
[(0, Quadratic form in 2 variables over Integer Ring with coefficients: 
[ 1 0 ]
[ * 3 ]), (1, Quadratic form in 2 variables over Integer Ring with coefficients: 
[ 12 0 ]
[ * 36 ])] [(0, 2, 12), (1, 2, 1728)] [4/5, 4/5] 32/125
32/125
Q=QuadraticForm(ZZ, 4, [2,2,2,2,2,0,2,4,4,4])
Q.Pall_mass_density_at_odd_prime(5)
[(0, Quadratic form in 2 variables over Integer Ring with coefficients: 
[ 2 0 ]
[ * 6 ]), (1, Quadratic form in 2 variables over Integer Ring with coefficients: 
[ 24 0 ]
[ * 72 ])] [(0, 2, 48), (1, 2, 6912)] [4/5, 4/5] 32/125
32/125
Q = DiagonalQuadraticForm(ZZ, [1,1,1])
Q.Kitaoka_mass_at_2()
1/2
Q = DiagonalQuadraticForm(ZZ, [1,1,1])
Q.Watson_mass_at_2()
384
Print("a")
---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-15-b33d3f360c16> in <module>()
----> 1 Print("a")

NameError: name 'Print' is not defined