home *** CD-ROM | disk | FTP | other *** search
/ Antennas / Antennas_CD-ROM_Walnut_Creek_September_1996.iso / mathaids / misc / simpsoni.asc < prev   
Text File  |  1996-06-30  |  1KB  |  53 lines
  1. 1 N$="INTEG1"
  2. 2 REM**********************************
  3. 3 REM**   SIMPSON'S RULE OF         ***
  4. 4 REM**   NUMERICAL INTEGRATION (C) ***
  5. 5 REM**   K.VANSLETTE - JUNE 1984   ***
  6. 6 REM**   32 PGMS FOR C-64, P259    ***
  7. 7 REM**********************************
  8. 8 S$=" *** "
  9. 9 PRINT "";TAB(16-LEN(N$)/2);S$;N$;S$
  10. 10 GOTO 90:REM**************************
  11. 20 REM
  12. 30 REM ENTER EQUATION AT LINE 50
  13. 40 REM
  14. 50 Y=4/(1+X^2):REM PI
  15. 60 REM
  16. 70 REM
  17. 80 RETURN:REM**************************
  18. 90 PRINT "  THIS PROGRAM WILL INTEGRATE EQUATIONS USING SIMPSON'S RULE."
  19. 100 PRINT "  ENTER THE EQUATION AT LINE 50 AS":PRINT"  Y = F(X) - - - Y/N"
  20. 110 INPUT "                     ";X$
  21. 120 IF X$="Y" THEN LIST 30-50
  22. 130 PRINT "======================================="
  23. 140 INPUT "ENTER LOWER LIMIT OF X ----";LX
  24. 150 INPUT "ENTER UPPER LIMIT OF X ----";UX
  25. 160 PRINT "======================================="
  26. 170 PRINT " HIT ANY KEY TO HALT..."
  27. 180 N=2: PRINT
  28. 190 PRINT "=SEGMENTS=   =INTEGRAL="
  29. 200 GET R$:IF R$=""THEN 200
  30. 210 DX=(UX-LX)/N:T=0
  31. 220 X=LX:GOSUB50:T=T+Y
  32. 230 X=UX:GOSUB50:T=T+Y
  33. 240 M=N/2:Z=0
  34. 250 FOR J=1 TO M
  35. 260 X=LX+DX*(2*J-1)
  36. 270 GOSUB 50:Z=Z+Y:NEXT J
  37. 275 GOSUB 50:K0=D
  38. 280 T=T+4*Z:M=M-1:IF M=0 THEN 320
  39. 290 Z=0:FOR J=1 TO M
  40. 300 X=LX+DX*2*J
  41. 310 GOSUB 50:Z=Z+Y:NEXT J:T=T+2*Z
  42. 320 A=DX*T/3
  43. 330 PRINT "  ";N,"  ";A
  44. 340 N=N*2:GOTO 210
  45. 400 END
  46. 1000 GET R$
  47. 1010 IF R$="" THEN RETURN
  48. 1020 GET R$
  49. 1030 IF R$="" THEN 1020
  50. 1040 RETURN
  51. 1050 END
  52.  
  53.