#--------------------------------
# tacchino esponenziale
#   

m0 =  1    # kg
v0 = 0.01  # kg/g  (10 g/h)

alpha = v0/m0   #  1/h
tau = 1/alpha   #  h
cat(sprintf("alpha = %.3f h^(-1); tau = %.1f h \n", alpha, tau))

par(mfrow=c(2,1)) # 2 plot insieme
old.mar = par("mar")
par(mar=c(4.0,4.0,0.1,1.0))   # margini ottimizzati

tM = 300
curve(m0*exp(t/tau), 0, tM, xname='t', xlab='t [h]', ylab='m(t) [kg]',
      col='blue', lwd=2, xaxs="i")
curve(v0*exp(t/tau), 0, tM, xname='t', xlab='t [h]', ylab='v(t) [kg/h]',
      col='red', lwd=2, xaxs="i")

par(mar=old.mar)
par(mfrow=c(1,1)) # resettiamo il parametro mfrow