X0 = BetaDistr(2, 2) y = X0 * exp(A) y.summary() y.plot(label="Y0*exp(A)") Y[-1].plot('r') figure() err = y.get_piecewise_pdf() - M3.as1DDistr().get_piecewise_pdf() err.plot() show() stop print "---", [K] + Y print M M.varschange(A, K) print M for i in xrange(n): M.varschange(Y[i], Y[i+1]) print M M.varschange(K, A) M.plot() #M.condition(Y[n], 2) print M M.eliminate(K) print M for i in xrange(n-1,-1,-1): M.eliminate(Y[i]) print M M.eliminate(A) print M
X0 = BetaDistr(2, 2) y = X0 * exp(A) y.summary() y.plot(label="Y0*exp(A)") Y[-1].plot('r') figure() err = y.get_piecewise_pdf() - M3.as1DDistr().get_piecewise_pdf() err.plot() show() stop print "---", [K] + Y print M M.varschange(A, K) print M for i in xrange(n): M.varschange(Y[i], Y[i + 1]) print M M.varschange(K, A) M.plot() #M.condition(Y[n], 2) print M M.eliminate(K) print M for i in xrange(n - 1, -1, -1): M.eliminate(Y[i]) print M M.eliminate(A) print M