home *** CD-ROM | disk | FTP | other *** search
/ 8bitfiles.net/archives / archives.tar / archives / genie-commodore-file-library / C64Educational / ROOT-GUESSER (.txt) < prev    next >
Commodore BASIC  |  2019-04-13  |  1KB  |  56 lines

  1. 1 PRINT"":PRINT"[147]":POKE 53280,0:POKE 53281,0
  2. 10 PRINT" THIS PROGRAM HELPS DETERMINE THE ROOTS":PRINT
  3. 15 PRINT"OF AN (N) DEGREE POLYNOMIAL THAT":PRINT
  4. 20 PRINT"EQUALS ZERO.":PRINT
  5. 21 PRINT
  6. 22 PRINT"A(1)+A(2)(X)+...+A(N+1)(X**N)=0":PRINT
  7. 25 PRINT"ENTER DEGREE OF POLYNOMIAL:":INPUT N:IF N=0 GOTO 200
  8. 30 PRINT:PRINT"ENTER COEFFICIENTS:":PRINT""
  9. 32 DIM A(N+1)
  10. 35 FOR I=1 TO N+1
  11. 40 PRINT" A("I")=:":INPUT A(I)
  12. 45 NEXT I
  13. 46 PRINT:PRINT"DO YOU WANT TO PRINT OUTPUT(Y/N)?"
  14. 47 GET AA$:IF AA$="" THEN 47
  15. 48 IF AA$="Y" THEN CT=1.0:OPEN4,4:A=1
  16. 49 PRINT"[147]":GOSUB 500
  17. 50 PRINT:PRINT"THE LAST GUESS WAS X="X",F(X)="Y
  18. 51 PRINT:PRINT"ENTER X(THE POSSIBLE ROOT):":INPUT X:Y=0
  19. 65 Y=A(1)
  20. 70 FOR J=1 TO N
  21. 75 Y=Y+(A(J+1)*(X^J))
  22. 80 NEXT J
  23. 85 PRINT"X=" X
  24. 90 PRINT"":PRINT"F(X)=" Y:PRINT:PRINT"FOR X TO BE A ROOT Y=0."
  25. 91 IF CT=1.0 THEN PRINT#4,"X=" X:PRINT#4,"F(X)="Y:PRINT#4
  26. 95 PRINT:PRINT"DO YOU WISH TO TRY ANOTHER VALUE?(Y/N)
  27. 97 [161] X$:[139] X$[178]"" [167] 97
  28. 101 [139] X$[178]"N" [167] [137] 109
  29. 103 [153]"LOAD":CT[178]0.0:[141]500
  30. 104 [137] 50
  31. 109 [153]"LOAD":[153]"DO YOU WISH TO:"
  32. 110 [153]:[153]"      1)RE-RUN PROGRAM."
  33. 111 [153]:[153]"      2)EXIT TO BASIC."
  34. 112 [153]:[153]:[153]"PICK A NUMBER:"
  35. 113 [133] BG$
  36. 114 [139] BG$[178]"1"[167] [138]
  37. 115 [139] BG$[178]"2"[167] [137]117
  38. 116 [137] 109
  39. 117 [160]4:[162]:[128]
  40. 200 [153]:[153]" NOT A VALID INPUT":[153]:[137]25
  41. 500 [153]"THE POLYNOMIAL IS OF DEGREE"N:[153]
  42. 501 [139] CT[178]1.0 [167] [152]4,"THE POLYNOMIAL IS OF DEGREE"N:[152]4,
  43. 505 [153]"THE POL    
  44. 16713 L IS:"
  45. 506 [139] CT[178]1.0 [167] [152]4,"THE POLYNOMIAL IS:"
  46. 507 [139] CT[178]1.0 [167] [152]4," F(X)="
  47. 508 [153]" F(X)="
  48. 509 [129] J[178]1 [164] N[170]1
  49. 510 [153] A(J)"*(X**"(J[171]1)")"
  50. 511 [139] CT[178]1.0 [167] [152]4,A(J)"*(X**"(J[171]1)")"
  51. 515 [130] J
  52. 516 [153]" =0.0"
  53. 517 [139] CT[178]1.0 [167] [152]4," =0.0":[152]4
  54. 520 [139] A[178]1 [167] CT[178]1
  55. 600 [142]
  56.