code: def simple_2model():
      mu = -2.1
      tau = 1.3
      p = 0.4
      with Model() as model:
          x = pm.Normal("x", mu, tau=tau, initval=0.1)
          pm.Deterministic("logx", pt.log(x))
          pm.Bernoulli("y", p)
      return model.initial_point(), model
  
imports: from itertools import product
  import numpy
  import pytensor
  import pytensor.tensor
  from pytensor.compile.ops import as_op
  import pymc
  from pymc import Categorical, Metropolis, Model, Normal
  from pymc.pytensorf import floatX_array
label_desc: Use image node_id 7 for calling a global function with example usage: simple_2model() and returns: model,  model