UnscaledLobattoZeros

Returns the zeros of UnscaledLobatto polynomials.

Interface

܀ Interface
️܀ See example
↢

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

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

Interface​

Interface