10 REM **** OBOE ACOUSTIC LENGTH - FROM CONES ****                  USE RUN 1000 TO CALCULATE U-BEND
20 NS = 6: REM NO. OF SECTIONS
30 H = 0: REM MOUTHHOLE HEIGHT, RELATIVE TO DIAMETER THERE *

50 REM THE FOLLOWING DATA ARE THE LENGTHS OF THE NS SECTIONS
55  DATA 22.9,22,20,5.1,.2,241.2
REM 60 DATA 124.8,40,.0001,.0001,.0001,241.2
70 REM  THE FOLLOWING ARE THE DIAMETERS OF THE NS+1 CROSS
REM 80 DATA 1E-10,3.074,4.06,4.06,4.06,4.06,10
81 DATA 5.5 ,2.37,3.37,3.95,4.06,4.06,10

82 SCREEN 11: CLS
83 INPUT "DO YOU RESULTS PRINTED (Y OR N)"; RP$: CLS
84 IF RP$ = "; Y; " THEN LPRINT "; MOUTH; HEIGHT = "; H
90 DIM X(NS + 1): DIM D(NS + 1): DIM R(NS + 1)
93 PI = 3.14159
95 POKE 16553, 255
100 FOR N = 1 TO NS: READ X(N): NEXT N
110 FOR N = 1 TO NS + 1: READ D(N): NEXT N
120 FOR N = 1 TO NS
125 IF D(N + 1) = D(N) THEN R(N) = 1000000!: GOTO 140
130 R(N) = X(N) * D(N) / (D(N + 1) - D(N))
140 NEXT N
141 VIEW (0, 0)-(639, 240)
142 WINDOW (0, 0)-(410, 10.5)
143 LINE (0, 0)-(410, 0)
144 FOR N = 1 TO NS
145 DIST(N) = DIST(N - 1) + X(N): PRINT N, DIST(N)
146 LINE (DIST(N - 1), D(N))-(DIST(N), D(N + 1))
147 LINE (DIST(N - 1), 0)-(DIST(N - 1), D(N))
152 NEXT N
190 LOCATE 19, 1: PRINT " N  DIAM LENGTH DISTANCE CONE APEX": PRINT
192 IF RP$ = "Y" THEN LPRINT CHR$(30); : LPRINT " N", "DIAM", "LENGTHS", "DISTANCE", "CONE APEX": LPRINT
194 FOR N = 1 TO NS + 1
196 PRINT USING "##  ##.##  ###.#  ###.#  ####.#"; N; D(N); X(N); DIST(N); R(N)
197 IF RP$ = "Y" THEN LPRINT N, D(N), X(N), DIST(N), R(N)
198 NEXT N
199 WHILE INKEY$ = "": WEND
200 SCREEN 11: CLS 0
202 VIEW (280, 440)-(600, 200)
204 WINDOW (46, 370)-(78, 420)
206 REM  FL IS FREQUENCY LEVEL WHERE C4 =48, A4=57
208 PRINT : PRINT "NOTE   FREQ   LENGTH, MM"
210 IF RP$ = "Y" THEN LPRINT : LPRINT "NOTE   FREQ   LENGTH, MM"
211 CLS
212 LINE (46, 370)-(78, 370): LINE (46, 370)-(46, 420)
218  FOR FL = 48 TO 76 STEP 4
220 F = 16.352 * EXP(FL / 17.3123)
230 K = F / 55386!: REM *  ASSUME VELOCITY OF SOUND =348,000 MM/SEC
300 B = 0
310 FOR N = 1 TO NS
320 X = ATN(1 / (1 / K / R(N) - B)) / K
324 IF X < 0 THEN X = X + PI / K: GOTO 324
330 X = X + X(N)
340 B = 1 / K / (R(N) + X(N)) - 1 / TAN(K * X)
350 NEXT N
360 L = X
365 IF L < 350 THEN L = L + PI / K: GOTO 365
370 IN = INT(FL / 12)
380 P = (FL - 12 * IN) / 4
390 IF P = 0 THEN C$ = "C"
391 IF P = 1 THEN C$ = "E"
392 IF P = 2 THEN C$ = "G#"
400 PRINT C$; : PRINT TAB(7); USING "####.#   ###.#"; F; L
402 IF RP$ = "Y" THEN LPRINT C$; : LPRINT TAB(7); USING "####.#   ###.#"; F; L
404 PSET (FL, L)
410 NEXT FL
420 FOR FL = 48 TO 76 STEP 4: LINE (FL, 370)-(FL, 372)
421 IN = INT(FL / 12)
422 P = (FL - 12 * IN) / 4
424 IF P = 0 THEN C$ = "C"
426 IF P = 1 THEN C$ = "E"
428 IF P = 2 THEN C$ = "G#"
452 LOCATE 29, 5 * FL / 4 - 22: PRINT C$;
454 NEXT FL
500 FOR L = 370 TO 420 STEP 10: LINE (46, L)-(47, L)
510  LOCATE 139 - L * .3, 31: PRINT L; : NEXT L

700 END
2000 REM - FORMULA FOR EFFECTIVE CORK DISTANCE
2010 INPUT "DIAM AT MOUTH HOLE, DIAM AT CORK "; DH, DC
2020 INPUT " ACTUAL CORK DISTANCE"; CD
2030 DD = DC / DH
2040 CN = CD * (1 + DD + DD * DD) / 3
2050 PRINT "EFFECTIVE CORK DISTANCE FOR CONICAL FRUSTRUM ="; CN
2060 END