code: def polyvander(x, polybase, order=5):
      polyarr = np.column_stack([polybase(i)(x) for i in range(order)])
      return polyarr
  
imports: from scipy import stats, integrate, special
  import numpy
label_desc: Use image node_id 1 for calling a global function with example usage: polyvander(x, polybase, order) and returns: polyarr