source: PSPA/parmelaPSPA/trunk/eint.f @ 427

Last change on this file since 427 was 12, checked in by lemeur, 12 years ago

parmela pspa initial

File size: 515 bytes
Line 
1      subroutine eint(a,ee,ek)
2c     evaluate elliptic integrals
3c-----------------------------------------------------------------------
4c
5      include 'ucom.h'
6c
7c-----------------------------------------------------------------------
8c*
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
16c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
Note: See TracBrowser for help on using the repository browser.