\START\
\COMMENT=by James Marshall, calculates flux for stellar model
\NAME=KAPPA
\FILE=kappa.85p
Disp "Find F(0) for model."
Disp "Edit T,Pe,KAPR,TAUR,"
Disp "XH in program."
Input "\LC-lambda\ [Ang]: ",\LC-lambda\
{4271.0,4365.2,4394.6,4422.6,4454.4,4491.3,4533.0,4578.1,4624.8,4671.3,\#\
4716.3,4759.4,4802.5,4850.2,4908.7,4984.1,5082.8,5213.4,5390.7,5494.9,\#\
5618.2,5759.5,5920.4,6103.3,6311.0,6545.4,6813.9,7130.9,7451.1,7732.6,\#\
7975.0,8185.8,8373.3,8544.0,8701.9}\->\T
5040/T\->\\LC-theta\
\10^\{-1.3602,-1.0453,-0.9392,-0.8339,-0.7278,-0.6250,-0.5118,-0.4022,\#\
-0.2922,-0.1826,-0.0738,0.0340,0.1415,0.2504,0.3628,0.4811,0.6079,0.7472,\#\
0.9092,1.0012,1.1078,1.2317,1.3753,1.5401,1.7259,1.9299,2.1525,2.3985,\#\
2.6295,2.8188,2.9732,3.1012,3.2107,3.3069,3.3932}\->\P
0.72\->\XH
1.67\E\-24\->\mH
XH/mH\->\Hpg
1.044\E\-26\->\\LC-alpha\0
log (e^1)\->\loge
13.60\->\I
0.0010968\->\R
1.99654\->\a0
-1.18267\E\-5\->\a1
2.64243\E\-6\->\a2
-4.40524\E\-10\->\a3
3.23992\E\-14\->\a4
-1.39568\E\-18\->\a5
2.78701\E\-23\->\a6
-2.2763-1.6850*log \LC-lambda\+0.76661*(log \LC-lambda\)\^2\-0.0533464*\#\
(log \LC-lambda\)^3\->\f0
15.2827-9.2864*log \LC-lambda\+1.99381*(log \LC-lambda\)\^2\-0.142631*\#\
(log \LC-lambda\)^3\->\f1
-197.789+190.266*log \LC-lambda\-67.9775*(log \LC-lambda\)\^2\+10.6913*\#\
(log \LC-lambda\)^3-0.625151*(log \LC-lambda\)^4\->\f2
(a0+a1*\LC-lambda\+a2*\LC-lambda\\^2\+a3*\LC-lambda\^3+a4*\LC-lambda\^4+\#\
a5*\LC-lambda\^5+a6*\LC-lambda\^6)*1\E\-18\->\\LC-alpha\bf
1\E\-26*\10^\(f0+f1*log \LC-theta\+f2*(log \LC-theta\)\^2\)\->\\LC-alpha\ff
1.2398\E\4/\LC-lambda\\->\X\LC-lambda\
1-\10^\(-X\LC-lambda\*\LC-theta\)\->\SE
(h*c/(\LC-lambda\*1\E\-10))\J\\>\\eV\\->\E\LC-lambda\
1\->\n0
While (E\LC-lambda\<(13.60-(13.60*(1-1/n0\^2\))))
n0+1\->\n0
End
I*(1-1/(n0+3)\^2\)\->\X3
1-(0.3456*(\LC-lambda\*R)^(-1/3)*(\LC-lambda\*R/(n0)\^2\-0.5))\->\gn0
1-(0.3456*(\LC-lambda\*R)^(-1/3)*(\LC-lambda\*R/(n0+1)\^2\-0.5))\->\gn0p1
1-(0.3456*(\LC-lambda\*R)^(-1/3)*(\LC-lambda\*R/(n0+2)\^2\-0.5))\->\gn0p2
1+0.3456*(\LC-lambda\*R)^(-1/3)*(loge/(\LC-theta\*X\LC-lambda\)+0.5)\->\gf
13.60*(1-1/(n0)\^2\)\->\Xn0
13.60*(1-1/(n0+1)\^2\)\->\Xn0p1
13.60*(1-1/(n0+2)\^2\)\->\Xn0p2
(gn0/(n0)^3)*\10^\(-\LC-theta\*Xn0)+(gn0p1/(n0+1)^3)*\10^\(-\LC-theta\*\#\
Xn0p1)+(gn0p2/(n0+2)^3)*\10^\(-\LC-theta\*Xn0p2)\->\gSUM
\LC-alpha\0*\LC-lambda\^3*(gSUM+((loge/(2*\LC-theta\*I))*(\10^\(-X3*\#\
\LC-theta\)-\10^\(-I*\LC-theta\))))\->\HIbf
\LC-alpha\0*\LC-lambda\^3*gf*(loge/(2*\LC-theta\*I))*\10^\(-\LC-theta\*I)\#\
\->\HIff
4.158\E\-10*\LC-alpha\bf*P*\LC-theta\^(5/2)*\10^\(0.754*\LC-theta\)\->\HMbf
P*\LC-alpha\ff\->\HMff
0.30103\->\c0
-0.00001\->\c1
\10^\(c0+c1*log \LC-theta\)\->\u0
1\->\u1
(\10^\(-0.1762)*\10^\(-\LC-theta\*I)*T^(5/2)*(u1/u0))\->\\LC-phi\
(1+\LC-phi\/P)\^-1\\->\PHI
(((HIbf+HIff+HMbf)*SE)+HMff)*PHI*Hpg\->\KAP
{.2812\E\-2,.5207\E\-2,.6399\E\-2,.7868\E\-2,.9670\E\-2,.1188\E\-1,\#\
.1460\E\-1,.1793\E\-1,.2204\E\-1,.2708\E\-1,.3329\E\-1,.4093\E\-1,\#\
.5032\E\-1,.6185\E\-1,.7602\E\-1,.9348\E\-1,.1152,.1429,.1814,.2081,\#\
.2446,.2971,.3751,.4940,.6780,.9655,1.4286,2.2198,3.3876,4.8292,6.5899,\#\
8.3341,1.0365\E\1,1.2598\E\1,1.5047\E\1}\->\KAPR
{0,.0001,.0002,.0003,.0004,.0006,.0010,.0016,.0025,.0040,.0063,.0100,\#\
.0158,.0251,.0398,.0631,.1000,.1585,.2512,.3162,.3981,.5012,.6310,\#\
.7943,1.000,1.2589,1.5849,1.9953,2.5119,3.1623,3.9811,5.0119,6.3096,7.9433,\#\
10.0000}\->\TAUR
0\->\TAU(1)
For(Q,2,35,1)
(KAP(Q)/KAPR(Q))*(TAUR(Q)-TAUR(Q-1))+TAU(Q-1)\->\TAU(Q)
End
1\->\CNT1
1\->\CNT2
While (0.397>TAU(CNT1))
CNT1+1\->\CNT1
End
While (2.723>TAU(CNT2))
CNT2+1\->\CNT2
End
inter(TAU(CNT1-1),T(CNT1-1),TAU(CNT1),T(CNT1),0.397)\->\Tp3
inter(TAU(CNT2-1),T(CNT2-1),TAU(CNT2),T(CNT2),2.723)\->\T2p7
(1.1904397\E\27/\LC-lambda\^5)/(e^(1.438769\E\8/(\LC-lambda\*Tp3))-1)\->\Bp3
(1.1904397\E\27/\LC-lambda\^5)/(e^(1.438769\E\8/(\LC-lambda\*T2p7))-1)\->\B2p7
0.8839*Bp3+0.1161*B2p7\->\F
\pi\*F\->\FLUX
Disp "F(0)=",FLUX
Disp "KAP,TAU also found."
\STOP\

    Source: geocities.com/dronak/progs

               ( geocities.com/dronak)