Cocal: /home/uryu/Cocal_svn_checkinout/Cocal/code/Function/lagint

00001 function lagint_2nd(x,y,v)
00002   use phys_constant, only : long
00003   implicit none
00004   real(long) :: lagint_2nd
00005   real(long) :: x(2),y(2), v
00006   real(long) :: dx12, dx21
00007   real(long) :: xv1, xv2, wex1, wex2
00008 !
00009       dx12 = x(1) - x(2)
00010       dx21 = - dx12
00011       xv1 = v - x(1)
00012       xv2 = v - x(2)
00013       wex1 = xv2/dx12
00014       wex2 = xv1/dx21
00015 !
00016       lagint_2nd = wex1*y(1) + wex2*y(2)
00017 !
00018 end function lagint_2nd

/home/uryu/Cocal_svn_checkinout/Cocal/code/Function/lagint_2nd.f90