#!/usr/bin/env python
# -*- coding: utf8 -*-
# Soubor:  20111004-Tue-stridavy-vykon.py
# Datum:   04.10.2011 08:21
# Autor:   Marek Nožka, nozka <@t> spseol <d.t> cz
# Licence: GNU/GPL 

import pylab as pl
from pylab import pi, e, sin

print pi
print e    

f=77      # frekvence
fi = 50   # fázový posun
fi = input("zadej fázový posun: ")
##############################################
# vytvořím si časovou osu
t = []
krok = 0.2
n = 0 # číslo vzorku
vzorek = 0
while vzorek<=20:
    t.append(vzorek)
    n = n+1
    vzorek = krok * n
###############################################
# vytvořím vzorky namětí
u = []
n = 0
while n < len(t):
    u.append( sin( 2*pi*f*t[n]/1000) )
    n = n + 1

###############################################
# vytvořím vzorky proudu
i = []
n = 0
while n < len(t):
    i.append( sin( 2*pi*f*t[n]/1000 + fi*pi/180) )
    n = n + 1
    
###############################################
# vypočet výkonu
p = []
n = 0
while n < len(t):
    p.append( u[n] * i[n] )
    n = n+1

pl.plot(t,u)
pl.plot(t,i)
pl.plot(t,p)
pl.grid()
pl.show()