msta1 Function

public function msta1(x, Mp) result(nn)

Determine the starting point for backward recurrence such that the magnitude of Jn(x) at that point is about 10^(-MP)

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: x

Argument of Jn(x)
integer, intent(in) :: Mp

Value of magnitude

Return Value integer

Starting point


Calls

proc~~msta1~~CallsGraph proc~msta1 specfun_module::msta1 proc~envj specfun_module::envj proc~msta1->proc~envj

Called by

proc~~msta1~~CalledByGraph proc~msta1 specfun_module::msta1 proc~cikna specfun_module::cikna proc~cikna->proc~msta1 proc~ciknb specfun_module::ciknb proc~ciknb->proc~msta1 proc~cikva specfun_module::cikva proc~cikva->proc~msta1 proc~cikvb specfun_module::cikvb proc~cikvb->proc~msta1 proc~cjynb specfun_module::cjynb proc~cjynb->proc~msta1 proc~cjyva specfun_module::cjyva proc~cjyva->proc~msta1 proc~cjyvb specfun_module::cjyvb proc~cjyvb->proc~msta1 proc~ikna specfun_module::ikna proc~ikna->proc~msta1 proc~iknb specfun_module::iknb proc~iknb->proc~msta1 proc~ikv specfun_module::ikv proc~ikv->proc~msta1 proc~jyna specfun_module::jyna proc~jyna->proc~msta1 proc~jynbh specfun_module::jynbh proc~jynbh->proc~msta1 proc~jyv specfun_module::jyv proc~jyv->proc~msta1 proc~lamn specfun_module::lamn proc~lamn->proc~msta1 proc~lamv specfun_module::lamv proc~lamv->proc~msta1 proc~rctj specfun_module::rctj proc~rctj->proc~msta1 proc~sphi specfun_module::sphi proc~sphi->proc~msta1 proc~sphj specfun_module::sphj proc~sphj->proc~msta1 proc~ch12n specfun_module::ch12n proc~ch12n->proc~ciknb proc~ch12n->proc~cjynb proc~jynb specfun_module::jynb proc~jynb->proc~jynbh proc~jyndd specfun_module::jyndd proc~jyndd->proc~jynbh proc~rmn1 specfun_module::rmn1 proc~rmn1->proc~sphj proc~jyzo specfun_module::jyzo proc~jyzo->proc~jyndd proc~mtu12 specfun_module::mtu12 proc~mtu12->proc~jynb proc~rmn2so specfun_module::rmn2so proc~rmn2so->proc~rmn1 proc~rswfo specfun_module::rswfo proc~rswfo->proc~rmn1 proc~rswfo->proc~rmn2so proc~rswfp specfun_module::rswfp proc~rswfp->proc~rmn1

Source Code

   function msta1(x,Mp) result(nn)

      real(wp),intent(in) :: x !! Argument of Jn(x)
      integer,intent(in) :: Mp !! Value of magnitude
      integer :: nn !! Starting point

      real(wp) :: a0 , f , f0 , f1
      integer :: it , n0 , n1

      a0 = abs(x)
      n0 = int(1.1_wp*a0) + 1
      f0 = envj(n0,a0) - Mp
      n1 = n0 + 5
      f1 = envj(n1,a0) - Mp
      do it = 1 , 20
         nn = n1 - (n1-n0)/(1.0_wp-f0/f1)
         f = envj(nn,a0) - Mp
         if ( abs(nn-n1)<1 ) exit
         n0 = n1
         f0 = f1
         n1 = nn
         f1 = f
      enddo

   end function msta1