Example #1
0
sys.h(Matrix([x[0]]))

# HG!!!
M, J, J_inv = sys.hg()

print "Mapa="
print(M)
print "\n\n"

print "Jacobiano="
print(J)
print "\n\n"

print "Inversa del Jacobiano="
print(J_inv)
print "\n\n"

print "Para modelica:"

for i in xrange(J_inv.lines):
    for j in xrange(J_inv.cols):
        print "Jac_inv[%d,%d] = "%(i+1,j+1),J_inv[i,j],";"
    


I   = 0.023
pprint(calcular_H([-0.2+0.4*I,-0.2-0.4*I,-1]))
pprint(calcular_H([-2+4*I,-2-4*I,-10]))
pprint(calcular_H([-20+40*I,-20-40*I,-100]))
pprint(calcular_H([-200+400*I,-200-400*I,-1000]))