100 ' save "hme-1.bas",a
110 CLS 3: SCREEN 3,0,0,1:COLOR ,,,,2:CONSOLE 0,24,0,1,0
120 LINE (80,5)-(560,90),5,BF
130 LINE (82,7)-(558,88),15,BF
140 PRINT
150 PRINT "               ハーモニックミーン方程式　による表面エネルギーの計算"
160 PRINT "               ジオメトリックミーン方程式による表面エネルギーの計算"
170 PRINT "               ヤング・デュプレ方程式　　による界面エネルギーの計算"
180 PRINT "                      wu p.615-618    96/2/6     by T.Tokoro"
190 PRINT
200 DEF FNCG(X1,X2)=((X1*Q2-X2*Q1)/A)^2                 'FOR GME
210 DEF FNCQ(G,T)=.5*G*(1+COS(T/180*3.14159))           '    HME
220 COLOR 6
230 PRINT "         水とＭＩのコンタクト角を入力して下さい。Water,MI(deg.)   ";
240 COLOR 7
250 INPUT "",T1,T2
260 G1=72.8:G2=50.8      'surface tentions of water and MI
270 D1=22.1:D2=44.1      'dispersion component of water and MI
280 P1=G1-D1:P2=G2-D2
290 Q1=.5*FNCQ(G1,T1)
300 Q2=.5*FNCQ(G2,T2)
310 S1=D1+P1-Q1
320 S2=D2+P2-Q2
330 R1=P1*(D1-Q1)
340 R2=P2*(D2-Q2)
350 U1=D1*(P1-Q1)
360 U2=D2*(P2-Q2)
370 V1=Q1*D1*P1
380 V2=Q2*D2*P2
390 A=1/S1/S2
400 W=R1/S1-R2/S2
410 X=A*(R1*U2-R2*U1)+V2/S2-V1/S1
420 Y=A*(U1*V2-U2*V1)
430 X=-.5*X/W
440 Y=(X^2)-Y/W
450 IF Y<0 THEN PRINT "y<0",Y:GOTO 110
460 Y=SQR(Y)
470 D3=X-Y
480 P3=(V1-R1*D3)/(S1*D3+U1)
490 G3=D3+P3
500 WSL=G1*(1+COS(T1*3.14159/180))
510 GSL=G3-G1*COS(T1*3.14159/180)
520 '
530 COLOR 3
540 PRINT "ＨＭＥによる計算結果"
550 GOSUB *GAMEN
560 '
570 D3=X+Y
580 P3=(V2-R2*D3)/(S2*D3+U2)
590 G3=D3+P3
600 WSL=G1*(1+COS(T1*3.14159/180))
610 GSL=G3-G1*COS(T1*3.14159/180)
620 '
630 GOSUB *GAMEN
640 '
650 D1=21.8:D2=49.5
660 P1=G1-D1:P2=G2-D2
670 Q1=FNCQ(G1,T1)
680 Q2=FNCQ(G2,T2)
690 R1=SQR(D1)
700 R2=SQR(D2)
710 S1=SQR(P1)
720 S2=SQR(P2)
730 A=S1*R2-S2*R1
740 D3=FNCG(S1,S2)
750 P3=FNCG(R1,R2)
760 G3=D3+P3
770 WSL=G1*(1+COS(T1*3.14159/180))
780 GSL=G3-G1*COS(T1*3.14159/180)
790 '
800 COLOR 3
810 PRINT "ＧＭＥによる計算結果"
820 GOSUB *GAMEN
830 '
840 INPUT "          Hit Any Key!",DAMMY$
850 GOTO 100
860 '
870 *GAMEN
880 COLOR 7
890 PRINT "    表面エネルギー　γs＝γsd ＋ γsp  ：";:COLOR 4
900 PRINT USING "  ####.# = ####.#  ＋  ####.#";G3;D3;P3
910 PRINT: COLOR 7
920 PRINT "    界面エネルギー　ＷSL，γSL         ：";:COLOR 5
930 PRINT USING "  ####.# ,  ####.#     ";WSL,GSL
940 PRINT:COLOR 7
950 RETURN