code: def _(dist: GammaDistribution, size, seed):
      # https://stackoverflow.com/questions/42150965/how-to-plot-gamma-distribution-with-alpha-and-beta-parameters-in-python
      return scipy.stats.gamma.rvs(
          a=float(dist.k),
          scale=float(dist.theta),
          size=size,
          random_state=seed,
      )
  
imports: from functools import singledispatch
  from sympy.core.symbol import Dummy
  from sympy.functions.elementary.exponential import exp
  from sympy.utilities.lambdify import lambdify
  from sympy.external import import_module
  from sympy.stats import DiscreteDistributionHandmade
  from sympy.stats.crv import SingleContinuousDistribution
  from sympy.stats.crv_types import ChiSquaredDistribution, ExponentialDistribution, GammaDistribution, LogNormalDistribution, NormalDistribution, ParetoDistribution, UniformDistribution, BetaDistribution, StudentTDistribution, CauchyDistribution
  from sympy.stats.drv_types import GeometricDistribution, LogarithmicDistribution, NegativeBinomialDistribution, PoissonDistribution, SkellamDistribution, YuleSimonDistribution, ZetaDistribution
  from sympy.stats.frv import SingleFiniteDistribution
label_desc: Use image node_id 5 for calling a global function with example usage: _(dist, size, seed) and returns: scipy