**TI86**
 Program file dated 06/23/04, 19:55      
 I
 9
Diff2Eqg9
7
 AxesOnÖFuncÖ{2,3}dim CÖ{2,2}dim C»ÖLbl TÖ[0,0]VEC:[0,0]COFÖreal (0*C)C:0*C»C»Öreal 0ra:rarb:real 0c1:c1c2Ö0t:0f1ÖClLCDÖ"aQ''+bQ'+cQ=0"s0ÖOutpt(1,1,"2nd DIFFEQ.W/DAMPING"+s0ÖOutpt(3,4,"UNKNOWN(S)        F1 Q(t) EQ."ÖOutpt(5,1,"F2 SERIES LC CIRCUIT    W/ RESISTANCE"ÖMenu(1,"F1",A1,2,"F2",B1,5,QU,QÖLbl A1ÖClLCDÖDisp "USE INIT. CONDS(Y/N)?"ÖZGETKEYÖ0f:0f1ÖIf K102Ö1fÖIf f==1:Goto A2ÖZLCINITÖLbl A2ÖClLCDÖDisp s0:Disp ""ÖInput "a? ",a0ÖInput "b0? ",b0ÖInput "c? ",c0ÖLbl T8Öb0/a0b0:c0/a0c0Öb0-4c0chÖIf ch<0:Goto IMAGINÖIf ch==0:Goto T4Ö(b0+ch)/2raÖ(b0-ch)/2rbÖIf f==1ÖThenÖ1c1:1c2Öy1=c1*e^(ra*x)+c2*e^(rb*x)ÖIf ra==0Öy1=c1+c2*e^(rb*x)ÖEndÖIf f==1:Goto AXÖe^(ra*t)C(1,1)Öe^(rb*t)C(1,2)ÖqC(1,3)Öra*e^(ra*t)C(2,1)Örb*e^(rb*t)C(2,2)ÖiC(2,3)ÖFor(x,1,2)ÖIf C(x,1)0ÖThenÖ1/C(x,1)dÖFor(j,1,3)Öd*C(x,j)C(x,j)ÖEndÖEndÖEndÖ(C(1,3)-C(2,3))/(C(1,2)-C(2,2))c2ÖC(1,3)-c2*C(1,2)c1Öy1=c1*e^(ra*x)+c2*e^(rb*x)ÖGoto AXÖLbl T4Öb0/2ra:rarbÖ1c1:1c2ÖIf f1ÖThenÖqc2Ö(i-ra*c2)/rbc1ÖEndÖy1=(c1*x+c2)e^(ra*x)ÖGoto AXÖLbl IMAGINÖb0/2»Öch/2ra:rarbÖIf f1ÖGoto T5ÖIf b0==0ÖThenÖy1=c1*cos (ra*x)+c2*sin (ra*x)ÖElseÖy1=c1*e^(b0/2*x)cos (ra*x)+c2*e^(b0/2*x)sin (ra*x)ÖEndÖGoto AXÖLbl T5ÖIf q0:Goto T7Öi/rac2Öy1=c2*e^(t+»*x)sin (ra*x)ÖGoto AXÖLbl T7ÖIf t0:Goto T9Öqc1Ö(q*»)/rac2Öy1=c1*e^(b0/2*x)cos (ra*x)+c2*e^(b0/2*x)sin (ra*x)ÖGoto AXÖLbl T9Öch/2ra:rarbÖ1C»(1,1)Ö1C»(1,2)Ö[q,i]VECÖe^((»+ra)t)(»+ra)C»(2,1)Öe^((»+rb)t)(»+rb)C»(2,2)Ösimult(C»,VEC)COFÖCOF(1)c1ÖCOF(2)c2Öy1=c1*e^((»+ra)x)+c2*e^((»+rb)x)ÖClLCD:EqSt(y1,YÖIf f11:ThenÖOutpt(1,1,"y1="+YÖElse:Outpt(1,1,"Q(t)="+YÖEndÖOutpt(3,1,"c1=":Outpt(3,4,round(c1,3)ÖOutpt(4,1,"c2=":Outpt(4,4,round(c2,3)ÖOutpt(5,1,"ra=":Outpt(6,1,"rb="ÖOutpt(5,4,round(ra,3):Outpt(6,4,round(rb,3)ÖOutpt(5,14,"»=":Outpt(5,16,real »ÖOutpt(6,12,"PER=":Outpt(6,16,imag round(2Ä/ra,2)ÖGoto T6ÖLbl AXÖClLCD:EqSt(y1,YÖIf f11:ThenÖOutpt(1,1,"y1="+YÖElse:Outpt(1,1,"Q(t)="+YÖEndÖIf f1ÖThenÖOutpt(4,1,"c1=":Outpt(4,4,round(c1,3)ÖOutpt(4,11,"c2=":Outpt(4,15,round(c2,3)ÖOutpt(5,12,"t=":Outpt(5,15,round(t,3)ÖEndÖOutpt(5,1,"ra=":Outpt(6,1,"rb="ÖOutpt(5,4,round(ra,3):Outpt(6,4,round(rb,3)ÖIf ch<0:ThenÖOutpt(6,12,"»=":Outpt(6,15,»ÖOutpt(7,1,"PER=":Outpt(7,5,round(2Ä/ra,2)ÖOutpt(7,11,"b0/2=":Outpt(7,17,round(b0/2,2)ÖEndÖMenu(1,MO,TX,2,PR,T,5,QU,QÖLbl TXÖClLCDÖLbl T6ÖZASKGRAFÖIf K102:Goto TÖZStdÖTraceÖGoto TÖLbl B1ÖZLCINITÖ0f:1f1ÖClLCDÖOutpt(1,1,"LQ''+RQ'+Q/C=0, L IS HENRYS, R OHMS, C IS"ÖOutpt(3,1,"FARADS. Q IS CHARGE  ON CAPACITOR."Ö4lf:ZLINFEEDÖInput "L? ",a0ÖInput "R? ",b0ÖInput "C? ",xÖ1/xc0ÖGoto T8ÖLbl QÖZMAKÖDelVar(CÖDelVar(C»ÖÐ@