from numpy import zeros, linspace

import matplotlib.pyplot as plt

from scipy import stats

N$\_$t $=6$

B $=[0\mathrm{.}022,0\mathrm{.}036,0\mathrm{.}06,0\mathrm{.}101,0\mathrm{.}169,0\mathrm{.}266]$

dB $=$ zeros$($N$\_$t $-1)$

t $=$ linspace$(0,$ N$\_$t$-1,$ N$\_$t$)$

for i in range$(0,$ N$\_$t$-1)$:

dB$[$i$]=$ B$[$i$+1]-$ B$[$i$]$

print$("$B $=",$ B$,"$

$","$dB $=",$ dB$)$

copy$\_$dB $=$ dB$[$:$($N$\_$t $-1)]$

copy$\_$B $=$ B$[$:$($N$\_$t $-1)]$

copy$\_$t $=$ t$[$:$($N$\_$t $-1)]$

slope, intercept, r$\_$value, p$\_$value, std$\_$err $=$ stats.linregress $($copy$\_$B$,$ copy$\_$dB$)$

print$(\text{'}$slope $=\text{'},\text{'}\%\mathrm{.}3$f$\text{'}\%$ slope, 'intercept $=\text{'},\text{'}\%\mathrm{.}4$f$\text{'}\%$ intercept$)$

plt$.$figure$(1)$

plt$.$plot$($copy$\_$t$,$ copy$\_$dB$,\text{'}$bo$\text{'})$

plt$.$xlabel$(\text{'}$t$\text{'})$; plt$.$ylabel$(\text{'}$dB$\text{'})$

plt$.$figure$(2)$

plt$.$plot$($copy$\_$B$,$ copy$\_$dB$,\text{'}$bo$\text{'})$

plt$.$xlabel$(\text{'}$B$\text{'})$; plt$.$ylabel$(\text{'}$dB$\text{'})$

plt$.$figure$(3)$

plt$.$plot$($t$,$ B$,\text{'}$bo$\text{'},$ t$,$ B$[0]*(1+$slope$)**$t$)$

plt$.$legend$([\text{'}$data$\text{'},$ 'equation'$],$ loc $=$ 'upper left'$)$

plt$.$xlabel$(\text{'}$t$\text{'})$; plt$.$ylabel$(\text{'}$dB$\text{'})$