def OnFeedData(self, evt):
     pipeInData = gen()
     self.nerfModel.SendPipe(pipeInData)
    x_1 = x_1 + (x_1_k1 + 2*x_1_k2 +2*x_1_k3 + x_1_k4) / 6.0  
    x_2 = x_2 + (x_2_k1 + 2*x_2_k2 +2*x_2_k3 + x_2_k4) / 6.0  

    #dx_0_plot.append((dx_0 + ddx_0) / 2.0)
    #dx_1_plot.append((dx_1 + ddx_1) / 2.0)
    #dx_2_plot.append((dx_2 + ddx_2) / 2.0)
    x_0_plot.append(x_0)
    x_1_plot.append(x_1)
    x_2_plot.append(x_2)
    ## End of Heun's Method
    Ia_fiber45 = max(GI*(lce_i-x_1-LSR0), 0.0)
    Ia_fiber_plot45.append(Ia_fiber45)
    return [x_0, x_1, x_2]

Lce = []
Lce = gen(SAMPLING_RATE)

#plot(Lce)

Ia_fiber=[]
Ia_fiber2=[]
Ia_fiber45=[]
x_0_plot = []
x_1_plot = []
x_2_plot = []
Ia_fiber_plot = []
Ia_fiber_plot2= []
Ia_fiber_plot45= []
dx_0_plot = []
dx_1_plot = []
dx_2_plot = []