Line | |
---|
1 | subroutine eint(a,ee,ek) |
---|
2 | c evaluate elliptic integrals |
---|
3 | c----------------------------------------------------------------------- |
---|
4 | c |
---|
5 | include 'ucom.h' |
---|
6 | c |
---|
7 | c----------------------------------------------------------------------- |
---|
8 | c* |
---|
9 | b=1.0-a |
---|
10 | c=log(b) |
---|
11 | ee=1.0+b*(.4630106-0.2452740*c+b*(0.1077857-0.04125321*c)) |
---|
12 | ek=1.38629436-.5*c+b*(0.1119697-0.1213486*c+b*(.07253230-.02887472 |
---|
13 | 1*c)) |
---|
14 | return |
---|
15 | end |
---|
16 | c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
---|
Note: See
TracBrowser
for help on using the repository browser.