derivative4 Function

  function derivative4(f,dh,order)  result(df)
    real(8),dimension(:),intent(in) :: f
    real(8),intent(in)              :: dh
    integer,intent(in),optional     :: order
    real(8),dimension(size(f))      :: df
    integer                         :: i,L,order_
    L=size(f)
    order_=4;if(present(order))order_=order
    if(L < order_ + 1) stop "derivative4: L < order+1."
    select case(order_)
    case(2)
       df = derivF4_n222(f,dh)
    case(4)
       df = derivF4_n444(f,dh)
    case(6)
       df = derivF4_n666(f,dh)
    case default
       stop "derivative4: order is not 2,4 or 6"
    end select
  end function derivative4