Skip to main content

ChebyshevFirst1D example 8

This example shows the usage of [[ChebyshevFirst1D_]] class. We test GaussRadauQuadrature function in this routine, which returns the GaussRadauQuadrature points for Chebyshev1 polynomial.

Modules and classes

  • [[ChebyshevFirst1D_]]

Usage

PROGRAM main
use easifemBase
use easifemClasses
implicit none
type(ChebyshevFirst1D_) :: obj
real( dfp ), allocatable :: x( :, : )
integer( i4b ) :: n
real( dfp ), parameter :: a = -1.0_DFP

n=1

  n = 1
obj=ChebyshevFirst1D(varname="x", n=n)
x = obj%GaussRadauQuadrature(a)
call display( x, "pt | wt for n="//tostring(n) )

n=2

  n = 2
obj=ChebyshevFirst1D(varname="x", n=n)
x = obj%GaussRadauQuadrature(a)
call display( x, "pt | wt for n="//tostring(n) )

n=3

  n = 3
obj=ChebyshevFirst1D(varname="x", n=n)
x = obj%GaussRadauQuadrature(a)
call display( x, "pt | wt for n="//tostring(n) )

n=4

  n = 4
obj=ChebyshevFirst1D(varname="x", n=n)
x = obj%GaussRadauQuadrature(a)
call display( x, "pt | wt for n="//tostring(n) )
END PROGRAM main