# w = 200
# e = exp(1)
# v = 1/60
# T = 50

# def law(e,w,v,t):
#     return 50+e**(w*v*((t**2)/2))

# def cooling(e,w,v,t,T):
#     ts = odeint(partial(law,v,T),T,t)
#     print('temp',ts)
#     print('time',t)
#     return ts

# time = np.linspace(0,20,21)
# a = cooling(e,w,v,time,T)
# pyplot.plot(time, a)
# pyplot.show()