ratprint:false; wxplot_size:[500,400]; kill(y,t,r,f,c,%omega,m); eq:y(t)=-r*c*'diff(y(t),t)+m*sin(%omega*t); /* this expression must be subtracted at time zero to set the equation's initial value */ init_val:-(c*m*r*(%e^-(t/r*c))*%omega)/(c^2*r^2*%omega^2+1); assume(%omega > 0); atvalue(y(t),t=0,init_val); sol:desolve(eq,y(t)); f(t,r,c,%omega,m) := ev(rhs(sol),fullratsimp,factor); display(f(t,r,c,%omega,m)); f:440; %omega:2*%pi*f; r:1000; c:0.5e-6; m:1; wxplot2d( [sin(%omega*t), f(t,r,c,%omega,m)], [t,0,.005], [gnuplot_preamble, "unset key;set title 'f = 440 Hz, r = 1000, c = 0.5 uf'"]);