The charge on the capacitor $q(t)$ as a function of time can be computed as
$q(t) = q_0e^{-Rt/(2L)}cos(\sqrt{\frac{1}{LC} - (\frac{R}{2L})^2}t)$
where $t = $ time, $q_0 = $ the initial charge, $R = $ the resistance, $L = $ inductance, and $C = $ capacitance.
Use MATLAB to generate a plot of this function from $t = 0$ to 0.8, given that $q_0 = 10$, $R = 60$, $L = 9$, and $C = 0.00005.$