kolpeex

require 'complex.rb'

J = Complex.new(0,5e-3)
E1 = Complex.new(0,50)
E2 = 40
R1 = 5e+3
R2 = 6e+3
L1 = 318*1e-3
L2 = 318*1e-3
C1 = 3.18*1e-9
C2 = 6.36*1e-9
f = 5e+3
w = f

Zl1 = Complex.new(0, w*L1)
Zl2 = Complex.new(0, w*L2)
Zc1 = Complex.new(0, -1/(w*C1))
Zc2 = Complex.new(0, -1/(w*C2))

G12 = 1/Zc1
G23 = (1/(R1+Zl1)) + (1/Zc2)
G13_ = 1/(Zl2+R2)
G13 = 0 + G13_

V1 = 0
V3 = (E2*G13-J-G23*E1*G12/(G12+G23)) / (G13+G23-G23*G23/(G12+G23))
V2 = (V3*G23-E1*G12)/(G12+G23)

Zo1 = Zc1+Zl2+R2+Zc2*(R1+Zl1)/(Zc2+R1+Zl1)
Uo2 = J*(Zl2+R2)*(Zc1+Zc2*(R1+Zl1)/(Zc2+R1+Zl1))/(Zl2+R2+Zc1+Zc2*(R1+Zl1)/(Zc2+R1+Zl1))
Zo3 = Zl2+R2+Zc1+(R1+Zl1)*Zc2/(R1+Zl1+Zc2)

i1 = -E1/Zo1 + Uo2/(Zc1+(Zc2*(R1+Zl1))/(Zc2+R1+Zl1)) - E2/Zo3
i2 = E1/Zo1 + Uo2/(Zl2+R2) + E2/Zo3

i3 = -(E1-(Zc1+Zl2+R2)*E1/Zo1)/(R1+Zl1) +
  (Uo2-Zc1*Uo2/(Zc1+(Zc2*(R1+Zl1))/(Zc2+R1+Zl1)))/(R1+Zl1) -
  (E2-(Zc1+Zl2+R2)*E2/Zo1)/(R1+Zl1)
i4 = ((V2-V3)/Zc2)

potr= i3*i3*R1 +
      i2*i2*R2 +
      i3*i3*Zl1 +
      i2*i2*Zl2 +
      i1*i1*Zc1 +
      i4*i4*Zc2

print   i3*i3*R1
print "\n"
print        i2*i2*R2
print "\n"
print        i3*i3*Zl1
print "\n"
print        i2*i2*Zl2
print "\n"
print        i1*i1*Zc1
print "\n"
print        i4*i4*Zc2
print "\n"
print "\n"

isp = -E1*i1 + E2*i2 + J*(-V3)

print potr
print "\n"
print isp
print "\n"