def discovery(t, a, c):
sqt3 = sqrt( 3.0 )
halfsqt3 = sqt3 / 2.0
arg = a * halfsqt3 * t
return (c / (3.0 * a * a)) * (exp(-a*t) - exp(a*t/2.0) * (cos(arg) - sqt3 * sin(arg)))
@interact
def f(a=(1..100), c=(1..100)):
var("t")
show(plot(discovery(t,a/100,c),0,10,rgbcolor=(0.0,1.0,0.0)))
