home *** CD-ROM | disk | FTP | other *** search
/ 64'er 1985 August / 64er_Magazin_85-08_1985_Markt__Technik_de.d64 / newea2 (.txt) < prev    next >
Commodore BASIC  |  2022-10-26  |  11KB  |  287 lines
  1. 10 rem*********************************
  2. 20 rem*                               *
  3. 30 rem*             newea2            *
  4. 40 rem*                               *
  5. 50 rem*        von  hans haberl       *
  6. 60 rem*                               *
  7. 70 rem*********************************
  8. 100 rem  steuerzeichen:
  9. 101 rem   "[147]" = shift + clr/home
  10. 102 rem   "" = crsr down
  11. 103 rem   "[145]" = crsr up
  12. 104 rem   "[157]" = crsr back
  13. 105 rem   "" = ctrl + 9
  14. 106 rem   "[146]" = ctrl + 0
  15. 107 rem   "[159]" = ctrl + 4
  16. 108 rem   "[158]" = ctrl + 8
  17. 109 rem   "[150]" = comDrore + 3
  18. 110 rem   "[155]" = comDrore + 8
  19. 130 D6,6:print"[147][155]":D1,0,40,23,81,7:D0,1,38,25,81,7:D3,3,34,19,32
  20. 150 printD(11,9)"ac-netzwerkanalyse":printD(12,15)"von  hans haberl"
  21. 160 DT1,"goto300"+chr$(13)
  22. 170 dimi,j,r,s,f,u,v,x,y,p,q,k,n,m,o,c,l,j1,j2,kf
  23. 180 dimx1,x2,y1,y2,xf,xs,yf,ys,tl,tp,t0,z,d,e,w,c$,x$,n$,l%,q%,i%
  24. 190 dimr(20,21),s(20,21),m%(20,21),k(20),y(2,20,20),f(9,20),u(9,20),w(9,20)
  25. 200 dimi(90),j(90),b(90),a(90),x(2,3,5),n(10),p(11)
  26. 210 dimc$(9),d$(9),e$(9),k$(9),b$(9)
  27. 220 tl=10/log(10):tp=180/(NULL):t0=2*(NULL):fori=0to9:readc$(i),d$(i),e$(i):nexti
  28. 230 dataf,knoten ,anzahl der ,a,widerstaende,r,d,kapazitaeten,c
  29. 240 datao,induktivitaeten,l,n,transistoren,c b e,e,y-vierpole,1 1' 2 2'
  30. 250 datal,op's,+ - a,s,frequenz,frequenzfile,0,0,r,0,0,w
  31. 280 i=0:D :getc$:i=i+1:D)c$<>""ori=300:print"[147]":gosub500
  32. 300 close2:close4:close15
  33. 400 D::print">";:gosub5010
  34. 410 forc=0to7:ifc$<>c$(c)thennextc:gosub500:D6
  35. 420 print:onc+1gosub1000,1000,1500,2000,3000,3300,4050,4000:D6
  36. 500 print"menue:":print"n=netzwerk eingeben":print"e=elemente eingeben"
  37. 510 print"f="d$(7)"gang":print"a=ausgabe":print"d=diagramm"
  38. 520 print"l=daten laden":print"s=daten speichern":return
  39. 1000 c$="":gosub5650:gosub5700:ifq=0thenp=0:return
  40. 1010 i=q:ifc<1andu<9thenprintd$(7)"en von ";:gosub5700
  41. 1020 Df:ifq>0andf(q,0)>0thenforj=0tof(q,0)+1:f(i,j)=f(q,j):nextj
  42. 1030 Df:ifc=0thenprint"ohne aenderung ?";:gosub5010
  43. 1040 q=i:l=0:input"notiz";k$(q):ifc<1thenw(q,0)=c:ifc$<>"j"thenl=1
  44. 1080 gosub1400:ifp<3orc<0thenreturn
  45. 1100 c=1:D::print"file#[146] , n[146]ormierung oder d[146]rucker?";
  46. 1110 gosub5010:j=asc(c$)-48:l=0
  47. 1120 ifc$="n"thengosub1300:D6
  48. 1130 ifc$="d"thengosub7200:D6
  49. 1140 D5 ifj<1orj>9:D5 iff(j,0)=0
  50. 1150 gosub1200:D6:return
  51. 1200 gosub5900:p=1:l=1:D :f=f(q,p):iff(j,l)<fthenf=f(j,l)
  52. 1220 c$=str$(f):printright$(c$,len(c$)-1):print"[145]";
  53. 1230 iff=f(q,p)thenprinttab(7)"[159]";:x=u(q,p)
  54. 1240 Df:gosub5950:printtab(23):x=w(q,p):gosub5950:print:print"[145][155]";:p=p+1
  55. 1250 iff=f(j,l)thenprinttab(15):x=u(j,l)
  56. 1260 Df:gosub5950:printtab(31):x=w(j,l):gosub5950:l=l+1
  57. 1270 print:D)p>f(q,0)orl>f(j,0)
  58. 1280 printtab(10)"[159]file"q"[155]file"j:return
  59. 1300 input"bezugsfrequenz";f:input"db-wert";s:fori=1tof(q,0)
  60. 1310 iff=f(q,i)thenr=s-u(q,i):fori=1tof(q,0):u(q,i)=u(q,i)+r:nexti
  61. 1320 Df:gosub1400:return
  62. 1330 nexti:return
  63. 1400 gosub5900:forp=1to20
  64. 1410 iff(q,p)=0andc=0thenl=1
  65. 1420 ifl=1thenprintf(q,p);:inputf(q,p):print"[145][145]"
  66. 1430 iff(q,p)>0thenonc+2gosub6000,4400,4700:nextp
  67. 1440 f(q,0)=p-1:return
  68. 1500 gosub5750:q=p:ifq=0thenreturn
  69. 1600 gosub1700:fori=1toq:forj=1tof(p(i),0):x1=x2:y1=y2
  70. 1610 x2=log(f(p(i),j))*xf+xs:y2=u(p(i),j)*yf+ys
  71. 1620 ifj>1thenDx1,y1,x2,y2,1:Dx1+1,y1,x2+1,y2,1
  72. 1630 nextj:ifq>1thenDx2-18,y2-6,p(i)+48,1,1
  73. 1640 nexti:gosub5050:DG:print"ausgabe auf drucker ?";:gosub5010
  74. 1650 ifc$="j"thengosub7300
  75. 1660 return
  76. 1700 x1=f(p(1),1):x2=x1+1:y2=u(p(1),1):y1=y2-1
  77. 1710 fori=1toq:p=p(i):r=f(p,f(p,0)):s=f(p,1):ifr>x2thenx2=r
  78. 1720 ifs<x1thenx1=s
  79. 1730 forj=1tof(p,0):r=u(p,j):ifr>y2theny2=r:Dg:ifr<y1theny1=r
  80. 1740 nextj,i:yf=185/(y1-y2):ys=192-y1*yf:xf=288/log(x2/x1):xs=32-log(x1)*xf
  81. 1750 D6,4:s=(y2-y1)/10:r=10^int(log(s)/log(10)+.2):s=r*int(s/r+1)
  82. 1760 i=0:D :i=i+1:y=s*int(y1/s+i)
  83. 1770 ify<y2thenv=y*yf+ys:DO0,v-3,str$(y),1,1,8:D32,v,319,v,1
  84. 1780 D)y>y2
  85. 1800 ifx2/x1>1.1thengosub1820:Dg:gosub1860
  86. 1810 return
  87. 1820 s=log(x2/x1)/10:v=x2-x1:x=x1
  88. 1830 D :r=10^int(log(x*v/(x+v))/log(10)):x=r*int(x/r+1):gosub1900
  89. 1840 x=x*exp(s):D)x>x2:return
  90. 1860 s=(x2-x1)/6:r=10^int(log(s)/log(10)+.2):s=r*int(s/r+1)
  91. 1870 i=0:D :i=i+1:x=s*int(x1/s+i):gosub1900:D)x>x2:return
  92. 1900 u=x:D::D5 ifu<1000:u=u/1000:D6:c$=right$(str$(u),5)
  93. 1910 u=log(x)*xf+xs:DOu-4-4*len(c$),193,c$,1,1,8:Du,0,u,191,1:return
  94. 2000 c=-1:gosub1000:gosub5000:ifc$="j"orp<2thenreturn
  95. 2010 l%=p-1:c=6:gosub5660:ifu>0thenp(10)=p(1):Dg:p(10)=9-int(q/9)
  96. 2020 D :print"optimierung von max. 4 elementen"
  97. 2030 p=1:D :printp".element "b$(p);:inputb$(p):gosub6050:D)p>4
  98. 2040 gosub5000:D)c$<>"j":p(0)=9-p:ifp=9thenreturn
  99. 2050 p(9)=q:p=p(10):f(p,0)=l%:w(p,0)=1:k$(p)="ist":z=1:q%=10:print"x=abbruch"
  100. 2100 D::l=2:q=p(q%):forp=1tol%:f(q,p)=f(p(9),p):gosub4400:D5 ifr=0
  101. 2110 w(q,p)=u(q,p)-u(p(q%-1),p):nextp
  102. 2130 forj=10toq%:u=-r:r=0:forp=1tol%:r=r+w(p(q%),p)*w(p(j),p)/w(p(9),p)^2
  103. 2140 nextp,j:ifq%=11thens=(1+atn(u/r/6)/3)^2:gosub3500:q%=10:D6:z=0:return
  104. 2150 print:print"zyklus"z"fehler:"r:forp=1top(0):printb$(p),b(p(p));
  105. 2160 iflen(b$(p))>4thenprintb(p(p+4)):Dg:print
  106. 2170 nextp:getc$
  107. 2180 ifc$="x"orr<1thenprint"abbruch?";:gosub5010:D5 ifc$="j"orc$="x"
  108. 2200 q%=11:i%=Dr(z,p(0))+1:z=z+1:s=1.1:d=b(p(i%)):e=b(p(i%+4)):gosub3500
  109. 2210 D6:z=0:return
  110. 3000 j=j1:fori=1to5:p(i)=n(i):nexti:D :print""e$(0);d$(0)"(ohne 0)",n;
  111. 3010 inputn:fori=1to6:printe$(0);d$(i)tab(30)n(i);:inputn(i):n(i)=abs(n(i))
  112. 3020 nexti:gosub5000:gosub5600
  113. 3030 D)j1<91andn>0andn<21andn(6)<6andc$<>"j"
  114. 3040 m=1:j=j+j1:fori=1to5:p=m+p(i)*(1+int(i/4)):m=m+n(i)*(1+int(i/4))
  115. 3050 ifp>mthenfork=mtoj:gosub5620:nextk
  116. 3060 ifp<mthenfork=jtomstep-1:gosub5620:nextk
  117. 3070 nexti:D :k=1:fori=1to6:ifn(i)>0thenprint""d$(i)" an "d$(0)
  118. 3080 Df:forj=1ton(i):onigosub5100,5100,5100,5150,5200,5150:gosub5250:nextj
  119. 3090 nexti:D :print"ausgangs"d$(0),n(9);:inputn(9)
  120. 3100 D)n(9)>0andn(9)<=n
  121. 3110 gosub5000:D)c$<>"j":print
  122. 3300 D :k=0:j(0)=1:fori=1to6:ifn(i)>0ori=1thenprint:printd$(i)" "e$(i)
  123. 3310 Df:forj=sgn(i-1)ton(i):onigosub5300,5300,5300,5350,5450,5400
  124. 3320 Df:k=k+1+int(i/4):nextj
  125. 3330 nexti:gosub5000:D)c$<>"j"
  126. 3500 ifz>0thenb(p(i%))=d*s:b(p(i%+4))=e*s^sgn(49-p(i%+4))
  127. 3600 gosub6200:u=-1:m=n+1:fori=1ton:forj=1tom:m%(i,j)=0:nextj,i
  128. 3610 k=1:j2=1:forl=0to2:fori=1ton:forj=1ton:y(l,i,j)=0:nextj,i
  129. 3620 j2=j2+n(l+1):D::D5 ifk=j2:i=i(k):j=j(k):gosub5800:k=k+1
  130. 3630 D6:nextl:y(0,1,1)=y(0,1,1)+1/b(0):m%(1,m)=u
  131. 3650 j2=j2+2*n(4):D::D5 ifk=j2:l=0:s=abs(b(k))/.026:r=s*5e-4:p=i(k+1):o=i(k)
  132. 3660 i=p:j=o:gosub5820:j=j(k):r=s:gosub5820:o=j:r=s/a(k):gosub5820
  133. 3670 l=1:r=s/t0/b(k+1):gosub5820:i=i(k):p=i:r=r*5e-4:gosub5820:k=k+2:D6
  134. 3700 s=1:j2=j2+2*n(5):D::D5 ifk=j2:i=i(k):j=j(k):p=i:o=j:kf=0:gosub5850
  135. 3710 i=i(k+1):j=j(k+1):kf=1:gosub5850:p=i:o=j:kf=3:gosub5850
  136. 3720 i=i(k):j=j(k):kf=2:gosub5850:s=s+1:k=k+2:D6
  137. 3750 j2=j2+2*n(6):D::D5 ifk=j2:l=i(k+1):fori=0to2:forj=1ton:y(i,l,j)=0
  138. 3760 nextj,i:y(0,l,l)=1:m%(l,l)=u:m%(l,i(k))=u:m%(l,j(k))=u
  139. 3770 a(k)=t0*b(k+1)/b(k):k=k+2:D6:j2=j2-2*n(6)
  140. 3780 fori=1ton:m%(i,0)=i:nexti
  141. 3790 p=n(9):ifp<nthenq=n:gosub6300
  142. 3800 fori=1ton-1:forj=i+1ton:m%(i,j)=m%(i,j)orm%(j,i):m%(j,i)=m%(i,j):nextj,i
  143. 3810 ifn<3then3920
  144. 3820 forp=1ton-2:f=n*n:fork=pton-1:v=0:fori=pton-1
  145. 3830 ifm%(i,k)thenforj=i+1ton:w=m%(k,j)andnotm%(i,j):v=v-w:nextj
  146. 3840 nexti:v=2*v-m%(k,m):ifv=fthengosub6350:ifr>sthens=r:q=k
  147. 3850 ifv<fthenf=v:q=k:gosub6350:s=r
  148. 3860 nextk:ifq>pthengosub6300
  149. 3870 iff=0then3910
  150. 3880 fori=p+1ton-1
  151. 3890 ifm%(i,p)thenforj=i+1ton:m%(i,j)=m%(i,j)orm%(p,j):m%(j,i)=m%(i,j):nextj
  152. 3900 nexti
  153. 3910 nextp
  154. 3920 fori=1ton:m%(i,m)=m%(i,m)orm%(i-1,m):k(m%(i,0))=i:nexti:gosub6250
  155. 3930 ifz>0thenreturn
  156. 3950 print"ausgabe auf drucker?";:gosub5010:ifc$="j"thengosub7000
  157. 3960 return
  158. 4000 f=0:open15,8,15:print"schaltung speichern ?";:gosub5010
  159. 4010 ifc$="j"andn*j1>0thengosub4100
  160. 4020 print:iff=0thenprinte$(7)"s speichern ?";:gosub5010:ifc$="j"thengosub4300
  161. 4030 close2:close15:return
  162. 4050 open15,8,15:gosub5500
  163. 4060 iff=0theninput#2,c$:ifright$(c$,5)="newea"thenonval(c$)gosub4110,4300
  164. 4070 close15:close2:return
  165. 4100 gosub5500:iff>0thenreturn:Dg:print#2,"1newea":input"notiz";n$
  166. 4110 print:ifc=7thenprint#2,chr$(34)+n$+chr$(34):Dg:input#2,n$:printn$
  167. 4120 n(0)=n:fori=0to9:ifc=7thenprint#2,n(i):Dg:input#2,n(i)
  168. 4130 nexti:n=n(0):gosub5600:fori=0toj1
  169. 4140 ifc=7thenprint#2,i(i):print#2,j(i):print#2,b(i):print#2,a(i):Dg:input#2,i(i),j(i),b(i),a(i)
  170. 4150 nexti:ifn(5)=0then4190
  171. 4160 fork=1ton(5):fori=0to2:forj=0to3
  172. 4170 ifc=7thenprint#2,x(i,j,k):Dg:input#2,x(i,j,k)
  173. 4180 nextj:nexti:nextk
  174. 4190 gosub5550:close2:iff=0andc=6thengosub3300
  175. 4200 return
  176. 4300 onc-5gosub5660,5750:ifp=0orc+u=16thenreturn
  177. 4310 ifc=7thengosub5500:iff=0thenprint#2,"2newea":print#2,p:Dg:iff>0thenreturn
  178. 4330 ifc=6theninput#2,p:ifp>uthenprint"nur"u"von"p:p=u
  179. 4340 fori=1top:q=p(i)
  180. 4345 ifc=6theninput#2,k$(q):printq;k$(q):Dg:print#2,chr$(34)+k$(q)+chr$(34)
  181. 4350 j=0:D :ifc=6theninput#2,f(q,j),u(q,j),w(q,j)
  182. 4360 ifc=7thenprint#2,f(q,j):print#2,u(q,j):print#2,w(q,j)
  183. 4370 j=j+1:D)j>f(q,0):f(q,j)=0:nexti:gosub5550:return
  184. 4400 gosub6200:s=t0*f(q,p):ifs=0orn=0thengosub6400:return
  185. 4410 fori=1ton:u=k(i):forj=1ton:v=k(j):s(u,v)=y(1,i,j)*s-y(2,i,j)/s
  186. 4420 r(u,v)=y(0,i,j):nextj:s(i,m)=0:r(i,m)=0:nexti:r(k(1),m)=1/b(0)
  187. 4450 k=j2:D::D5 ifk=j2+2*n(6):i=k(i(k+1)):ifi=0then4490
  188. 4460 r=1+s*s/a(k)^2:u=b(k)/r:v=u*s/a(k)
  189. 4470 j=k(i(k)):r(i,j)=r(i,j)-u:s(i,j)=s(i,j)+v
  190. 4480 j=k(j(k)):r(i,j)=r(i,j)+u:s(i,j)=s(i,j)-v
  191. 4490 k=k+2:D6
  192. 4500 fork=1ton-1:u=r(k,k):v=s(k,k):f=u*u+v*v
  193. 4510 iff=0thengosub6400:return
  194. 4550 fori=k+1ton:ifnotm%(i,k)then4600
  195. 4560 x=r(i,k):y=s(i,k):r=(x*u+y*v)/f:s=(y*u-x*v)/f:forj=k+1tom
  196. 4580 ifm%(i,j)then:x=r(k,j):y=s(k,j):r(i,j)=r(i,j)-r*x+s*y:s(i,j)=s(i,j)-s*x-r*y
  197. 4590 nextj
  198. 4600 nexti,k:u=r(n,n):v=s(n,n):f=u*u+v*v:iff=0thengosub6400:return
  199. 4610 x=r(n,m):y=s(n,m):r=(x*x+y*y)/f
  200. 4620 gosub6250:ifr=0thenprint"ua=0":p=p-1:l=1:return
  201. 4630 u(q,p)=tl*log(r):ifl=2thenreturn
  202. 4640 r=x*u+y*v:s=y*u-x*v
  203. 4650 ifs=0thenw(q,p)=(1-sgn(r))*90:Dg:w(q,p)=90*sgn(s)-tp*atn(r/s)
  204. 4700 x=f(q,p):printx;:printspc(12-len(str$(x)));:x=u(q,p):gosub5950
  205. 4710 x=w(q,p):printspc(7):gosub5950
  206. 4720 print:ifl=0orp<2thenreturn
  207. 4750 f=f(q,p):u=u(q,p):v=w(q,p):forj=pto2step-1
  208. 4760 iff(q,j-1)>fthenf(q,j)=f(q,j-1):u(q,j)=u(q,j-1):w(q,j)=w(q,j-1):nextj
  209. 4770 f(q,j)=f:u(q,j)=u:w(q,j)=v:return
  210. 5000 print"korrektur? ";
  211. 5010 D7,20:print"[158] [155]";:gosub5050:print"[157] [157]"c$:DI:return
  212. 5050 poke198,0:wait198,1:getc$:return
  213. 5100 print" "e$(i);j,i(k);j(k);:inputi(k),j(k):return
  214. 5150 print" "left$(d$(i),i-3);j,e$(i);i(k);j(k);i(k+1);
  215. 5160 inputi(k),j(k),i(k+1):j(k+1)=0:return
  216. 5200 D :print" vp"j,e$(i):printi(k);j(k);i(k+1);j(k+1);
  217. 5210 inputi(k),j(k),i(k+1),j(k+1):D)j(k+1)<=n:return
  218. 5250 forq=ktok+i/4
  219. 5260 ifi(q)<0orj(q)<0orn<i(q)orn<j(q)thenj=j-1:return
  220. 5270 nextq:k=q:return
  221. 5300 D :gosub6500:inputb(k):D)b(k)<>0:return
  222. 5350 D :gosub6550:inputa(k),b(k),b(k+1):D)b(k+1)>0anda(k)>0:return
  223. 5400 D :gosub6600:inputb(k),b(k+1):D)b(k+1)>0:return
  224. 5450 gosub6650:forl=0to3:gosub6660
  225. 5460 inputx(0,l,j),x(1,l,j),x(2,l,j):nextl:return
  226. 5500 D :print:input"filename ($=directory)";n$:ifn$="$"thenprint:Dj"$:*=s"
  227. 5510 print:gosub5550:iff>0thenreturn
  228. 5520 D)n$<>"$":open2,8,2,n$+",s,"+e$(c+2)
  229. 5550 DF1:input#15,f,c$,r,s:iff>0thenprintf;c$;r;s
  230. 5560 return
  231. 5600 j1=n(1)+n(2)+n(3)+2*(n(4)+n(5)+n(6)):return
  232. 5620 q=k+p-m:i(k)=i(q):j(k)=j(q):b(k)=b(q):a(k)=a(q):return
  233. 5650 print"":printe$(7)"belegung:"
  234. 5660 u=0:fori=1to9:k=f(i,0):ifk=0thenu=u+1:p(u)=i
  235. 5670 ifk>0andc<>6thenprinti":"k$(i):print" "k"freq :"f(i,1)"-"f(i,k)
  236. 5680 nexti:print:p=u:return
  237. 5700 D :D :printe$(7)" (1-9)"q;:inputq
  238. 5710 D)q>=0andq<=9:D)f(q,0)>0orc<1orq=0:return
  239. 5750 gosub5650:D :fori=1to9-u:printi".";:q=p(i):gosub5700:p(i)=q
  240. 5760 ifq>0thennexti
  241. 5770 p=i-1:gosub5000:D)c$<>"j":return
  242. 5800 p=i:o=j:ifl=1thenr=b(k):Dg:r=1/b(k)
  243. 5820 ifi=jorp=othenreturn
  244. 5830 y(l,p,i)=y(l,p,i)+r:y(l,o,j)=y(l,o,j)+r:m%(p,i)=u:m%(o,j)=u
  245. 5840 y(l,p,j)=y(l,p,j)-r:y(l,o,i)=y(l,o,i)-r:m%(p,j)=u:m%(o,i)=u:return
  246. 5850 forl=0to2:r=x(l,kf,s):ifr<>0thengosub5820
  247. 5860 nextl:return
  248. 5900 print:printd$(7);spc(5)"betrag (db)   ";
  249. 5910 ifw(q,0)=0thenprint"phase (grad)":Dg:print"toleranz(db)"
  250. 5920 return
  251. 5950 x$=str$(sgn(x)*int(abs(x)*1e3+.5))
  252. 5960 D::D5 iflen(x$)>4:x$=D[("0",x$,1):D6
  253. 5970 x$=D[(".",x$,len(x$)-3):D,"####.###",x$:return
  254. 6000 printDu(" ",10)"[145]":printf(q,p)tab(8):x=u(q,p):gosub5950:inputu(q,p)
  255. 6010 D :print"[145]"tab(23):x=w(q,p):gosub5950
  256. 6020 inputw(q,p):D)w(q,p)>0:gosub4750:return
  257. 6050 r=1:c$=b$(p):l=len(c$):ifc$="0"thenp=10-p:return
  258. 6060 fori=1to3:forj=rtor+2:ifmid$(c$,j,1)<>e$(i)thennextj,i:return
  259. 6070 k=val(right$(c$,l-j)):ifk<int(i/2)ork>n(i)thenreturn
  260. 6080 forj=1toi:k=k+n(j):nextj:k=k-n(i):p(p+4)=k:ifr=1then6120
  261. 6100 fori=3to4:n$=mid$(c$,i,1)
  262. 6110 ifn$="*"thenp(p+4)=k+50:Dg:ifn$<>"/"thennexti:return
  263. 6120 ifr=1thenp(p)=k:r=4:ifl>4then6060:Dg:p=p+1:return
  264. 6200 print"[150]***  busy  ***[155]":return
  265. 6250 print"[145]"Du(" ",14)"[145]":return
  266. 6300 fori=1ton:r=m%(i,p):m%(i,p)=m%(i,q):m%(i,q)=r:nexti
  267. 6310 forj=0tom:r=m%(p,j):m%(p,j)=m%(q,j):m%(q,j)=r:nextj:return
  268. 6350 r=u:forj=1ton:r=r-m%(k,j):nextj:return
  269. 6400 gosub6250:r=0:print"nicht loesbar":p=p-1:l=1:return
  270. 6500 print" "e$(i);j"("i(k);j(k)")",b(k);:return
  271. 6550 print" t"j"("i(k);j(k);i(k+1)") b0,ic,ft":printa(k);b(k);b(k+1);:return
  272. 6600 print" op"j"("i(k);j(k);i(k+1)") v0,ft":printb(k);b(k+1);:return
  273. 6650 print" vp"j"("i(k);j(k);i(k+1);j(k+1)") 1/r  c  1/l":return
  274. 6660 print"  y"11+8*int(l/2)+l;x(0,l,j);x(1,l,j);x(2,l,j);:return
  275. 6700 gosub6650:forl=0to3:gosub6660:print:nextl:return
  276. 7000 input"notiz";n$:gosub7400:print"   "n$
  277. 7020 k=0:fori=1to6:ifn(i)>0ori=1thenprint:printd$(i)" "e$(i)
  278. 7030 Df:forj=sgn(i-1)ton(i):onigosub6500,6500,6500,6550,6700,6600
  279. 7040 Df:print:k=k+1+int(i/4):nextj
  280. 7050 nexti:print:print"ausgangsknoten:",n(9):print#4:close4:return
  281. 7100 gosub7400:print"optimierungsverlauf":print#4:close4:return
  282. 7200 gosub7400:print"   "k$(q):gosub1400:print#4:close4:return
  283. 7300 gosub7400:fori=1toq:ifq>1thenprint"   "p(i);
  284. 7310 print"   "k$(p(i)):nexti:print:print#4:close4:DW:return
  285. 7400 open4,4:cmd4:print:printchr$(14)"ac-netzwerkanalyse"chr$(15)
  286. 7410 print:return
  287.