#--------------------------------------------------
# potenza in sinusoidale: circuito puramente ohmico
#
#  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((1+cos(2*om*x*T))/2, lwd=3, col='red', add=TRUE)
abline(h=1/2, lty=2, col=2)