]>
Commit | Line | Data |
---|---|---|
4e9e3152 | 1 | function eksarp(x,v2,x1,x116,ptm,indpi) |
2 | IMPLICIT DOUBLE PRECISION (A-H,O-Z) | |
3 | dimension ptm(10) | |
4 | ||
5 | z=v2 | |
6 | xx=x-x1 | |
7 | if (x.le.x116) then | |
8 | eksarp=ptm(1)+ptm(2)*z+ptm(3)*z**2+ptm(4)*xx+ptm(5)*xx**2 | |
9 | else | |
10 | z1= log(x116/x1) | |
11 | xx1=x116-x1 | |
12 | ff0=ptm(1)+ptm(2)*z1+ptm(3)*z1**2+ptm(4)*xx1+ptm(5)*xx1**2 | |
13 | xx16=x-x116 | |
14 | z16= log(x/x116) | |
15 | qexp=19.d0+(indpi-1)*(indpi-2)*8.d0/2.d0 | |
16 | eksarp=ff0+ptm(6)*xx16+ptm(7)*xx16**2+ptm(8)*xx16**3 | |
17 | + +ptm(9)*xx16**qexp+ptm(10)*z16 | |
18 | endif | |
19 | 11 return | |
20 | end | |
21 | ||
22 |