Python Model.eliminate示例

示例#1

文件： diffeq.py 项目： dkasak/pacal

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
figure()
M.plot()
print M.nddistr.pdf(linspace(0,2.5,100))
X0 = BetaDistr(1, 1)
y = X0 * exp(A)
y.summary()
y.plot(label="Y0*exp(A)")
figure()
err = y.get_piecewise_pdf() - M.as1DDistr().get_piecewise_pdf()

示例#2

文件： diffeq.py 项目： ntankovic/pacal

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
figure()
M.plot()
print M.nddistr.pdf(linspace(0, 2.5, 100))
X0 = BetaDistr(1, 1)
y = X0 * exp(A)
y.summary()
y.plot(label="Y0*exp(A)")
figure()
err = y.get_piecewise_pdf() - M.as1DDistr().get_piecewise_pdf()