[Prev][Next][Index][Thread]
Medhurts K formula [Qbasic listing]
----------
From: Edward V. Phillips [SMTP:ed-at-alumni.caltech.edu]
Sent: Thursday, April 09, 1998 2:25 PM
To: tesla-at-stic-dot-net
Subject: Re: Medhurts K formula
Rick:
Here is the file for calculating winding self and mutual
inductance. It runs in QBASIC on my Mac and, with perhaps some
tuning, it should run on a PC; I haven't tried that. It is a
text file. Please let me know how it comes across.
Ed
' ********** THIS SEEMS TO BE THE RIGHT ONE ON 6/24/90 **********
' THIS IS GOOD MUTUAL, SAVED 3:45 PM, 6/24/90
' IT HAS BEEN EDITED FOR SHORT LINES ON 9/4/96, AND HOPEFULLY CAN BE SENT
' BY EMAIL. REMAINS TO BE TESTED.
10 ' ********** MUTUAL **********
12 '
14 ' THIS IS A MUTUAL INDUCTANCE PROGRAM WHICH SEEMS TO WORK.
16 '
18 ' IT IS BASED ON MATERIAL IN THE BOOK:
20 '
22 ' "INDUCTANCE CALCULATIONS
24 ' WORKING FORMULAS AND TABLES"
26 ' 'BY FREDERIC W. GROVER
28 '
30 ' DOVER PUBLICATION S974
32 '
34 ' MUTUAL INDUCTANCE IS CALCULATED BY FORMULA (108),
36 ' PAGE 123, OF CHAPTER 15, USING DOUBLE DIFFERENCE
38 ' INTERPOLATION OF TABLES 29 AND 30.
40 '
42 ' INDUCTANCE IS CALCULATED BY FORMULAS (118) & (120)
44 ' OF CHAPTER 16, USING DOUBLE DIFFERENCE INTERPOLATION
46 ' OF TABLES 36 AND 37. TABLES 38 AND 39 ARE APPROXIMATED
48 ' BY FORMULAS FROM "THE RADIOTRON DESIGNERS HANDBOOK".
50 '
52 ' WIRE TABLES ARE INCLUDED AND GIVE DIAMETER AND
54 ' BOTH RESISTANCE AND WEIGHT PER 1000 FEET FOR
56 ' B&S SIZES FROM 10 TO 40.
58 ' '
60 ' CORRECTED (1600,1610) TO USE CORRECT TABLE ENTRIES
70 ' AND SAVED 1756 Z, 02/07/86
100 CLS
110 OV$="OLD ="
150 DIM W(3,40),BN(21,21),LN(11,11),KN(2,101),B(21)
160 GOSUB 6000
180 GOTO 400
200 ON ERROR GOTO 400
400 CLS:PRINT
405 OV=D1:LOCATE 2,40:PRINT OV$;OV
410 LOCATE 2,1:INPUT "DIAMETER OF SMALLER COIL (INCHES)";D1:LOCATE 2,40:PRINT OV
$;OV
415 IF D1=0 THEN D1=OV
420 PRINT
430 OV=LS:LOCATE 4,40:PRINT OV$;OV
440 LOCATE 4,1:INPUT " LENGTH OF SMALLER COIL (INCHES)";LS:LOCATE 4,40:PRINT OV$
;OV
450 IF LS=0 THEN LS=OV
460 PRINT
470 OV1=NS:OV2=WS:LOCATE 6,40:PRINT OV$;NS","WS
480 LOCATE 6,1:INPUT " TURNS, WIRE SIZE OF SMALLER COIL";NS:LOCATE 6,35:INPUT WS
:LOCATE 6,40:PRINT OV$;OV1;","OV2
490 IF NS=0 THEN NS=OV1:WS=OV2
500 PRINT
510 OV=D2:LOCATE 8,40:PRINT OV$;D2
520 LOCATE 8,1:INPUT " DIAMETER OF LARGER COIL (INCHES)";D2:LOCATE 8,40:PRINT OV
$;OV
530 IF D2=0 THEN D2=OV
540 PRINT
550 OV=LL:LOCATE 10,40:PRINT OV$;LL
560 LOCATE 10,1:INPUT " LENGTH OF LARGER COIL (INCHES)";LL:LOCATE 10,40:PRINT OV
$;OV
570 IF LL=0 THEN LL=OV
580 PRINT
590 OV1=NL:OV2=WL:LOCATE 12,40:PRINT OV$;NL","WL
600 LOCATE 12,1:INPUT " TURNS, WIRE SIZE OF LARGER COIL";NL:LOCATE 12,35:INPUT W
L:LOCATE 12,40:PRINT OV$;OV1;","OV2
610 IF NL=0 THEN NL=OV1:WL=OV2
620 PRINT
630 OV=SE:LOCATE 14,40:PRINT OV$;SE
640 LOCATE 14,1:INPUT "DISTANCE BETWEEN CENTERS OF COILS";SE:LOCATE 14,40:PRINT
OV$;OV
650 IF SE=0 THEN SE=OV
660 PRINT
680 PRINT "CHANGE ANY OF THESE QUANTITIES?" '
700 T$=INKEY$:IF T$="" THEN 700
720 IF T$ <>"Y" AND T$<>"N" THEN 700
740 IF T$="N" THEN 980
760 GOTO 400
980 PRINT:PRINT " CALCULATING "
;
1000 DS=D1+W(1,WS) ' AVERAGE DIAMETER OF SMALLER COIL
1020 KS=NS*W(1,WS)/(LS-W(1,WS)) 'WINDING FACTOR OF SMALLER COIL
1030 IF KS<.999 THEN 1040
1035 LOCATE 15,1:PRINT "NO MORE THAN";LS/W(1,WS)"; TURNS ON SMALLER COIL!":GOTO
1100
1040 DB=D2+W(1,WL) ' AVERAGE DIAMETER OF LARGER COIL
1060 IF D2>=DS THEN 1160
1080 LOCATE 15,1:PRINT "'SMALL' COIL HAS BIGGER DIAMETER THAN 'BIG' ONE! !"
1085 PRINT "MINIMUM DIAMETER OF OUTER COIL =";DS+W(1,WL)"INCHES"
1100 PRINT " TRY AGAIN "
1110 T$=INKEY$:IF T$="" THEN 1110
1120 CLS
1140 GOTO 760
1160 KL=NL*W(1,WL)/(LL-W(1,WL)) 'WINDING FACTOR OF LARGER COIL
1170 IF KL<.999 THEN 1180
1175 LOCATE 15,1:PRINT "NO MORE THAN";LL/W(1,LL);"TURNS ON LARGER COIL!":GOTO 1
100
1180 RS=1.27*DS ' RADIUS OF SMALLER COIL (CM)
1200 M1=1.27*LS ' HALF LENGTH OF SMALLER COIL (CM)
1220 N1=.5*NS/M1 ' WINDING PITCH OF SMALLER COIL
1240 AL=1.27*DB ' RADIUS OF LARGER COIL (CM)
1260 M2=1.27*LL ' HALF LENGTH OF LARGER COIL (CM)
1280 N2=.5*NL/M2 ' WINDING PITCH OF LARSER COIL
1300 S=2.54*SE ' CENTER TO CENTER COIL SPACING (CM)
1320 X(1)=S+(M1+M2)
1340 X(2)=S+(M1-M2)
1360 X(3)=S-(M1-M2)
1380 X(4)=S-(M1+M2)
1400 A2=AL^2
1420 FOR U=1 TO 4
1440 R(U)=SQR(A2+X(U)^2)
1460 O(U)=A2/(A2+X(U)^2)
1480 NEXT U
1500 AF=DS/DB ' "ALPHA"
1520 FOR U=1 TO 4
1540 IF AF>=.9 AND O(U)>=.9 THEN 1600
1560 IF AF<=.05 OR O(U)<=.05 THEN 1700
1580 GOTO 1800
1600 X=100*(AF-.9)
1620 Y=100*(O(U)-.9)
1640 GOSUB 4000
1660 GOTO 1860
1700 X=20*AF
1720 Y=20*O(U)
1740 GOSUB 4240
1760 GOTO 1860
1800 X=20*AF
1820 Y=20*O(U)
1840 GOSUB 4360
1860 NEXT U
1900 BR=R(1)*B(1)-R(2)*B(2)-R(3)*B(3)+R(4)*B(4)
1920 M=.019739208#*RS^2*N1*N2*BR ' CALCULATE MUTUAL AT LAST! !
1940 GOSUB 5000
3000 CLS
3005 PRINT
3010 PRINT TAB(15)"DIAMETER"TAB(28)"LENGTH"TAB(41)"TURNS"TAB(54)"B&S #"
3020 PRINT
3030 PRINT "SMALL COIL"TAB(15)D1 TAB(28)LS TAB(41)NS TAB(54)WS
3040 PRINT "LARGE COIL"TAB(15)D2 TAB(28)LL TAB(41)NL TAB(54)WL
3050 PRINT
3060 PRINT "CENTER TO CENTER SPACING =";SE
3070 PRINT
3080 PRINT "INDUCTANCE OF SMALL COIL ="SI" MICROHENRIES"
3090 PRINT
3100 PRINT "INDUCTANCE OF LARGE COIL ="IL" MICROHENRIES"
3110 PRINT
3120 PRINT "MUTUAL INDUCTANCE ="M" MICROHENRIES"
3130 PRINT
3140 PRINT "COUPLING COEFFICIENT ="M/SQR(SI*IL)
3150 PRINT
3190 TT$="N"
3200 PRINT "PRINT RESULTS? ";
3202 TT$=INKEY$:IF TT$<>"Y" AND TT$<>"N" THEN 3202
3210 IF TT$="Y" THEN 3300
3220 CLS
3230 PRINT "ANOTHER RUN?"
3240 TT$=INKEY$:IF TT$<>"Y" AND TT$<>"N" THEN 3240
3250 IF TT$="Y" THEN 400
3260 GOTO 10000
3300 CLS
3310 TS$="N"
3320 INPUT "PRINT INTERNAL VALUES";ZZ
3500 LPRINT CHR$(27)"C"CHR$(33);
3510 LPRINT
3520 LPRINT TAB(16)"DIAMETER"TAB(30)"LENGTH"TAB(44)"TURNS"TAB(58)"B&S #"
3530 LPRINT
3540 LPRINT "SMALL COIL"TAB(16)D1TAB(30)LSTAB(44)NSTAB(58)WS
3550 LPRINT "LARGE COIL"TAB(16)D2TAB(30)LLTAB(44)NLTAB(58)WL
3560 LPRINT
3570 LPRINT "CENTER TO CENTER SPACING =";SE
3580 LPRINT
3590 LPRINT "INDUCTANCE DF SMALL COIL ="SI" MICROHENRIES"
3600 LPRINT
3610 LPRINT "INDUCTANCE OF LARGE COIL ="IL" MICROHENRIES"
3620 LPRINT
3630 LPRINT "MUTUAL INDUCTANCE ="M" MICROHENRIES"
3640 LPRINT
3650 LPRINT "COUPLING COEFFICIENT ="M/SQR(SI*IL)
3660 LPRINT
3670 LPRINT
3680 IF TS$="N" THEN 3770
3690 LPRINT
3700 LPRINT "X1 ="X(1),"X2 ="X(2),"X3 ="X(3)," X4 ="X(4)
3710 LPRINT
3720 LPRINT "ALPHA ="AF
3730 LPRINT
3740 FOR U=1 TO 4
3750 LPRINT "RO(";U;")="O(U)" ","B(";U;")=";B(U)
3760 NEXT U
3770 LPRINT CHR$(140);
3780 CLS
3790 GOTO 400
4000 XI=INT(X)
4002 GOTO 4020
4010 CLS: PRINT 4010; U,X:STOP
4020 YI=INT(Y)
4040 YF=Y-YI
4060 XF=X-XI
4080 YA=YF*(YF-1)*.5
4100 YB=1-YF*YF
4120 YC=YF*(YF+1)*.5 '
4140 Y1=LN(XI-1,YI-1)*YA+LN(XI-1,YI)*YB+LN(XI-1,YI+1)*YC
4160 Y2=LN(XI,YI-1)*YA+LN(XI,YI)*YB+LN(XI,YI+1)*YC
4180 Y3=LN(XI+1,YI-1)*YA+LN(XI+1,YI)*YB+LN(X,YI+1)*YC
4200 B(U)=Y1*XF*(XF-1)*.5+Y2*(1-XF*XF)+Y3*XF*(XF+1)*.5
4220 RETURN
4240 XI=INT(X)
4242 GOTO 4260
4250 CLS:PRINT 4230,U,X:STOP
4260 YI=INT(Y)
4280 XF=(X-XI)
4300 YF=(Y-YI)
4320 B(U)=BN(XI,YI)*(1-XF)*(1-YF)+BN(XI,YI+1)*YF*(1-XF)+BN(XI+1,YI)*XF*(1-YF)+BN
(XI+1,YI+1)*XF*YF
4340 RETURN
4360 XI=INT(X)
4362 GOTO 4380
4370 CLS:PRINT 4350,U,X:STOP
4380 YI=INT(Y)
4400 YF=Y-YI
4420 XF=X-XI
4440 YA=YF*(YF-1)*.5
4460 YB=1-YF*YF
4480 YC=YF*(YF+1)*.5
4500 Y1=BN(XI-1,YI-1)*YA+BN(XI-1,YI)*YB+BN(XI-1,YI+1)*YC
4520 Y2=BN(XI,YI-1)*YA+BN(XI,YI)*YB+BN(XI,YI+1)*YC
4540 Y3=BN(XI+1,YI-1)*YA+BN(XI+1,YI)*YB+BN(XI+1,YI+1)*YC
4560 B(U)=Y1*XF*(XF-1)*.5+Y2*(1-XF*XF)+Y3*XF*(XF+1)*.5
4580 RETURN
5000 DT=DS
5020 LE=LS
5040 T=NS
5060 K=KS
5080 GOSUB 5260
5100 SI=L
5120 DT=DB
5140 LE=LL
5160 T=NL
5180 K=KL
5200 GOSUB 5260
5220 IL=L
5240 RETURN
5260 PS=LE/DT
5280 PL=DT/LE
5300 GOSUB 5500
5320 L=.025069*(DT*T)^2*KN/LE
5340 G=.556853+LOG(K) '
5360 IF T<10 THEN 5420
5380 H=.336*(1-2.5/T+3.8/(T*T))
5400 GOTO 5440
5420 H=H(T)
5440 DL=.0159593*DT*T*(G+H)
5460 L=L-DL
5480 RETURN
5500 IF LE>DT THEN 5620
5520 PS=100*PS
5540 PI=INT(PS)
5560 P=PS-PI
5580 LT=1
5600 GOTO 5700
5620 PL=100*PL
5640 PI=INT(PL)
5660 P=PL-PI
5680 LT=2
5700 KN=.5*KN(LT,PI-1)*P*(P-1)+KN(LT,PI)*(1-P*P)+.5*KN(LT,PI+1)*P*(P+1)
5720 RETURN
6000 CLS
6020 PRINT " LOADING WIRE TABLES"
6030 FOR X=10 TO 40
6050 FOR Y=1 TO 3
6060 READ W(Y,X)
6070 NEXT Y
6080 NEXT X
6090 ' THIS TABLE GIVES, FOR WIRE SIZES FROM 10 TO 40, THE
6100 ' DIAMETER AND THE RESISTANCE AND WEIGHT PER 1000 FEET.
6110 DATA .104,.9989,31.8
6111 DATA .0927,1.26,25.25
6112 DATA .0828,1.588,20.05
6113 DATA .074,2.003,15.90
6114 DATA .0661,2.525,12.60
6115 DATA .0591,3.184,10
6116 DATA .0528,4.016,7.93
6117 DATA .0471,5.064,6.275
6118 DATA .0421,6.385,4.980
6119 DATA .0377,8.051,3.955
6120 DATA .0338,10.15,3.135
6121 DATA .0302,12.8,2.490
6122 DATA .027,16.14,1.970
6123 DATA .0241,20.36,1.565
6124 DATA .0215,25.67,1.245
6125 DATA .0192,32.37,.988
6126 DATA .0171,40.81,.7845
6127 DATA .0153,51.47,.6220
6128 DATA .0136,64.9,.4940
6129 DATA .0122,81.83,.3915
6130 DATA .0108,103.2,.3105
6131 DATA .0097,130.1,.2465
6132 DATA .0087,164.1,.1960
6133 DATA .0077,206.9,.1550
6134 DATA .0069,260.9,.1230
6135 DATA .0062,329,.0980
6136 DATA .0055,414.8,.0776
6137 DATA .0049,532.1,.0616
6138 DATA .0044,659.6,.0488
6139 DATA .0039,831.8,.0387
6140 DATA .0035,1049,.0307
6150 W(1,0)=1E-08
6160 ' RETURN
7000
7020 PRINT " LOADINB BN, TABLE 29"
7030 FOR A9=20 TO 0 STEP -1
7035 GOTO 7050
7050 FOR R9=20 TO 0 STEP-1
7060 READ BN(A9,R9)
7062 GOTO 7070
7065
7070 NEXT R9
7080
7090 NEXT A9
7100
7110 ' RETURN
7120 ' LINES 7140-7550 ARE TABLE 29, PP125-126 FROM GROVER.
7130 ' BEGINNING OF TABLE 29
7140 DATA .84833,.87727,.89552,.91020,.92264,.93345,.94298,.95144,.95900,.96576
7150 DATA .97180,.97718,.98194,.98612,.98974,.99282,.99535,.99735,.99880,.99969
DATA 1
7160 DATA .86783,.88982,.90561,.91859,.92971,.93944,.94805,.95573,.96261,.96877
7170 DATA .97428,.97919,.98354,.98736,.99066,.99346,.99577,.99759,.99891,.99972
DATA 1
7180 DATA .88418,.90175,.91531,.92666,.93655,.94524,.95298,.95990,.96612,.97169
7190 DATA .97668,.98114,.98509,.98855,.99155,.99409,.99618,.99783,.99902,.99975
DATA 1
7200 DATA .89870,.91296,.92456,.93444,.94314,.95085,.95774,.96393,.96951,.97452
7210 DATA .97901,.98302,.98658,.98970,.99240,.99469,.99657,.99805,.99912,.99978
DATA 1
7220 DATA .91176,.92344,.93329,.94185,.94944,.95622,.96231,.96781,.97276,.97723
7230 DATA .98124,.98483,.98801,.99080,.99322,.99526,.99695,.99827,.99922,.99980
DATA 1
7240 DATA .92356,.93318,.94150,.94885,.95542,.96132,.96668,.97151,.97588,.97983
7250 DATA .98338,.98656,.98938,.99185,.99399,.99581,.99730,.99847,.99931,.99983
DATA 1
7260 DATA .93426,.94217,.94917,.95543,.96107,.96618,.97082,.97503,.97884,.98230
7270 DATA .98541,.98820,.99068,.99285,.99473,.99633,.99764,.99866,.99940,.99985
DATA 1
7280 DATA .94394,.95045,.95629,.96157,.96637,.97074,.97472,.97835,.98164,.98464
7290 DATA .98732,.98975,.99190,.99380,.99543,.99682,.99795,.99884,.99948,.99987
DATA 1
7300 DATA .95270,.95803,.96286,.96727,.97130,.97499,.97837,.98146,.98427,.98683
7310 DATA .98913,.99121,.99306,.99467,.99608,.99727,.99825,.99901,.99956,.99989
DATA 1
7320 DATA .96060,.96492,.96888,.97252,.97586,.97894,.98176,.98435,.98672,.98887
7330 DATA .99082,.99257,.99413,.99550,.99669,.99770,.99852,.99916,.99963,.99990
DATA 1
7340 DATA .96769,.97115,.97434,.97730,.98003,.98256,.98488,.98702,.98807,.99076
7350 DATA .99237,.99383,.99512,.99626,.99725,.99809,.99877,.99931,.99969,.99992
DATA 1
7360 DATA .97400,.97673,.97927,.98163,.98382,.98584,.98772,.98945,.99103,.99248
7370 DATA .99379,.99498,.99603,.99696,.99776,.99844,.99900,.99944,.99975,.99994
DATA 1
7380 DATA .97958,.98169,.98366,.98550,.98721,.98880,.99028,.99164,.99289,.99404
7390 DATA .99508,.99601,.99685,.99759,.99823,.99877,.99921,.99955,.99980,.99995
DATA 1
7400 DATA .98444,.98603,.98751,.9889,.9902,.99142,.99254,.99358,.99454,.99542
7410 DATA .99622,.99694,.99758,.99815,.99864,.99905,.99939,.99966,.99985,.99996
DATA 1
7420 DATA .98862,.98976,.99084,.99186,.9928,.99369,.99451,.99527,.99598,.99662
7430 DATA .99721,.99774,.99822,.99864,.99900,.99930,.99955,.99975,.99989,.99997
DATA 1
7440 DATA .99212,.99291,.99365,.99435,.995,.99561,.99618,.99671,.9972,.99765
7450 DATA .99806,.99843,.99876,.99905,.99930,.99952,.99969,.99982,.99992,.99998
DATA 1
7460 DATA .99498,.99547,.99594,.99638,.9968,.99719,.99755,.99789,.99821,.99849
7470 DATA .99875,.99899,.99920,.99939,.99955,.99969,.99980,.99989,.99995,.99999
DATA 1
7480 DATA .99718,.99746,.99772,.99797,.99820,.99842,.99862,.99881,.99899,.99915
7490 DATA .99930,.99943,.99955,.99966,.99975,.99982,.99989,.99994,.99997,.99999
DATA 1
7500 DATA .99875,.99887,.99899,.9991,.9992,.9993,.99939,.99947,.99955,.99962
7510 DATA .99969,.99975,.99980,.99985,.99989,.99992,.99995,.99997,.99999,1,1
7520 DATA .99969,.99972,.99975,.99977,.99980,.99982,.99985,.99987,.99989,.99991
7530 DATA .99992,.99994,.99995,.99996,.99997,.99998,.99999,.99999,1,1,1
7540 DATA 1,1,1,1,1,1,1,1,1,1
7550 DATA 1,1,1,1,1,1,1,1,1,1,1
7560 ' END OF TABLE 29
8000
8020 PRINT " LOADING BN, TABLE 30"
8030 FOR A9=10 TO 0 STEP -1
8050 FOR R9=10 TO 0 STEP -1
8060 READ LN(A9,R9)
8070 NEXT R9
8080
8090 NEXT A9
8100
8110 ' RETURN
8120 ' LINES 8140-8350 ARE TABLE 30, PASE 125 FROM GROVER'
8130 ' START OF TABLE 30
8140 DATA .84833,.85698,.86298,.86820,.87292,.87727
8150 DATA .88133,.88515,.88877,.89222,.89552
8160 DATA .85294,.86035,.86606,.87107,.87562,.87982
8170 DATA .88376,.88747,.89100,.89436,.89757
8180 DATA .85686,.86366,.86910,.87391,.87829,.88236
8190 DATA .88617,.88978,.89320,.89647,.89960
8200 DATA .86063,.86693,.87210,.87672,.88094,.88487
8210 DATA .88857,.89207,.89539,.89858,.90162
8220 DATA .86428,.87014,.87506,.87949,.80356,.88736
8230 DATA .89094,.89433,.89757,.90066,.90362
8240 DATA .86783,.87329,.87798,.88223,.88615,.88982
8250 DATA .89329,.89658,.89972,.90273,.90561
8260 DATA .87127,.87639,.88086,.88494,.88872,.89226
8270 DATA .89562,.89881,.90186,.90478,.90759
8280 DATA .87462,.87944,.88370,.88761,.89125,.89468
8290 DATA .89792,.90102,.90397,.90681,.90954
8300 DATA .87788,.88242,.88649,.89024,.89375,.89706
8310 DATA .90020,.90320,.90607,.90883,.91148
8320 DATA .88107,.88536,.88924,.89285,.89622,.89942
8330 DATA .90246,.90536,.90815,.91082,.91340
8340 DATA .88418,.88824,.89195,.89541,.89866,.90175
8350 DATA .90469,.90750,.91020,.91280,.91531
8360 ' END OF TABLE 30.
9010 PRINT " LOADING NAGAOKA'S CONSTANT"
9020 PRINT " FOR SHORT COILS"
9030 FOR X=1 TO 2
9040 FOR Y=0 TO 100
9060 READ KN(X,Y)
9070 NEXT Y
9080 IF X=2 THEN 9130
9090
9110 PRINT " LOADING NAGAOKA'S CONSTANT"
9120 PRINT " FOR LONG COILS"
9130 NEXT X
9140 FOR X=1 TO 10
9150 READ H(X)
9160 NEXT X
9170 CLS
9180 RETURN
9190 ' 9230 TO 9330 ARE TABLE 36, PAGES 144-5, FROM GROVER.
9200 ' 9360 TO 9460 ARE TABLE 37, PASES 146-7, FROM GROVER.
9210 ' 9490 IS LINES 1-10 OF TABLE 39 FROM GROVER.
9220 ' BEGINNING OF TABLE 36
9230 DATA 0,.03496,.061098,.083907,.104562,.123615,.141395,.158119,.173942
DATA .18898
9240 DATA .203314,.217044,.230200,.242842,.255011,.266744,.27807,.259019
DATA .299164,.309876
9250 DATA .319825,.329479,.338852,.34796,.356816,.365432,.373818,.381986
DATA .389944,.397703
9260 DATA .405269,.412650,.419856,.42589,.433762,.440474,.447036,.453450
DATA .459724,.46586
9270 DATA .471865,.477742,.483496,.489129,.494646,.500052,.505348,.510359
DATA .515628 ,.520617
9280 DATA .52551,.530310,.535018,.539637,.544170,.548620,.552988,.557278
DATA .561491,.565628
9290 DATA .569691,.573683,.577606,.581462,.585252,.588976,.592683,.596239
DATA .599780 ,.603263
9300 DATA .606689,.61006,.613376,.616639,.619850,.623011,.626122,.629185
DATA .6322,.63517
9310 DATA .638094,.640974,.643811,.646605,.649358,.65207,.654743,.657376
DATA .659972,.662532
9320 DATA .665054,.66754,.669991,.672408,.674792,.677142,.679460,.681747
DATA .684003,.686228
9330 DATA .688423
9340 ' END OF TABLE 36
9350 ' BEGINNING OF TABLE 37
9360 DATA 1.0000,.995769,.991562,.987381,.983224,.979092,.974985,.970903
DATA .966847,.962815
9370 DATA .958807,.954825,.950868,.946935,.943025,.939141,.935284,.93145
DATA .927639,.923854
9380 DATA .920093,.916356,.912643,.908954,.90529,.901649,.898033,.894440
DATA .890971,.887325
9390 DATA .883803,.880305,.876829,.873377,.869948,.866542,.863158,.859799
DATA .856461 ,.853146
9400 DATA .849853,.846583,.843335,.84011,.836906,.833723,.830563,.827424
DATA .824307,.821211
9410 DATA .818136,.815082,.812049,.809037,.806046,.803075,.800125,.797195
DATA .794285 ,.791395
9420 DATA .788525,.785675,.782844,.780032,.77724,.774467,.771713,.768978
DATA .766262,.763565
9430 DATA .760886,.758225,.755582,.752958,.750351,.747762,.745191,.742637
DATA .74010 ,.737581
9440 DATA .735079,.732593,.730126,.727675,.72524,.722821,.720419,.718033
DATA .715663,.713308
9450 DATA .710969,.708647,.706339,.704047,.70177,.699509,.697262,.695030
DATA .692813,.690611
9460 DATA .688423
9470 ' END OF TABLE 37
9480 ' START OF TABLE 39
9490 DATA 0,.1137,.1663,.1973,.2180,.2329,.2443,.2532,.2604,.2664
9500 ' END OF TABLE 39
10000 END
OOPS. See a couple of line wraps here which may be bad news. Perhaps
you can fix those, or I can try to edit it to eliminate them. Let
me know what you can do. By the way, this thing gives results which
are accurate to one part in 100,000 IF you know the input parameters
that well.
Ed