Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=8), | intent(in), | dimension(:) | :: | f | ||
real(kind=8), | intent(in) | :: | dh | |||
integer, | intent(in) | :: | n |
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