Skip to main content

UnscaledLobattoMonomialExpansion

Returns all the monomial expansion of a UnscaledLobatto polynomials.

Interface

܀ Interface
️܀ See example
↢

INTERFACE
  MODULE PURE FUNCTION UnscaledLobattoMonomialExpansion(n) RESULT(ans)
    INTEGER(I4B), INTENT(IN) :: n
    REAL(DFP) :: ans(1:n + 1)
  END FUNCTION UnscaledLobattoMonomialExpansion
END INTERFACE

program main
  use easifembase
  implicit none
  integer( i4b ) :: n
  real( dfp ), allocatable :: coeff(:, :)
  type(string) :: astr
  n = 5
  astr = MdEncode( LobattoMonomialExpansion( n=n ) )
  call display( astr%chars(), "" )
end program main

x0	x1	x2	x3	x4	x5
0	0.7955	-0	-2.6517	0	1.8562

n : is the order of the polynomial
ans(:) contains the coefficient of monomials for polynomial order=n

Interface