
function xdot=hw1_p2_circuit(t,x)

u = sign(t)*15;
xdot = -x-1.5*x^3+u;

