/* wxMaxima 0.7.2 http://wxmaxima.sourceforge.net
Maxima 5.12.0 http://maxima.sourceforge.net
Using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (aka GCL)
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
This is a development version of Maxima. The function bug_report()
provides bug reporting information.

(%i1) diff(x,t,2);

Result

(%i2) EdoOscAmort: m*'diff(x,t,2)+b*'diff(x,t)+k*x=0;

Result

(%i3) assume(k>0);

Result

(%i4) solEdoOscAmort: ode2(EdoOscAmort,x,t);

Result

(%i5) rhs(solEdoOscAmort);

Result

(%i6) ev(rhs(solEdoOscAmort),%k1=0,%k2=xm);

Result

(%i7) solxOscAmort(xm,k,m,b,t):=ev(rhs(solEdoOscAmort),%k1=0,%k2=xm);

Result

(%i8) wxplot2d(solxOscAmort(10,1,4,1,t),[t,0,30]);

Result

(%i9) wxplot2d(solxOscAmort(10,1,4,0.4,t),[t,0,120]);

Result

(%i10) wxplot2d(solxOscAmort(10,1,4,0.4,t),[t,0,30]);

Result

(%i11) grafPosicao(tf,b):=wxplot2d(solxOscAmort(10,1,4,b,t),[t,0,tf]);

Result

(%i12) grafPosicao(30,1);

Result

(%i13) grafPosicao(30,0.3);

Result

(%i14) for b:0 step 0.4 thru 2 do (display(b),grafPosicao(30,b));

Result

(%i20) makelist(b^2,b,[0.,0.6,0.8,1.1,1.6,2]);

Result

(%i21) makelist(solxOscAmort(10,1,4,b,t),b,[0,0.6,1.4,2]);

Result

(%i22) wxplot2d(makelist(solxOscAmort(10,1,4,b,t),b,[0,0.6,1.4,2]),[t,0,30],[legend," "]);

Result

(%i23) grafPosicao2(tf,listab):=wxplot2d(makelist(solxOscAmort(10,1,4,b,t),b,listab),[t,0,tf],[legend,"
  "]);

Result

(%i24) grafPosicao2(30,[0.,0.4,0.8,1.2,1.6,2]);

Result

(%i25) solxb(b,t):=solxOscAmort(10,1,4,b,t);

Result

(%i26) map(lambda([b],solxOscAmort(10,1,4,b,t)),[0,0.6,1.4,2]);

Result

(%i27) grafPosicao3(tf,listab):=wxplot2d(map(lambda([b],solxOscAmort(10,1,4,b,t)),listab),[t,0,tf],[legend,"
  "]);

Result

(%i28) grafPosicao3(30,[0.,0.4,0.8,1.2,1.6,2]);

Result

(%i29) diff(solxOscAmort(xm,k,m,b,t),t);

Result

(%i30) solvOscAmort(xm,k,m,b,t):=diff(solxOscAmort(xm,k,m,b,t),t);

Result

(%i31) grafVelocidade3(tf,listab):=wxplot2d(map(lambda([b],solvOscAmort(10,1,4,b,t)),listab),[t,0,tf],[legend,"
  "]);

Result

(%i32) grafVelocidade3(30,[0.,0.4,0.8,1.2,1.6,2]);

Result

(%i33) grafK3(tf,listab):=wxplot2d(map(lambda([b],(4/2)*solvOscAmort(10,1,4,b,t)^2),listab),[t,0,tf],[legend,"
  "]);

Result

(%i34) grafK3(30,[0.,0.4,0.8,1.2,1.6,2]);

Result

(%i35)


Created with wxMaxima.