home *** CD-ROM | disk | FTP | other *** search
/ 64'er 1991 November / 64er_Magazin_91-11_1991_Markt__Technik_de_Side_A.d64 / gauss (.txt) < prev    next >
Commodore BASIC  |  2022-10-26  |  712b  |  21 lines
  1. 1 print"[147][155]":poke53280,0:poke53281,0:poke19,64
  2. 2 input"anzahl der unbekannten/gleichungen: ";n:dimgl(n+1,n):poke19,0
  3. 3 print"":gosub15:fora=1ton:ifgl(a,a)=0thengosub7:ifs=0then5
  4. 4 gosub9:gosub13:nexta
  5. 5 gosub13:print"loesung:":gosub6:print:end
  6. 6 fori=1ton:forj=1ton:printgl(j,i);:nextj:print"[194]"gl(j,i):nexti:return
  7. 7 fori=aton:ifgl(a,i)=0thennexti:ifi>nthens=0:return
  8. 8 s=1:forj=1ton+1:h=gl(j,a):gl(j,a)=gl(j,i):gl(j,i)=h:nextj:return
  9. 9 fori=1ton:ifi=athen12
  10. 10 f1=gl(a,a):f2=gl(a,i):forj=1ton+1
  11. 11 gl(j,i)=f2*gl(j,a)-f1*gl(j,i):nextj
  12. 12 nexti:return
  13. 13 fori=1ton:kf=gl(i,i):forj=1ton+1:ifkf<>0thengl(j,i)=gl(j,i)/kf
  14. 14 nextj:nexti:return
  15. 15 fori=1ton:printi"[157].gl.: ";:forj=1ton+1:ifj=n+1thenprint"[194] ";
  16. 16 poke19,64:inputh$:poke19,0:ifpos(0)=40andj<>n+1thenprint
  17. 17 forl=1tolen(h$):ifmid$(h$,l,1)<>"/"thennextl
  18. 18 h=val(h$):ifl<len(h$)thenh=val(left$(h$,l-1))/val(right$(h$,len(h$)-l))
  19. 19 gl(j,i)=h:ifj<>n+1thenprint" ";
  20. 20 nextj:print:nexti:return
  21.