SUBROUTINE LagrangeLambda(x,i,xvals,Nvals,Val_Lambda) USE Numeric_Kinds_Module IMPLICIT NONE INTEGER k INTEGER, INTENT(IN):: i, Nvals REAL(KIND=WP_Kind), INTENT(IN):: x, xvals(Nvals) REAL(KIND=WP_Kind), INTENT(OUT):: Val_Lambda VAl_Lambda=1.0_WP_Kind !IF(i.ne.1)THEN DO k=1,Nvals IF(k.ne.i)THEN Val_Lambda=Val_Lambda*(x-xvals(k))/(xvals(i)-xvals(k)) ENDIF ENDDO !ELSE ! DO k=2,Nvals ! Val_Lambda=Val_Lambda*(x-xvals(k))/(xvals(i)-xvals(k)) ! ENDDO !ENDIF RETURN END SUBROUTINE LagrangeLambda