/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 18.10.1 ] */
/* [wxMaxima: input   start ] */
kill(all);reset();
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]
proste operacje arytmetyczne
   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
suma:1+1;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
dzialania:1+2-3/2*4^3;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
1/3,numer;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
(1/3)/(1/7);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
0/1;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]
liczba urojna
   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
(-2)^(1/2);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
%^2;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
(-1)^(n);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sqrt(2);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]
równanie liniowe
   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
eq1:a*x+b-c=d;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sol:solve(eq1,x);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sol;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sol[1];
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x:rhs(sol[1]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
kill(x);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
eq2:a*x^2+b*x+c;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sol2:solve(eq2,x);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x1:rhs(sol2[1]);
x2:rhs(sol2[2]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
a:-1$ b:2$ c:-10$
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x1,ev;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
ev(x2);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]

równanie 3 stopnia

   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
kill(a,b,c);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
eq3:a*x^3+b*x^2+c*x+d;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sol3:solve(eq3,x);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x1:rhs(sol3[1])$
x2:rhs(sol3[2])$
x3:rhs(sol3[3]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
a:3$ b:-1$  c:10$  d:2;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x1,ev;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x1,float;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x2,numer;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x3,float;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
eq3,ev;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
wxplot2d(eq3,[x,-1,1]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]



   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
kill(x,a,b,c,d);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
eq4:a*x^4+b*x^3+c*x^3+d*x+e;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sol4:solve(eq4,x)$
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x1:rhs(sol4[1])$
x2:rhs(sol4[2])$
x3:rhs(sol4[3])$
x4:rhs(sol4[4])$
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
a:3$ b:-1$  c:10$  d:2$ e:10;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x1,ev,numer;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x2,ev,numer;
x3,ev,numer;
x4,numer;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
eq4,ev;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
wxplot2d([eq4],[x,-4,4],[y,-20,50]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]

równania  stopnia 5 i  większych --  tylko numeryczne 

   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
eq5:x^5+x^2+10;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
algsys([eq5],[x]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
wxplot2d([eq5],[x,-3,2]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]

układy  równań liniowych

   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
kill(a,b,c,d,e);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
eq6:a*x+b*y=d;
eq7:d*x-a*y=0;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sol6_7:solve([eq6,eq7],[x,y]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x:rhs(sol6_7[1][1]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
y:rhs(sol6_7[1][2]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
kill(x,y);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sol6_7:linsolve([eq6,eq7],[x,y]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sol6_7[2];
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]

układ 2 równań  nieliniowych

   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
eq8:y^2+x^2+c=0;
eq9:a*x^2+b*x+c=0;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sol9:solve([eq8,eq9],[x,y]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x1:rhs(sol9[1][1]);
y1:rhs(sol9[1][2]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x2:rhs(sol9[2][1]);
y2:rhs(sol9[2][2]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x3:rhs(sol9[3][1]);
y3:rhs(sol9[3][2]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x4:rhs(sol9[4][1]);
y4:rhs(sol9[4][2]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]

Rachunek różniczkowy

   [wxMaxima: comment end   ] */


/* [wxMaxima: comment start ]
pochodna z wyrażenia
   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
w1:x^2*sin(x^2);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
p_w1:diff(w1,x,1);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]
pochodna z funkcji 
 definicja funkcji
   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
y(x):=x^2*sin(x^2);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
y(2);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]
pochodna
   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
p_yx1:diff(y(x),x);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]
definicja  pochodnej  jako  funkcji
   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
py(x):=''(diff(y(x),x,2));
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
py(2);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
wxplot2d([y(x),p_yx1,py(x)],[x,-%pi/2,%pi/2],[y,-10,10]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]




   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
kill(all)$reset()$
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
y1(x):=sin(1/x);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
wxplot2d([y1(x)],[x,-1,1]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
y1dot(x):=''(diff(y1(x),x));
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
y1ddot(x):=''(diff(y1(x),x,2));
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
wxplot2d([y1(x),y1dot(x),y1ddot(x)],[x,-1,1],[y,-20,20]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]

 Rachunek całkowy

   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
kill(all)$reset();
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
y(x):=x^n;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
ydot(x):=''(diff(y(x),x));
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
assume(n>1);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
yc(x):=''(integrate(ydot(x),x));
/* [wxMaxima: input   end   ] */
/* [wxMaxima: answer  start ] */
yc(x):=''(integrate(ydot(x),x));y;
/* [wxMaxima: answer  end   ] */
/* [wxMaxima: answer  start ] */
n;
/* [wxMaxima: answer  end   ] */


/* [wxMaxima: input   start ] */
yc(x):=''(integrate(ydot(x),x));
/* [wxMaxima: input   end   ] */
/* [wxMaxima: answer  start ] */
y;
/* [wxMaxima: answer  end   ] */


/* [wxMaxima: input   start ] */
wxplot2d([x^2,x^0],[x,-1,1],[y,-1,2]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
'integrate(sin(x)*tan(x),x);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
integrate(exp(-1/x^2),x);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
integrate(x^2,x,0,10);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
wxplot2d([x^2],[x,-10,10]);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]



   [wxMaxima: comment end   ] */


/* [wxMaxima: caption start ]

   [wxMaxima: caption end   ] */
/* [wxMaxima: image   start ]
png
iVBORw0KGgoAAAANSUhEUgAAAb8AAAC3CAIAAAANVcf8AAAAA3NCSVQICAjb4U/gAAAetUlEQVR4nO3df1QTZ74/8Kd7rruJbn+EijkNu9RbRFaQBdGaG8GjtoAXiIXa74ogFruwiNatVPtLA+6tiK0/anGvRUqhW6sY456iVEgWRWMXYjZaqVwKfpFNVXaTNopJbVeTene/fv+YOkyTEEIgmUzyfh3/mJkMwwd6+uaZZ555nvvu3r1LAABghH7EdgEAAJyE9AQA8ATSEwDAE0hPAABPID0BADyB9AQA8ATSEwDAE0hPAABPID0BADyB9AQA8ATSEwDAE0hPAABPID0BADyB9AQA8ATSEwDAE0hPAABPID0BADyB9AQA8ATSEwDAE0hPAHDOZrP9+tcFlZWV1K7RaNTr9dS2Xq8XiyVdXV3sVcc+pCcAONfa2trd3ZOfn0/t1tbW5ebmUdsRERFFRYUyWZnNZmOvQJYhPQHACZvNVlNTW1RUKBAInJ6wePFik8nU2trq48L8B9ITIOjYbDa9Xk/fhlMsFoter7dYLNSuVqs1mUxpaWn0+Tdv3iSE6PV6o9FICBEIBElJiQ0NR3xevr9AegIEHR6PV1W1Nzc3r6mpiT64ZUtFbm6e1WqldpuammNiokUiESHEYDDk5ua1t2sIIbm5ebW1ddQ5UmlGd3cPFaZBCOkJEIzWr19HCCkvr6CyT6FQtLdryspkVFzabLb2ds3ChanUyRERETqdNiMjnRCi02k3bSqjjkdFRRFCOjo6WPkRWIf0BAhGIpGorExGCCktLevq6tq1qzIpKVEqlVKfGgwGQsjkyZOHvQghpKPjM+/W6q+QngBBSiqVJiUldnf3FBYWCYXC0lIZ/dHFixcJIRMnThz2IlSDNDghPQGCF3X/TghZtixnqGfrw2puVo5dRVyC9AQIXnSXZX29PJhHbnoG6QkQpIxGY3l5hVAoLCoqNJlM27fvoD+aMGGC+9cJ2pt3pCdAMLLZbKWlZYSQF19cW1BQkJSU2NysVKvV1Kfh4eGEkP7+fqdfyNwN2tt2gvQE8KUFkjlsl/C96urq7u6enJzsBQsWEEJWr15FCHn77d3UaPmwsDBCyBdffOH4hTJZqU6no7ap0U4JCTN8VrZfQXoCBJ2uri65XCEUCouLi6kjERER69aVmEymffv2EUJ4PF5SUqJGc4b5VVFRUwkh7e2alpbj1JHe3l5CSEJCgk+r9xv33b17l+0aAILFAskctfbM8Of5AbVa/dprG48c+Yga1Emx2WwGg4HP51MH169/yWKxvP9+HXtlsgltTwBwQiKRCIVClUrFPMjj8SIiIqjoNBqN7e2a5cvzWCqQfUhPAHCCx+MVFRXW1NTS84bYUalUQqFQIpH4uDD/gfQEAOeSk5OFQmFDQ4PjR0ajsaamVibbwOPxfF+Yn/g3tgsAAD/F4/H27fuD2Wx2/IjP5x88eCAiIsL3VfkPpCeA76Skp7FdwsgIBAKnb3AOdTyo4Jk7AIAn0PYEgFHR6XRarVYikcTFxQVVNyieGgH4jrKpme0Sxt7169flcsULL5TMm7egsrJSp9MFyYQjaHsC+M72ioqoab+wOxgWFhYwTTa5XCGXKwghOTnZAd8aRb8ngO8k/cec//XC/3ExMdHDzgPv2gMP3D9lyhTPvraj4zMXc4UEcIwiPQF8Z+5/zHH/sfvNmzephdgCw7p1JdnZ2WxXMZaQngC+4/v33O2WHfZAf3//rVu3XJzguu0ZExO9cGHq3Llzme/LBwb0ewIEstEPaHfnCo7pGcChSUN6AsCYCYbQpCE9AWC0kpISpdKM+Pj4oHoBCf2eAACewGh5AABPID0BADyB9ATwnSVZWUNNNgycg/QE8J34mTMtzqbLBC5CegIAeALpCQDgCaQnAIAnkJ4AAJ5AegIAeALpCeA7KakLBSEhbFcBYwNvavqIzWbTarXUTF/Tpk0L8qVcAQIAZgnxBaPRWFy82mQy0UcCb6ZYgGCDO3dfKC0tM5lMOTnZtbU1b765VSgU7tpVqdPp2K4LfO2PCsWXRiPbVcDYQNvT63Q6XXd3T05OdklJCXUkPDw8NzevpeW4WCxmtzbwsb5Ll2bOmsV2FTA2kJ5e9/nnnxNCnnzySfpIREREWZksNDSUvaIAYLSQnl5nMBgJIePHj2celEqlLJUDAGMD/Z4AAJ5AenrdAw/cTwi5ffs28+DmzeWbN5ezVBEAjAGkp9fFxcURQk6ePEkf0ev1zc3KK1eusFYTAIwa+j29TiKRCIVCuVxx//33z58/f2BgoKLiDULI4sVP251psVjMDpM/2q2mTQ8aTUhImDFjhjcLBwBX8K6RL3R1dclkZczR8klJiW+9tdPuNLVa/dprG928ZkZG+qZNZWNWIgCMENqevhAbG7tv3x80Gk1Hx2dhYaLp06c7HekZHh7u+9r8mUKh6O29xHYVY6mwsKCtrY27P1RU1FS8I0dDevqIQCCQSqUeD1TKyEinNs6fP//VVyZy72FUYOvtvdTcrGS7irGUmpJ8seei6k8tbBcCYwDp6UciIiJ0Oq3rczZvLqcCZcqUKT4pCsbS0Y8++he5j+0qYGwgPf2dWCyhNg4ePICZmQD8B9ITuKGsTBYAL2htLS9fmpu7Y+cOtgsZAb1en5ubx3YV/ogz4z2r91YbDAa2qwAA+B430lOr1e7ctm1eYtLzq1drtVqr1cr8dMrkf6f+sVUeAAQhbqRn54VOaqNFqVqekxs7Lbp6b3VfXx+7VQFAMONGehavKladOP7Sq6/SR3Zu25aWkpqTnd3a2spiYQAQtDjz1CgyMjIyMjJ/Rf6FCxcO7N/folQRQs7pzp7TnWW7NAB3YVW4QMKNtqfVas3Jzm48etRqtUokkneqqs52nH99S/nj4tlsl+Y7AwMDbJcAo/W4eLZAIGC7Chgb3EhPjUZzTnd2fcmLsxNmUj2eISEhy/Ly5ArFR41H2a7Ou4RCIbVx/fp1disBACZupCchJHzyo9QG3eNJPXyn5n8LYLNmzWS7BBgzWBUukHAjPZOTk5tVqv3ygwvT06gj53Rnl+fkZqSlVe+tZrc2APf1XbpkN94OuIsbT43MZnNvb298fLxEIunr6zvZenLntm2EkP4rV6mNkAcfHDdu3NZy7s3WHj8jwfUJV69cpTauMSa4AwDWcSM9VUrl70rLwic/uiR76ZLsJcWrivNX5B9vaTkkl1PP3MeNG/fjceNOKFVsVzpiw9b83b0JWKnV5Wh/0f7lR6zON3Hr1j/6Lnl3prWO8+fpbYvF4tXvBTBS3EjPutpacq+luXPbtqXLcjOk0sysrMysrM7Ozmcys9gukAWnTp7UnDrFdhXe9R1j6u7/vXOHxUoAHHGj37NZpaqufY8en3So/uDynNwn5s9vPHr05z//Obu1+VgwTOsJwAncaHvy+fzk5OTk5GS7Ts/1JS9SJ3z9zTc/+tF9jc3NXi3jb/39t27dHv68MfVRw5HPu7sJIWFhIoJpPQH8BjfS02w2d3R0JCQk0G8caTSauvfeo180+u7OHULIY16e/tLb13fq044OKj2/+fZb5vFp0THPOKwrFxhCQ0MfnvgwIaSqam9bu4Y6iFd0wN9wIz0PKw5T7c2Vq1clzZ0bHx/v2BQNVFFRU6nJ5L/55gfp+ejkR9OlGSwV5SMPPPggvT1u3DgWKxkrG8uwkF/g4Ea/5yen1dTGu1V76WGefX19kZGRxauK2a3N2yZMmMB2CQDgBDfS8/0PPtgvP7h0WS61Sz18p944qj9wgN3aANz3pdFos9nYrgLGBjfSk8/nR0VFbamoONtxnvnw/Zzu7O9KcSsEnPGHujojlkgIFNxIT4PBMDthZqlM1tvbm5ycLFcoPtG0v/Tqq/TL7wAAPsaN9Dz28THCGOZZf+AAn88vXlV86vTpgJ9jCQD8EzfSM39FPn3D3n/l6u9Ky6imqFarDfg5lgDAP3FjxJLT0fKH6g8eqj8YPDfvN2/eZLsEABjEjbZnX19f9d7qzs5OaohS18UeZlOU7eq8ix6x1H5v3DgA+ANutD3p9ubj4tnSRYvS0tODZ7R8eHg42yUAgBPcaHt2df0PtUENUZqdMJNa2P1nP/tZwI+Wh0DyXEGBKCyM7SpgbHCj7flOVZXBYDitVtfV1lK36i1KVYtSRc34yXZ13hUSElJWJmO7Chgbj4hEbJcAY4Yb6Wm1Wqll4Jbl5Wm12va2tner9hLG3PIBTCAQSKVSx+N4iATALm7cuV+4cCEjLa3+wAGz2SyRSF5+5RW7l46CEB4icRFWhQsk3Gh7Hti/nxrm+bvSMmpieYlEQj04MhgM8xKT2C5wLG3e/IPVmagJliAw9F26NHPWLLargLHBjfRMTU0137hBzeZJD/MsKCxMS08PC7g+eMQlACdw4849MytLrlB81Hh05epV1BHmG0fs1saWmJhotksACGrcaHsSQsxmc1xcXFxcXEFhoUqppB++H6o/yHZpY8zxCXtoaOjEiRPtDvL5fF9VBABOcCM9zWbz7ISZVI9nfHw8/fC9uakp8NLT6RN2APA33LhzVymV5N4cS/TE8hKJhJrxk+3qACAYcSM9p8fGMns86YnlG49iejrgksipU9HlEjC4cedO93h2dHTQS2me052l19QE4IRfZWezXQKMGW60PQkhBoMhJCQEE8sDgJ/gRnr29fXNS0x6fvXq1tZWq9UaFhZWvKq4WaViLhUHAOBL3EjPk60nCSEtSlVx4W9ip0Xv2L69s7OTz+dTD47Yrg4AghE30jN/Rf5++UH6rfZ3q/Y+k5lFLUdsNpvZrQ0AghM3nhpRzUyJRMKcDpl6aoQViSEg6fX6ixcvJicna7XatrZ2QsjcuUkSicRqtWo0mo6Oz+gjPB6P7WKDFDfSU6vVyjZsoF5sL15VnL8iX6PR0A/fAQLPxYsXy8sr1OrTfX1/nTVr5qefnm9uViYlJfb1/ZUQwjzy1ls72S42SHEjPZlzLK1cvSp14cLgWZkDgtzhw4d4PJ7RaHz66Wfa2zVJSYkVFVuYRywWi0AgYLvMYMSNfs/U1FTHTs/Go0cffvhhrMwBAWzJkl9RN+YikSgjI50QsmJFvt0RdP2zhRttz8ysrMysrM7OTu0ZLbPTkxDy0quvsl0dgLc4zg4zfvx4VioBR9xoe1qtVqvVGhcXRy9HvDA9jfoIt+3eZrPZbDYb21UA+B1upKdGo6EmBzEYDHw+Pzk5+Z2qKj9/48hisYwydGw22ydq9SgvsrW8XNPWNpqLnGo9+ezSpX9SKj/v6hpNJQABhhvpeaShgZocxPGNo1OnT7NdnRM2m23lc8/96qnMreVbPE6uY42N/7VRJk1OaWw44tlFPlGrTyhVpa+8+n8WPeXZ+ASbzfZB7XvXTde2lW9Zv+a3ly9f9uAiAAGJG+n59OLF9K06/cYRNU8du4UN5VD9weuma//49tsTSuXKXxd4cIUvv/yyqnI3IeRf//pX5Y4dOq12pFf4+uuv37i3RNKtf/zjT0pPFvyQ19dfN12jtu/cuaNqbvbgIgABiRvpSd+qv76lnH74Ts9Tx25tjq5cvvzh+3X07n0eXWTX9h3M3UlC4UivcKKl5TtGi/V/Lnw20gasxWL5sLaOeeSP9Qe/0OtHWglAQOJGelqt1h3btw8MDCzLy5MrFKoTx+keTz8cMK8+deru/7tL7169fFnZNLIm2xd6/ad/+QvzyO63do0o+ywWC9V0pQ1cuy6vrx9RGQf2feh4cPvWN0Z0EfBMaGhoRkY6czLQqKipwx4BX+JGemo0GmqY5xPz59cfOEAN8zx1+jRznTg/4dheI4R8UPuexWJx/yKHDtovN9Lb03Oq9aT7Vzja0OB48MPaOvdbjl/o9Q0KhePx3p6ekf4xAA+IxeJNm8pEIhF9JDs7e9gj4EvcSE9CCHXDzlxKU6vVTp069eVXXmG7tB/Yu2eP48HrpmtOjzt1Tnf2hFLleHxHRcWXRqM7V/jSaHRMcIr7LUdmgkdFR+fmP8usZER/DAACEjfSk5oUmbphp44w1zhitzamL/R6ZvClZz5Fb59Qqtwc8VNXU0NvP5GaEiqcRO/+oc55JtphnhYVHb26ZC2962bL0S7BC4qKlq9YwazE/T8GAIGKG+lJCOnr64uMjLQbLU8NY/L2t/5Erd5aXu7Omcz2Wkp62m9LSpiJs6dy97B9l8qm5t6eHnq3cOXKlzdspHdPKFWfqNWur2CX4GtK1i7KzIyKHlz83Z2WIzPBxYlzHhfP5vF4mxi/hBNKlR/2OAP4Ejfe1DQYDGkpqY+LZy/NyaGnCDEYDMc+PnZYcaj/ylXhxImEEGbG9V+5yoyh0TMaDHuqXbVz7dprS3NzqcT5bdFK6khvT8+xxkYXK9tQgyvp3WcLCx4RiR4RiVLS0+grv7N7t9jlpGTMe/OU9LTpsbGEkFc2bijIW04f37tnz8ayIWf2+0StZv7qilZ937M8PTaWWcmON7Z+eOgQpkeDoMWNtuexj48RQs7pzq4veZEe6Umvz0EI+fG4cT8eN+6EUkX/G9voJIT0Xfy/rrsdme21xdnZj0VEkHuJQx+vqtzt4iLHGhvpwZWEkKzFi6mNVWvW0G3Y66Zr71W/O9QVzunOMn/wpbnfL1vyWETEYkZqu2g52my2d3YPPqxPSU+jfhCnlRxrbByqEoCAx42255LsJfff/9O62tr+K1cJITu3bdu5bdvC9LS85cvj4+N9U8OdO3d279r15k7ncyna3XHnMZ6xrFqz5sL583QsDnURuzFGL8tk9LRjAoFgReFvdtxbg6RBochYJGWGGsVmszET/NnCAuY5vyle2XZaTZcxVMvxVOtJZoKvWrOG+aldJVWVuydNmnTr1m3HH2dMXL1y1UtXBhg9bqRnSEjIsry8ZXl5Wq22va3t3aq9hJAWpapFqfLle+46zRllU3O6NMPuuOMdN3O+RbvE0WnOfKJWz1uwwO4izMGVocJJTyQ/yfw0XZrx59NqneYMtbt96xvVdbV2VzjVepKZ4HTTlcLj8Z5fu/a/Nsqo3euma/L6+ucKfvAelMVicfGDOK2EvqA3fHf37vAnAbCEG+lpNptv3LgRGRlJrc9RUFjY9uc/H5LLz+nOUq3RAYuZENI4wvcI+Xz+I26MlXvtpZfosNhRUTEjYYbdVzHba6HCSTnLltldwS5x3tm9+5fx8cxgshtc+fzatY6twqJVq+grUI/OmTlul+CrS9Y6Bt+8BQuYHZcf1tbNmz+f2T492tDg+gehrF23LvdeJQBBixv9niqlknopk5ofJCQkJDMrS65QfNR4lBrDZLV9Z7V991hExIj+uROdhJC169Yxd18v28R8dG6xWHYw1vVcUfgbpw9S6GcvhJDrpmt2r/HYDa50bJkSQh6LiGCOPbJ7dM7sMw0VTlqUmen0Z7G7E6/Zu5fethslOtQPQgh5RCSKmBrp9COA4MGNtmfTse+fGtEzIj+Z/GRkZGRcXFxcXJy3By09IhK9LJPREWn36Jz5Vo/jHTeNyj66Z7NBoVjw5BPUA3HHwZVDVbIoM/Pk8RP07Tn96Nyuz9RF8AkEAubPwuyLsBslOtQPQnkiJeUr45eEkOlxv3zgwQeZH8XPSHDxhU79YtovnB6vqtrb1q4Z6dUAfIMb6SlXKLRabXNT06H6g4SNp0Z2t95VlbunRUdPj439Qq9nttde3rDRxQieRZmZf5TL6Rbinsrdle/s4fF4zEc9Kelp9DQojng83pqStfQQqBNKVUrqwsfFs5kt2WGDL12a8fGRI3QEf1D7niRxjsVstktw10ORcvPy0jIyvL2cjl0uA/gVbty5E0IkEsmWioqzHefp+UFalCrqdSPfFPCqTMYc+r65rMxms9ndcbsIPkKI3YBzqg1r97CeHmM0lOmxscyxRzve2Hqxp4fZZzps8BFCXtm4gd6m3iJl3sJTw+NdX4EQgpXIIMhxJj0pISEh1BjP/fKD9OtGvvnWAoHg+bWD3Y7XTdd2vPEms73GjKShOA7/rGWMwLcbXDmUvPxnmYMu39wy2OvqZvDZdaGeUKp0jKdAzC5aABgKx9KTwufzJRIJvTiHz74v9cya3j11/Di97WbwkR8OOCeEWG7cYH7kzhXscryfMd+7+8G3KDOTWQbN/R8EIMhxMj1p1OtGvvyOdtlHG/aOm0YN/3Q87nRw5VDscpwyouCz60aguZngAMDt9PQ9u3Yfxe6tnmGlSzPEiXOYR0ImThxqcOVQ7Aa6Oz3i2vTY2CdSU5hHpE9noTcTwE2BkJ5/vXKZ+uebbzdvwYKomGjmkUmTnLRGXXv6mWeYu3PmJo10uo1HRKIJP/0pvTvj8VluDl9lSl64kN7m8fkuBksBgB1ujFjyN2kZGb3d3z8oj5s588mUFNfnOxJLJOmZTykbPyaEPBQieP6FFzwoo3b/h6sKCr82mx8KEZS9/roHV5DMmfPSxg1/lB+6evnyb9ete+ihhzy4iG+YTCY9llRiQ39/P9sl+Kn77uJVYo9oNZrNZZt+8pOfHDrS4NksbTab7fy5T2/evBn+aDg1bN4D1OtGhr//3eMr0Nfh8/n+Nt3c5s3lzc2eLAUKXpKRkb5p05BzGwYbtD09JElMVJ0awUJDjng8XuLcpFGWQXVTjr6zEt2dACMVCP2eAAC+h7Yn+K+oqKlDfXTr1q0zZ7R37txx5zqPPfbv06ZNG7u6hvHPf/6zs7Pzq69Mrk+Lipo6ZcoU35R048aNjo7Phv11DVuSi/8iQQj9nsA9Nptt9ernu7sH33CNiYlevjwvPDycEHL79u2TJ0/K5T9YTrm2tiZ2dF3D7tcmk5W2uze5SVFRYcEIx5l5QKfTvfBCiZsn+6akwID05ACLxaLRaEJDQ8ViMdu1+IWmpqby8sH3U9etK8l2WC1Kr9e/+OJ6k+n7BmBMTPT777u1IukoOT7pSkpKnD378QkTJphMJo3mDDP0CSFvvrl1gbMJCceKxWLJz3+O/j2QH/6lGRgYaGk5blewt0sKGEhPDtDr9bm5eXjcSXvqqSw6DpxGJ8VoND799OC42t//vtLbf366uroKCwfHzMbERMtkGyN++CaFWq1+7bXBdVKFQuHhw15cXM8uzZ02Le0ap94uKWDgqRFwTFdXFx2dQqEwc4h5oAkhIpGorGxw4RCtVuvt2o4cOUpvC4XCqqp3IhxeQluwYMHvf19J75pMptbWVi/VYzQah41OQohYLPZZSYEE6Qkc08NcLbmo0HUTKTExkd626wkdcxaLhRlVb7/91lC1icXidesG23oNDUe8VFJbWxu9HRMT7aJDUywW5+QMNuG9V1IgQXpyUldXV1NTk1qtZrsQFvT2XqK3ExKGmcdeIBBkZKTTu0aXa0qP0oULF+jtnJxsx1YnU2pqKr3d3d3DXGRlDDF/V8uX57k+ecmSJT4oKZAgPblHp9MVFhbV1NRGRUWxXQvLRG682h8WNniO1Wr1XjHXrg2u5CyRSFyfLBAIYhizJZjNZm+UxGwLU4+JXBCJREKh0NslBRKkJ8fQHfwy2QZ3sgO+/fZb33wjZkNv4sSJw54/efJkevv27dveKInJdVuYMmvWTG+XEUiQnlxCR6cPHh9zgjvzhvztb3+nt0NCQrxZjl9zp9fi5s2bPqgkYCA9OePTT89T0SkUCqdODd5XPhISZtDbp0+fdn2y0WikB64LhUKvvs7PLGxgYGDY869cuUJvjx8/3hslMTsHbjBWMXDKZrO5OcgfKEhPzjCZTEKhMCkp0WQy7du3j+1yWMN8UlRTU+v64cbhw4fp7czMRV4si5AJEybQ28OOjrJYLMxh815qFP/yl4OvV509e9b1yZ2dncxdd+70gxzSkzOEQmF1ddX69esIIXK5oquri+2K2CESiZhNqv/+7z1DnanT6ZijlObPn+/VwpiLY8vlCtd3yhrNYCsvKSnRS41i5sOrmppaFyXZbLaKijfoXeboJRgK0pMzZs2aKRKJRCIRNVTw7bcrbTYb20WxY+XKwfd5mpuV69e/5NgCbWpqYr4/k5GR7u3GlN3oqLfe2jXUfyCj0ch801QqzfBSSXFxcczH6KWlZUMFaHV1NfNtzkWLvNtODwxIT+7JzMwUCoXd3T2NjY1s18IOsVjMzKn2ds1//md6ZWUlNQZWoVA89VQWM54IIYWFvpj5YtmywcUB29s1MlmpY6x3dXUVF6+md4VC4bDDmzzG4/FkssGFsru7e4qLV9vdtVgsls2by5mNdB/8pQkMeM+dAxzfc6cfvh858lFwjlsa0VRGvhyiUFlZafdSU1FR4fTp0ydOnNjf39/W1u77KTnq6upqamqZR2JiohcuTJ0wYYJafdrud4iX3N2H9OQAo9FYW1sXFTWVOR1GXV2dwWBMSJghlUpZrM0HxGJvNc3AAzqd16cL4AqkJ/g7pKdfQXrS0O8JAOAJpCcAgCeQngAAnkC/JwCAJ9D2BADwBNITAMATSE8AAE8gPQEAPIH0BADwBNITAMATSE8AAE8gPQEAPIH0BADwBNITAMAT/x9w69GyWXg2iAAAAABJRU5ErkJggg==
   [wxMaxima: image   end   ] */


/* [wxMaxima: input   start ] */
kill(all)$ reset();
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
depends(x,t);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
eq:m*diff(x,t,2)-k*x-c*diff(x,t)=0;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
eq1:diff(x,t,2)+K*x+ω_0^2*diff(x,t)=0;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: comment start ]
K=(k/m),  ω_0:sqrt(c/m)=2*π*f_0
   [wxMaxima: comment end   ] */


/* [wxMaxima: input   start ] */
/*assume(ω_0>0,K>0);*/;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
sol1:ode2(eq1,x,t);
/* [wxMaxima: input   end   ] */
/* [wxMaxima: answer  start ] */
neg;
/* [wxMaxima: answer  end   ] */


/* [wxMaxima: input   start ] */
sol1p:ic2(sol1,t=0,x=10,'diff(x,t)=0);
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
ω_0:0$ K:1$
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
x:rhs(sol1p),ev;
/* [wxMaxima: input   end   ] */


/* [wxMaxima: input   start ] */
wxplot2d([x],[t,0,100]);
/* [wxMaxima: input   end   ] */



/* Old versions of Maxima abort on loading files that end in a comment. */
"Created with wxMaxima 18.10.1"$