Ejemplo n.º 1
0
#end

dist = pc.Gamma(2)
print pc.orth_bert(2, dist)
#  [1.0, q0-2.0, q0^2-6.0q0+6.0]
#end

dist = pc.Uniform(-1, 1)
print dist.ttr([0, 1, 2, 3])
#  [[ 0.          0.          0.          0.        ]
#   [-0.          0.33333333  0.26666667  0.25714286]]
#end

dist = pc.Lognormal(0.01)
orths = pc.orth_ttr(2, dist)
print orths
#  [1.0, q0-1.00501252086, q0^2-2.04042818514q0+0.842739860094]
#end

dist = pc.Iid(pc.Gamma(1), 2)
orths = pc.orth_ttr(2, dist)
print orths
#  [1.0, q1-1.0, q0-1.0, q1^2-4.0q1+2.0, q0q1-q1-q0+1.0, q0^2-4.0q0+2.0]
#end

q = pc.variable()
dist = pc.Normal()
print pc.E(q, dist)
#  0.0
#end
Ejemplo n.º 2
0
dist = pc.Gamma(2)
print pc.orth_bert(2, dist)
#  [1.0, q0-2.0, q0^2-6.0q0+6.0]
#end

dist = pc.Uniform(-1,1)
print dist.ttr([0,1,2,3])
#  [[ 0.          0.          0.          0.        ]
#   [-0.          0.33333333  0.26666667  0.25714286]]
#end

dist = pc.Lognormal(0.01)
orths = pc.orth_ttr(2, dist)
print orths
#  [1.0, q0-1.00501252086, q0^2-2.04042818514q0+0.842739860094]
#end

dist = pc.Iid(pc.Gamma(1), 2)
orths = pc.orth_ttr(2, dist)
print orths
#  [1.0, q1-1.0, q0-1.0, q1^2-4.0q1+2.0, q0q1-q1-q0+1.0, q0^2-4.0q0+2.0]
#end

q = pc.variable()
dist = pc.Normal()
print pc.E(q, dist)
#  0.0
#end


Ejemplo n.º 3
0
from pylab import *
import current as pc

dist = pc.Iid(pc.Gamma(1), 3)
x, y, z = pc.variable(3)
poly = x + y + z + x * y * z
print poly
#  q0q1q2+q2+q1+q0
#end

print pc.E_cond(poly, z, dist)
#  2.0q2+2.0
#end

print pc.E_cond(poly, [x, y], dist)
#  q0q1+q1+q0+1.0
#end
Ejemplo n.º 4
0
from numpy import *
import current as pc

x, y = pc.variable(2)
print x
# q0
#end

polys = pc.Poly([1, x, x * y])
print polys
#  [1, q0, q0q1]
#end

print pc.basis(4)
#  [1, q0, q0^2, q0^3, q0^4]
#end

print pc.basis(1, 2, dim=2)
#  [q0, q1, q0^2, q0q1, q1^2]
#end

print pc.basis(1, [1, 2])
#  [q0, q1, q0q1, q1^2, q0q1^2]
#end

print pc.basis(1, 2, dim=2, sort="GRI")
#  [q0^2, q0q1, q1^2, q0, q1]
#end

poly = pc.Poly([1, x**2, x * y])
print poly(2, 3)
Ejemplo n.º 5
0
from numpy import *
import current as pc

x,y = pc.variable(2)
print x
# q0
#end

polys = pc.Poly([1, x, x*y])
print polys
#  [1, q0, q0q1]
#end

print pc.basis(4)
#  [1, q0, q0^2, q0^3, q0^4]
#end

print pc.basis(1, 2, dim=2)
#  [q0, q1, q0^2, q0q1, q1^2]
#end

print pc.basis(1, [1, 2])
#  [q0, q1, q0q1, q1^2, q0q1^2]
#end

print pc.basis(1, 2, dim=2, sort="GRI")
#  [q0^2, q0q1, q1^2, q0, q1]
#end

poly = pc.Poly([1, x**2, x*y])
print poly(2, 3)