home *** CD-ROM | disk | FTP | other *** search
/ Play and Learn 2 / 19941.ZIP / 19941 / EDUCTORY / GEOGRAPY / M15.BAS (.txt) < prev    next >
Encoding:
GW-BASIC  |  1994-02-04  |  1.0 KB  |  36 lines
  1. 0  CLS:PRINT
  2. 10   CLS:PRINT "INTEGRATION: GAUSSIAN QUADRATURE"
  3. 20   PRINT 
  4. 30  DEF FNC(X)=X^3:REM DEFINE FUNCTION HERE (DEFFNC(X)="FUNCTION")
  5. 40   DATA  .076526521,.15275339,.22778585,.14917299,.37370609
  6. 50   DATA  .14209611,.510867,.13168864,.6305368,.11819453
  7. 60   DATA  .74633191,.10193012,.83911697,.083276742,.91223443
  8. 70   DATA  .062672048,.96397183,.04060143,.9931286,.017614007
  9. 80   PRINT "INTEGRATION LIMITS (LOWER,UPPER)";
  10. 90   INPUT X,Y
  11. 100   PRINT "NUMBER OF INTERVALS";
  12. 110   INPUT N
  13. 120  S = (Y - X) / N / 2
  14. 130  T = X + S
  15. 140  R = 0
  16. 150   FOR I = 1 TO N
  17. 160  P = 0
  18. 170   FOR J = 1 TO 10
  19. 180   READ A,B
  20. 190  P = P + B * ( FNC(S * A + T) +  FNC(T - S * A))
  21. 200   NEXT J
  22. 210   RESTORE 
  23. 220  R = R + P * S
  24. 230  T = T + 2 * S
  25. 240   NEXT I
  26. 250   PRINT "INTEGRAL = ";R
  27. 260   PRINT 
  28. 270   PRINT "CHANGE DATA AND RECOMPUTE?"
  29. 280   PRINT "(0=NO,1=NEW INTEGRATION LIMITS, 2=NEW NO. OF INTERVALS)";
  30. 290   INPUT S
  31. 300   IF S = 1 THEN 80
  32. 310   IF S = 2 THEN 100
  33. 315  INPUT "Run again (Y/N) ";A$
  34. 317  IF A$="Y" OR A$="y" THEN 10
  35. 320   RUN "MATH"
  36.