import numpy as np
import scipy

a = np.concatenate(([3], [0]*5, np.arange(-1, 1.002, 2/9.0)))
print("a = ", a)

print("np.r_[3,[0]*5,-1:1:10j] = ", np.r_[3,[0]*5,-1:1:10j])

a =  [ 3.          0.          0.          0.          0.          0.
 -1.         -0.77777778 -0.55555556 -0.33333333 -0.11111111  0.11111111
  0.33333333  0.55555556  0.77777778  1.        ]
np.r_[3,[0]*5,-1:1:10j] =  [ 3.          0.          0.          0.          0.          0.
 -1.         -0.77777778 -0.55555556 -0.33333333 -0.11111111  0.11111111
  0.33333333  0.55555556  0.77777778  1.        ]

print("np.mgrid[0:5:4j,0:5:4j] = ", np.mgrid[0:5:4j,0:5:4j])

np.mgrid[0:5:4j,0:5:4j] =  [[[0.         0.         0.         0.        ]
  [1.66666667 1.66666667 1.66666667 1.66666667]
  [3.33333333 3.33333333 3.33333333 3.33333333]
  [5.         5.         5.         5.        ]]

 [[0.         1.66666667 3.33333333 5.        ]
  [0.         1.66666667 3.33333333 5.        ]
  [0.         1.66666667 3.33333333 5.        ]
  [0.         1.66666667 3.33333333 5.        ]]]

p = np.poly1d([3,4,5])
print("p = ", p)
print("p*p=", p*p)

print('p.integ(k=6) = ', p.integ(k=6))

p =     2
3 x + 4 x + 5
p*p=    4      3      2
9 x + 24 x + 46 x + 40 x + 25
p.integ(k=6) =     3     2
1 x + 2 x + 5 x + 6