UnscaledLobattoZeros

Returns the zeros of UnscaledLobatto polynomials.

Interface

܀ Interface

INTERFACE
  MODULE FUNCTION UnscaledLobattoZeros(n) RESULT(ans)
    INTEGER(I4B), INTENT(IN) :: n
    !! order of UnscaledLobatto polynomial, should be greater than equal to 2
    REAL(DFP) :: ans(n)
    !!
  END FUNCTION UnscaledLobattoZeros
END INTERFACE

️܀ See example

This example shows the usage of UnscaledLobattoZeros method.

• n should be greater than 2.

program main
  use easifembase
  implicit none
  integer( i4b ) :: n
  real( dfp ), allocatable :: pt( : )
  type(string) :: msg, astr

  n = 2; call callme

-1	1

  n = 3; call callme

-1	0	1

  n = 4; call callme

-1	-0.44721	0.44721	1

  n = 5; call callme

-1	-0.65465	4.29344E-17	0.65465	1

  contains
  subroutine callme
    pt = UnscaledLobattoZeros( n=n )
    msg = "Zeros of J(x), n = " &
        & // tostring( n )
    call display(msg%chars())
    astr = MdEncode( pt )
    call display( astr%chars(), "" )
  end subroutine

end program main

↢