run1

python
import numpy
import matplotlib.pyplot as plt

a = numpy.loadtxt('samples/120.raw',skiprows=1,delimiter=',')

x = a[:,0]

y = a[:,1]

num = int(len(a)*0.07)

xcrop, ycrop = x[0:num],y[0:num]
#shortening the data set to the Modulus Data

xfit,yfit = numpy.polyfit(xcrop,ycrop,1)
#fitting a first-order polynomial to the Modulus data

linfit = numpy.poly1d((xfit,yfit))

YMod = -numpy.diff(linfit)
#differentiating the polynomial provides the Modulus

print(YMod)