home *** CD-ROM | disk | FTP | other *** search
/ 64'er 1987 January / 64er_Magazin_87-01_1987_Markt__Technik_de.d64 / lineare gl. (.txt) < prev    next >
Encoding:
Commodore BASIC  |  1987-01-01  |  1.1 KB  |  36 lines
  1. 1 rem loesen eines lineraren gleichungsystems mit hilfe des horner-schemas
  2. 2 :
  3. 3 rem germano caronni greubstelstr.10 ch-5430 wettingen tel.0041/56262574
  4. 4 :
  5. 10 poke53280,0:poke53281,0:print"[147]"
  6. 20 input"anzahl unbekannte ";n:ifn<2orn>10thenprint"[145][145]":goto20
  7. 30 dima(n,n),b(n)
  8. 40 fori=1ton:print"zeile";i;:forj=1ton:printtab(12);"spalte ";j;
  9. 50 inputa(i,j):nextj
  10. 60 printtab(12);"konstante ";
  11. 70 inputa(i,0):nexti
  12. 80 print
  13. 90 :
  14. 100 fori=1ton-1:fl=0:dk=i
  15. 102 ifa(i,i)<>0then105
  16. 103 dk=dk+1:ifdk>nthenprint"tut mir leid, ueberbestimmt !":end
  17. 104 forj=0ton:a=a(i,j):a(i,j)=a(dk,j):a(dk,j)=a:next:print"*-":goto102
  18. 105 a=a(i,i)*-1:fork=iton:a(i,k)=a(i,k)/a:nextk:a(i,0)=a(i,0)/a
  19. 110 forj=i+1ton
  20. 119 ifa(j,i)=0thenprint"*+":goto140
  21. 120 a=a(j,i):fork=iton:a(j,k)=a(j,k)/a:nextk:a(j,0)=a(j,0)/a
  22. 130 fork=iton:a(j,k)=a(j,k)+a(i,k):nextk:a(j,0)=a(j,0)+a(i,0)
  23. 140 nextj,i
  24. 145 :
  25. 150 fori=nto1step-1
  26. 155 ifa(i,i)=0thenprint"die gleichung hat keine eindeutige loesung":end
  27. 160 a(i,0)=a(i,0)/a(i,i):a(i,i)=1:b(i)=a(i,0):printi,b(i):ifi=1then180
  28. 170 forj=1toi-1:a(j,0)=a(j,0)-a(j,i)*b(i):a(j,i)=0:nextj
  29. 180 nexti
  30. 190 print
  31. 199 :
  32. 500 fori=1ton:forj=1ton:printa(i,j);"[157]";:nextj:print"=";a(i,0):nexti
  33. 510 print""
  34. 520 input"name ";a$:ifa$=""thenend
  35. 530 open2,8,2,"horn"+str$(n)+"loesung,s,w":fori=1ton:print#2,b(i):next:close2
  36.