# Need help with fixing Python code to make output match with the corresponding graph

from numpy import zeros, linspace

import matplotlib.pyplot as plt

from scipy import stats

import math

dt $=5$

N$\_$t $=7$

B $=[0\mathrm{.}015,0\mathrm{.}021,0\mathrm{.}031,0\mathrm{.}04,0\mathrm{.}055,0\mathrm{.}075,0\mathrm{.}106]$

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

Bexp $=$ zeros $($N$\_$t$)$

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

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

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

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

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

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

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

tdouble $=$ math.log$(2)/$math$.$log$(1+$slope$)*$dt

mu$=$ math.log$(2)/$tdouble

Bexp$[0]=$ B$[0]$

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

Bexp$[$i$]=$ B$[0]*$ math.exp$($mu$*$t$\_$actual$[$i$])$

plt$.$figure$(2)$

plt$.$plot$($t$\_$actual, B$,"$bo$",$ t$\_$actual, B$[0]*$pow$($slope$+1,$t$),$ t$\_$actual, Bexp, $"$ro$")$

plt$.$legend$(["$data$","$B$0*(1+$c$)^($t$/5)","$B$0*$exp$($mu$*$t$)"],$ loc$=$ "upper left"$)$

plt$.$xlabel$("$time$")$

plt$.$ylabel$("$B$")$

dt$\_7=7$

t$\_$actual$\_7=$ linspace$(0,($N$\_$t$-1)*$dt$\_7,$ N$\_$t$)$

Bexp$\_7=$ zeros$($N$\_$t$)$

Bexp$\_7[0]=$ B$[0]$

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

Bexp$\_7[$i$]=$ Bexp$\_7[$i$-1]*$ math.exp$($mu$*$t$\_$actual$[$i$])$

plt$.$plot$($t$\_$actual, B$,"$bo$",$ t$\_$actual, B$[0]*$pow$($slope$+1,$t$),"$r$-",$ t$\_$actual$\_7,$ B$[0]*2\mathrm{.}718**($mu$*$t$\_$actual$\_7),"$gs$")$

plt$.$legend$(["$data$","$B$0*(1+$c$)^($t$/5)","$B$[0]*2\mathrm{.}718**($mu$*$t$\_$actual$\_7)"],$ loc$=$"upper left"$)$

plt$.$xlabel$("$time$")$

plt$.$ylabel$("$B$")$