home *** CD-ROM | disk | FTP | other *** search
/ HAM Radio 1 / HamRadio.cdr / tech / rf3 / bessel.bas (.txt) < prev    next >
Encoding:
GW-BASIC  |  1979-12-31  |  639 b   |  25 lines
  1. 10  PRINT, "THIS IS A PROGRAM TO CALCULATE THE BESSEL FUNCTION "
  2. 15  PRINT:PRINT,"K2UYH -------   MODIFIED IBM PC  ---- WA2TIF"
  3. 20  INPUT "ORDER";N
  4. 30  INPUT "VALUE";X
  5. 40  REM DEF FNFAC(I)
  6. 50  REM IF I=0 THEN FNRETURN 1
  7. 60  REM FNEND FNFAC(I-1)*I
  8. 65  GOSUB 150
  9. 70  J0=X^N/(2^N*FNFAC(N))
  10. 80  J1=(1-(X^2/(2^2*(N+1))))
  11. 85  GOSUB 150
  12. 90  J2=X^4/(2^4*FNFAC(2)*(N+1)*(N+2))
  13. 95  GOSUB 150
  14. 100  J3=X^6/(2^6*FNFAC(3)*(N+1)*(N+2)*(N+3))
  15. 105  GOSUB 150
  16. 110  J4=X^8/(2^8*FNFAC(4)*(N+1)*(N+2)*(N+3)*(N+4))
  17. 115  GOSUB 150
  18. 120  J5=X^10/(2^10*FNFAC(5)*(N+1)*(N+2)*(N+3)*(N+4)*(N+5))
  19. 130  J=J0*(J1+J2-J3+J4-J5)
  20. 140  PRINT "J";N;"(";X;")=";J
  21. 145  END
  22. 150  IF I=0 GOTO 160 ELSE GOTO 170
  23. 160  DEF FNFAC(I)=1:RETURN
  24. 170  DEF FNFAC(I)=(I-1)*I:RETURN
  25.