#------------------------------------------------------
# potenza in sinusoidale: circuito puramente capacitivo
#
#  GdA dicembre 2021
#------------------------------------------------------


f0 = 12   # V
R  = 50   # Ohm 
nu =1000  # Hz
om = 2*pi*nu
T = 1/nu

curve(cos(om*x*T), -0.5, 1.5, xlab='t / T',
      ylab='f(t)/f0;  P(t)/Pmax', lwd=3, col='blue')
abline(h=0)
grid()
curve(cos(om*x*T+pi/2), lwd=3, col='green', add=TRUE)
curve(-sin(2*om*x*T), lwd=3, col='red', n=201, add=TRUE)
abline(h=0, lty=2, col=2)