derivativeN Function

public function derivativeN(f, dh, n) result(df)

Arguments

Type IntentOptional Attributes Name
real(kind=8), intent(in), dimension(:) :: f
real(kind=8), intent(in) :: dh
integer, intent(in) :: n

Return Value real(kind=8), dimension(size(f))


Source Code

  function derivativeN(f,dh,n)  result(df)
    real(8),dimension(:),intent(in) :: f
    real(8),intent(in)              :: dh
    integer,intent(in)              :: n
    real(8),dimension(size(f))      :: df,tmp,dtmp
    integer                         :: i,L
    L=size(f)
    if(L < n + 1) stop "derivative4: L < order+1."
    tmp=f
    do i=1,n
       dtmp = derivF_n666(tmp,dh)
       tmp=dtmp
    enddo
    df=dtmp
  end function derivativeN