1 #!/usr/bin/python
   2 # Uloha COUL1
   3 # Vytvoril Radek Podgorny <radek@podgorny.cz>
   4 
   5 from math import *
   6 
   7 g = 9.806 #grav. konstanta
   8 E = 8.854187817 * 10**-12 #epsilon
   9 
  10 # Tady zadej svoje parametry
  11 L = 0.996 #delka zavesu
  12 d = 0.447 #vzdalenost od osy
  13 a = 0.00145 #polomer kulicky
  14 R = float(8520) #hustota kulicky
  15 
  16 
  17 V = float(4)/float(3) * pi * a**3 #objem
  18 m = V * R #hmotnost
  19 
  20 Fg = g * m #grav. sila
  21 v = (L**2 - d**2)**0.5 #vyska
  22 Fe = Fg * d/v #bocni sila
  23 
  24 Q = (Fe * float(4)*pi*E * (2*d)**2)**0.5 #naboj
  25 
  26 E1 = Q / (float(4)*pi*E * a**2) #intenz. na povrchu
  27 
  28 E2 = Q / (float(4)*pi*E * L**2) #intenz. v miste zavesu od jedne kulicky
  29 alfa = acos(d/L) #uhel
  30 E2y = float(2) * E2 * sin(alfa)
  31 E2x = 0 #tohle je jasny
  32 
  33 fi = float(2) * Q / (float(4)*pi*E * d) #potencial mezi kulickami (nula v nekon.)
  34 fi -= float(2) * Q / (float(4)*pi*E * L) #hodime nulu do mista zavesu
  35 
  36 
  37 print("Hmotnost je " + str(m) + " kg")
  38 print("Naboj je " + str(Q) + " C")
  39 print("Intenzita E1 je " + str(E1) + " V/m")
  40 print("Intenzita E2 a E2y jsou " + str(E2y) + " a Ex 0 V/m")
  41 print("Potencial mezi kulickami je " + str(fi))