Skip to main content

Polynomial1D example 6

  • This example shows the usage of [[Polynomial1D_]] class.
  • In this example, we test subtraction operators.

Modules and classes​

  • [[Polynomial1D_]]

Usage​

PROGRAM main
use easifemBase
use easifemClasses
implicit none
type(Polynomial1D_) :: f1, f2, f3
real(dfp), allocatable :: coeff( : )
integer(i4b), allocatable :: degree( : )

!!! note "Mono- Mono"

  f1=Monomial1D(1,"x") -  Monomial1D(2,"x")
  call f1%display( 'f(x)=' )

!!! example "result"

f(x)=x1−x2f(x)=x^1- x^2

!!! note "Mono- Scalar"

  f1 = Monomial1D(1,"x") -  1
  call f1%display("f(x)=")
  !!
  f1 = Monomial1D(1,"x") -  1.0
  call f1%display("f(x)=")
  !!

!!! note "Scalar- Mono"

  f1 = 1 -  Monomial1D(1,"x")
  call f1%display("f(x)=")
  !!
  f1 = 1.0 -  Monomial1D(1,"x")
  call f1%display("f(x)=")
  !!

!!! note "Poly- Mono"

  f2 = f1 -  Monomial1D(2,"x")
  call f2%display("f(x)=")
  !!

!!! note "Mono- Poly"

  f2 = Monomial1D(2,"x") -  f1
  call f2%display("f(x)=")
  !!

!!! note "Poly- Poly"

  f2 = f2 -  f1
  call f2%display("f(x)=")
  !!

!!! note "Poly- Scalar"

  f2 = f1 -  1
  call f2%display("f(x)=")
  !!
  f2 = f1 -  1.0
  call f2%display("f(x)=")
  !!

!!! note "Scalar- Poly"

  f2 = 1 -  f1
  call f2%display("f(x)=")
  !!
  f2 = 1.0 -  f1
  call f2%display("f(x)=")
  !!
END PROGRAM main