Example #1
0
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
Example #2
0
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