mirror of
https://github.com/PDP-10/stacken.git
synced 2026-03-01 01:19:17 +00:00
120 lines
3.0 KiB
QBasic
120 lines
3.0 KiB
QBasic
170 DIM Z(3),H(3),W(3)
|
||
190 V4=6.02E23
|
||
210 B4=.511
|
||
230 D4=2.8E-13
|
||
250 Z4=2
|
||
270 R4=1E-8
|
||
290 P4=1E-5
|
||
310 Z(1)=29
|
||
320 Z(2)=47
|
||
330 Z(3)=79
|
||
350 W(1)=63.5
|
||
360 W(2)=107.9
|
||
370 W(3)=197
|
||
390 H(1)=8.96
|
||
400 H(2)=10.5
|
||
410 H(3)=19.3
|
||
420 E4=5.3
|
||
430 S=3
|
||
440 T=1
|
||
450 A4=1
|
||
460 L4=10
|
||
470 PRINT"INSTRUCTIONS (0=NO,1=YES)";
|
||
480 INPUT I
|
||
490 IF I=0 THEN 600
|
||
500 IF I<>1 THEN 470
|
||
510 PRINT"SCATR2 PLOTS THEORETICAL PREDICTIONS OF THE FRACTION OF"
|
||
520 PRINT"ALPHA PARTICLES SCATTERED AT ANGLES FROM 10 TO 180"
|
||
530 PRINT"DEGREES BY THREE DIFFERENT HEAVY ATOM MODELS."
|
||
540 PRINT"INPUTS ARE: ALPHA PARTICLE ENERGY (BETWEEN 2 AND 8 MEV),"
|
||
550 PRINT"TARGET MATERIAL (1=COPPER,2=SILVER,3=GOLD),TARGET THICKNESS"
|
||
560 PRINT"(BETWEEN .5 AND 2 MICROMETERS), DETECTOR SIZE"
|
||
570 PRINT"(BETWEEN .5 AND 2 SQ CM), AND DISTANCE FROM TARGET TO"
|
||
580 PRINT"DETECTORS (BETWEEN 10 AND 40 CM)."
|
||
590 PRINT
|
||
600 PRINT"ENERGY=";E4;
|
||
610 INPUT I
|
||
620 IF (I-2)*(8-I)>=0 THEN 650
|
||
630 PRINT "2<=ENERGY<=8"
|
||
640 GOTO 600
|
||
650 E4=I
|
||
660 PRINT"MATERIAL=";S;
|
||
670 INPUT I
|
||
680 IF (I-1)*(I-2)*(I-3)=0 THEN 710
|
||
690 PRINT"1=CU,2=AG,3=AU"
|
||
700 GOTO 660
|
||
710 S=I
|
||
720 PRINT "THICKNESS=";T;
|
||
730 INPUT I
|
||
740 IF (I-.5)*(2-I)>=0 THEN 770
|
||
750 PRINT".5<=THICKNESS<=2"
|
||
760 GOTO 720
|
||
770 T=I
|
||
780 T4=1E-5*T
|
||
790 PRINT"SIZE=";A4;
|
||
800 INPUT I
|
||
810 IF (I-.5)*(2-I)>=0 THEN 840
|
||
820 PRINT".5<=SIZE<=2"
|
||
830 GOTO 790
|
||
840 A4=I
|
||
850 PRINT"DISTANCE=";L4;
|
||
860 INPUT I
|
||
870 IF (I-10)*(40-I)>=0 THEN 900
|
||
880 PRINT"10<=DISTANCE<=40"
|
||
890 GOTO 850
|
||
900 L4=I
|
||
920 N=T4*V4*H(S)/W(S)
|
||
940 D=Z4*Z(S)*B4*D4/E4
|
||
960 PRINT"MODEL (1=HARD SPHERES, 2=THOMSON, 3=RUTHERFORD)";
|
||
970 INPUT M
|
||
980 IF (M-1)*(M-2)*(M-3)<>0 THEN 960
|
||
1000 PRINT
|
||
1010 IF M>1 THEN 1050
|
||
1020 PRINT"HARD SPHERE MODEL"
|
||
1030 F=P4*N*R4^2*A4/(4*L4^2)
|
||
1040 F1=8+10*(LOG(F)/LOG(10)+9)
|
||
1050 IF M<>2 THEN 1080
|
||
1060 PRINT"THOMSON MODEL"
|
||
1070 F1=7
|
||
1080 IF M<3 THEN 1100
|
||
1090 PRINT "RUTHERFORD MODEL"
|
||
1100 PRINT E4;"MEV ALPHA'S ";T;"MICRON ";
|
||
1110 IF S>1 THEN 1130
|
||
1120 PRINT"COPPER";
|
||
1130 IF S<>2 THEN 1150
|
||
1140 PRINT"SILVER";
|
||
1150 IF S<3 THEN 1170
|
||
1160 PRINT"GOLD";
|
||
1170 PRINT" ";A4;"SQ CM COUNTERS ";L4;"CM AWAY"
|
||
1180 PRINT
|
||
1190 PRINT TAB(15);"FRACTION OF PARTICLES SCATTERED (LOG SCALE)"
|
||
1200 GOSUB 1310
|
||
1210 FOR I=10 TO 180 STEP 10
|
||
1220 IF M<3 THEN 1270
|
||
1240 S1=SIN(3.1416*I/360)^4
|
||
1250 F=N*D*D*A4/(16*L4^2*S1)
|
||
1260 F1=8+10*(LOG(F)/LOG(10)+9)
|
||
1270 PRINT I;TAB(7);"I";
|
||
1280 IF ABS(F1-38)>31 THEN 1370
|
||
1290 PRINT TAB(F1);"*"
|
||
1300 GOTO 1400
|
||
1310 PRINT " 1E-9 1E-8 1E-7 1E-6 1E-5";
|
||
1320 PRINT " 1E-4 1E-3"
|
||
1330 PRINT "ANGLE:";
|
||
1340 PRINT TAB(7);"I . . ..I . . ..I . . ..I . . ..I . . ..";
|
||
1350 PRINT "I . . ..I"
|
||
1360 RETURN
|
||
1370 IF F1<7 THEN 1390
|
||
1380 PRINT TAB(57);
|
||
1390 PRINT "(";F;")"
|
||
1400 NEXT I
|
||
1410 GOSUB 1340
|
||
1420 PRINT
|
||
1430 PRINT"NEXT PREDICTION (1=NEW MODEL, 2=NEW ENERGY, ETC.,";
|
||
1440 PRINT"3=END PROGRAM";
|
||
1450 INPUT K
|
||
1460 IF K=1 THEN 960
|
||
1470 IF K=2 THEN 600
|
||
1480 IF K<>3 THEN 1430
|
||
1490 END
|
||
*U*,>>B$$BB$$BB$$BL |