code: def log_(t, alpha=1, beta=1):
      """beta*np.log(alpha*(t-beta)+1)"""
      return beta * np.log(alpha * np.maximum(0, t - 10 * beta) + 1)
  
imports: import numpy
  import pandas
  from scipy import stats
  from scipy.optimize import newton
  from scipy.integrate import cumtrapz
label_desc: Use image node_id 4 for calling a global function with example usage: log_(t, alpha, beta) and returns: unknown