Skip to main content

Jacobi1D example 7

This example shows the usage of [[Jacobi1D_]] class. We test Zeros function in this routine, which returns the zeros of Jacobi polynomial.

Modules and classes​

  • [[Jacobi1D_]]

Usage​

PROGRAM main
use easifemBase
use easifemClasses
implicit none
type(Jacobi1D_) :: obj
real( dfp ), allocatable :: x( : )
integer( i4b ) :: n
real( dfp ), parameter :: alpha=0.0_DFP, beta=0.0_DFP

n=1​

  n = 1
  obj=Jacobi1D(varname="x", n=n, alpha=alpha, beta=beta)
  x = obj%zeros()
  call display( x, "zeros for n="//tostring(n), orient="ROW" )

n=2​

  n = 2
  obj=Jacobi1D(varname="x", n=n, alpha=alpha, beta=beta)
  x = obj%zeros()
  call display( x, "zeros for n="//tostring(n), orient="ROW" )

n=3​

  n = 3
  obj=Jacobi1D(varname="x", n=n, alpha=alpha, beta=beta)
  x = obj%zeros()
  call display( x, "zeros for n="//tostring(n), orient="ROW" )

n=4​

  n = 4
  obj=Jacobi1D(varname="x", n=n, alpha=alpha, beta=beta)
  x = obj%zeros()
  call display( x, "zeros for n="//tostring(n), orient="ROW" )
END PROGRAM main