| [x,rx,c,bi,dbi]=dyairy(x,rx,c,bi,dbi); |
function [x,rx,c,bi,dbi]=dyairy(x,rx,c,bi,dbi);
%***BEGIN PROLOGUE DYAIRY
%***SUBSIDIARY
%***PURPOSE Subsidiary to DBESJ and DBESY
%***LIBRARY SLATEC
%***TYPE doubleprecision (YAIRY-S, DYAIRY-D)
%***AUTHOR Amos, D. E., (SNLA)
% Daniel, S. L., (SNLA)
%***DESCRIPTION
%
% DYAIRY computes the Airy function BI(X)
% and its derivative DBI(X) for DASYJY
%
% INPUT
%
% X - Argument, computed by DASYJY, X unrestricted
% RX - RX=SQRT(ABS(X)), computed by DASYJY
% C - C=2.*(ABS(X)**1.5)/3., computed by DASYJY
%
% OUTPUT
% BI - Value of function BI(X)
% DBI - Value of the derivative DBI(X)
%
%***SEE ALSO DBESJ, DBESY
%***ROUTINES CALLED (NONE)
%***REVISION HISTORY (YYMMDD)
% 750101 DATE WRITTEN
% 890531 Changed all specific intrinsics to generic. (WRB)
% 891214 Prologue converted to Version 4.0 format. (BAB)
% 900328 Added TYPE section. (WRB)
% 910408 Updated the AUTHOR section. (WRB)
%***end PROLOGUE DYAIRY
%
persistent aa ax bb bjn bjp bk1 bk2 bk3 bk4 con1 con2 con3 cv d1 d2 daa dbb dbjn dbjp dbk1 dbk2 dbk3 dbk4 e1 e2 ex f1 f2 firstCall fpi12 i j m1 m1d m2 m2d m3 m3d m4d n1 n1d n2 n2d n3 n3d n4d rtrx s1 s2 spi12 t tc temp1 temp2 tt ; if isempty(firstCall),firstCall=1;end;
if isempty(i), i=0; end;
if isempty(j), j=0; end;
if isempty(m1), m1=0; end;
if isempty(m1d), m1d=0; end;
if isempty(m2), m2=0; end;
if isempty(m2d), m2d=0; end;
if isempty(m3), m3=0; end;
if isempty(m3d), m3d=0; end;
if isempty(m4d), m4d=0; end;
if isempty(n1), n1=0; end;
if isempty(n1d), n1d=0; end;
if isempty(n2), n2=0; end;
if isempty(n2d), n2d=0; end;
if isempty(n3), n3=0; end;
if isempty(n3d), n3d=0; end;
if isempty(n4d), n4d=0; end;
if isempty(aa), aa=zeros(1,14); end;
if isempty(ax), ax=0; end;
if isempty(bb), bb=zeros(1,14); end;
if isempty(bjn), bjn=zeros(1,19); end;
if isempty(bjp), bjp=zeros(1,19); end;
if isempty(bk1), bk1=zeros(1,20); end;
if isempty(bk2), bk2=zeros(1,20); end;
if isempty(bk3), bk3=zeros(1,20); end;
if isempty(bk4), bk4=zeros(1,14); end;
if isempty(con1), con1=0; end;
if isempty(con2), con2=0; end;
if isempty(con3), con3=0; end;
if isempty(cv), cv=0; end;
if isempty(daa), daa=zeros(1,14); end;
if isempty(dbb), dbb=zeros(1,14); end;
if isempty(dbjn), dbjn=zeros(1,19); end;
if isempty(dbjp), dbjp=zeros(1,19); end;
if isempty(dbk1), dbk1=zeros(1,21); end;
if isempty(dbk2), dbk2=zeros(1,20); end;
if isempty(dbk3), dbk3=zeros(1,20); end;
if isempty(dbk4), dbk4=zeros(1,14); end;
if isempty(d1), d1=0; end;
if isempty(d2), d2=0; end;
if isempty(ex), ex=0; end;
if isempty(e1), e1=0; end;
if isempty(e2), e2=0; end;
if isempty(fpi12), fpi12=0; end;
if isempty(f1), f1=0; end;
if isempty(f2), f2=0; end;
if isempty(rtrx), rtrx=0; end;
if isempty(spi12), spi12=0; end;
if isempty(s1), s1=0; end;
if isempty(s2), s2=0; end;
if isempty(t), t=0; end;
if isempty(tc), tc=0; end;
if isempty(temp1), temp1=0; end;
if isempty(temp2), temp2=0; end;
if isempty(tt), tt=0; end;
if firstCall, n1 =[20]; end;
if firstCall, n2 =[19]; end;
if firstCall, n3=[14]; end;
if firstCall, m1 =[18]; end;
if firstCall, m2 =[17]; end;
if firstCall, m3=[12]; end;
if firstCall, n1d =[21]; end;
if firstCall, n2d =[20]; end;
if firstCall, n3d =[19]; end;
if firstCall, n4d=[14]; end;
if firstCall, m1d =[19]; end;
if firstCall, m2d =[18]; end;
if firstCall, m3d =[17]; end;
if firstCall, m4d=[12]; end;
if firstCall, fpi12 =[1.30899693899575d+00]; end;
if firstCall, spi12 =[1.83259571459405d+00]; end;
if firstCall, con1 =[6.66666666666667d-01]; end;
if firstCall, con2 =[7.74148278841779d+00]; end;
if firstCall, con3=[3.64766105490356d-01]; end;
if firstCall, bk1(1) =[2.43202846447449d+00]; end;
if firstCall, bk1(2) =[2.57132009754685d+00]; end;
if firstCall, bk1(3) =[1.02802341258616d+00]; end;
if firstCall, bk1(4) =[3.41958178205872d-01]; end;
if firstCall, bk1(5) =[8.41978629889284d-02]; end;
if firstCall, bk1(6) =[1.93877282587962d-02]; end;
if firstCall, bk1(7)=[3.92687837130335d-03]; end;
if firstCall, bk1(8) =[6.83302689948043d-04]; end;
if firstCall, bk1(9) =[1.14611403991141d-04]; end;
if firstCall, bk1(10) =[1.74195138337086d-05]; end;
if firstCall, bk1(11) =[2.41223620956355d-06]; end;
if firstCall, bk1(12) =[3.24525591983273d-07]; end;
if firstCall, bk1(13) =[4.03509798540183d-08]; end;
if firstCall, bk1(14) =[4.70875059642296d-09]; end;
if firstCall, bk1(15) =[5.35367432585889d-10]; end;
if firstCall, bk1(16) =[5.70606721846334d-11]; end;
if firstCall, bk1(17) =[5.80526363709933d-12]; end;
if firstCall, bk1(18) =[5.76338988616388d-13]; end;
if firstCall, bk1(19) =[5.42103834518071d-14]; end;
if firstCall, bk1(20)=[4.91857330301677d-15]; end;
if firstCall, bk2(1) =[5.74830555784088d-01]; end;
if firstCall, bk2(2) =[-6.91648648376891d-03]; end;
if firstCall, bk2(3) =[1.97460263052093d-03]; end;
if firstCall, bk2(4) =[-5.24043043868823d-04]; end;
if firstCall, bk2(5) =[1.22965147239661d-04]; end;
if firstCall, bk2(6) =[-2.27059514462173d-05]; end;
if firstCall, bk2(7)=[2.23575555008526d-06]; end;
if firstCall, bk2(8) =[4.15174955023899d-07]; end;
if firstCall, bk2(9) =[-2.84985752198231d-07]; end;
if firstCall, bk2(10) =[8.50187174775435d-08]; end;
if firstCall, bk2(11) =[-1.70400826891326d-08]; end;
if firstCall, bk2(12) =[2.25479746746889d-09]; end;
if firstCall, bk2(13) =[-1.09524166577443d-10]; end;
if firstCall, bk2(14) =[-3.41063845099711d-11]; end;
if firstCall, bk2(15) =[1.11262893886662d-11]; end;
if firstCall, bk2(16) =[-1.75542944241734d-12]; end;
if firstCall, bk2(17) =[1.36298600401767d-13]; end;
if firstCall, bk2(18) =[8.76342105755664d-15]; end;
if firstCall, bk2(19) =[-4.64063099157041d-15]; end;
if firstCall, bk2(20)=[7.78772758732960d-16]; end;
if firstCall, bk3(1) =[5.66777053506912d-01]; end;
if firstCall, bk3(2) =[2.63672828349579d-03]; end;
if firstCall, bk3(3) =[5.12303351473130d-05]; end;
if firstCall, bk3(4) =[2.10229231564492d-06]; end;
if firstCall, bk3(5) =[1.42217095113890d-07]; end;
if firstCall, bk3(6) =[1.28534295891264d-08]; end;
if firstCall, bk3(7)=[7.28556219407507d-10]; end;
if firstCall, bk3(8) =[-3.45236157301011d-10]; end;
if firstCall, bk3(9) =[-2.11919115912724d-10]; end;
if firstCall, bk3(10) =[-6.56803892922376d-11]; end;
if firstCall, bk3(11) =[-8.14873160315074d-12]; end;
if firstCall, bk3(12) =[3.03177845632183d-12]; end;
if firstCall, bk3(13) =[1.73447220554115d-12]; end;
if firstCall, bk3(14) =[1.67935548701554d-13]; end;
if firstCall, bk3(15) =[-1.49622868806719d-13]; end;
if firstCall, bk3(16) =[-5.15470458953407d-14]; end;
if firstCall, bk3(17) =[8.75741841857830d-15]; end;
if firstCall, bk3(18) =[7.96735553525720d-15]; end;
if firstCall, bk3(19) =[-1.29566137861742d-16]; end;
if firstCall, bk3(20)=[-1.11878794417520d-15]; end;
if firstCall, bk4(1) =[4.85444386705114d-01]; end;
if firstCall, bk4(2) =[-3.08525088408463d-03]; end;
if firstCall, bk4(3) =[6.98748404837928d-05]; end;
if firstCall, bk4(4) =[-2.82757234179768d-06]; end;
if firstCall, bk4(5) =[1.59553313064138d-07]; end;
if firstCall, bk4(6) =[-1.12980692144601d-08]; end;
if firstCall, bk4(7)=[9.47671515498754d-10]; end;
if firstCall, bk4(8) =[-9.08301736026423d-11]; end;
if firstCall, bk4(9) =[9.70776206450724d-12]; end;
if firstCall, bk4(10) =[-1.13687527254574d-12]; end;
if firstCall, bk4(11) =[1.43982917533415d-13]; end;
if firstCall, bk4(12) =[-1.95211019558815d-14]; end;
if firstCall, bk4(13) =[2.81056379909357d-15]; end;
if firstCall, bk4(14)=[-4.26916444775176d-16]; end;
if firstCall, bjp(1) =[1.34918611457638d-01]; end;
if firstCall, bjp(2) =[-3.19314588205813d-01]; end;
if firstCall, bjp(3) =[5.22061946276114d-02]; end;
if firstCall, bjp(4) =[5.28869112170312d-02]; end;
if firstCall, bjp(5) =[-8.58100756077350d-03]; end;
if firstCall, bjp(6) =[-2.99211002025555d-03]; end;
if firstCall, bjp(7)=[4.21126741969759d-04]; end;
if firstCall, bjp(8) =[8.73931830369273d-05]; end;
if firstCall, bjp(9) =[-1.06749163477533d-05]; end;
if firstCall, bjp(10) =[-1.56575097259349d-06]; end;
if firstCall, bjp(11) =[1.68051151983999d-07]; end;
if firstCall, bjp(12) =[1.89901103638691d-08]; end;
if firstCall, bjp(13) =[-1.81374004961922d-09]; end;
if firstCall, bjp(14) =[-1.66339134593739d-10]; end;
if firstCall, bjp(15) =[1.42956335780810d-11]; end;
if firstCall, bjp(16) =[1.10179811626595d-12]; end;
if firstCall, bjp(17) =[-8.60187724192263d-14]; end;
if firstCall, bjp(18) =[-5.71248177285064d-15]; end;
if firstCall, bjp(19)=[4.08414552853803d-16]; end;
if firstCall, bjn(1) =[6.59041673525697d-02]; end;
if firstCall, bjn(2) =[-4.24905910566004d-01]; end;
if firstCall, bjn(3) =[2.87209745195830d-01]; end;
if firstCall, bjn(4) =[1.29787771099606d-01]; end;
if firstCall, bjn(5) =[-4.56354317590358d-02]; end;
if firstCall, bjn(6) =[-1.02630175982540d-02]; end;
if firstCall, bjn(7)=[2.50704671521101d-03]; end;
if firstCall, bjn(8) =[3.78127183743483d-04]; end;
if firstCall, bjn(9) =[-7.11287583284084d-05]; end;
if firstCall, bjn(10) =[-8.08651210688923d-06]; end;
if firstCall, bjn(11) =[1.23879531273285d-06]; end;
if firstCall, bjn(12) =[1.13096815867279d-07]; end;
if firstCall, bjn(13) =[-1.46234283176310d-08]; end;
if firstCall, bjn(14) =[-1.11576315688077d-09]; end;
if firstCall, bjn(15) =[1.24846618243897d-10]; end;
if firstCall, bjn(16) =[8.18334132555274d-12]; end;
if firstCall, bjn(17) =[-8.07174877048484d-13]; end;
if firstCall, bjn(18) =[-4.63778618766425d-14]; end;
if firstCall, bjn(19)=[4.09043399081631d-15]; end;
if firstCall, aa(1) =[-2.78593552803079d-01]; end;
if firstCall, aa(2) =[3.52915691882584d-03]; end;
if firstCall, aa(3) =[2.31149677384994d-05]; end;
if firstCall, aa(4) =[-4.71317842263560d-06]; end;
if firstCall, aa(5) =[1.12415907931333d-07]; end;
if firstCall, aa(6) =[2.00100301184339d-08]; end;
if firstCall, aa(7) =[-2.60948075302193d-09]; end;
if firstCall, aa(8) =[3.55098136101216d-11]; end;
if firstCall, aa(9) =[3.50849978423875d-11]; end;
if firstCall, aa(10) =[-5.83007187954202d-12]; end;
if firstCall, aa(11) =[2.04644828753326d-13]; end;
if firstCall, aa(12) =[1.10529179476742d-13]; end;
if firstCall, aa(13) =[-2.87724778038775d-14]; end;
if firstCall, aa(14)=[2.88205111009939d-15]; end;
if firstCall, bb(1) =[-4.90275424742791d-01]; end;
if firstCall, bb(2) =[-1.57647277946204d-03]; end;
if firstCall, bb(3) =[9.66195963140306d-05]; end;
if firstCall, bb(4) =[-1.35916080268815d-07]; end;
if firstCall, bb(5) =[-2.98157342654859d-07]; end;
if firstCall, bb(6) =[1.86824767559979d-08]; end;
if firstCall, bb(7) =[1.03685737667141d-09]; end;
if firstCall, bb(8) =[-3.28660818434328d-10]; end;
if firstCall, bb(9) =[2.57091410632780d-11]; end;
if firstCall, bb(10) =[2.32357655300677d-12]; end;
if firstCall, bb(11) =[-9.57523279048255d-13]; end;
if firstCall, bb(12) =[1.20340828049719d-13]; end;
if firstCall, bb(13) =[2.90907716770715d-15]; end;
if firstCall, bb(14)=[-4.55656454580149d-15]; end;
if firstCall, dbk1(1) =[2.95926143981893d+00]; end;
if firstCall, dbk1(2) =[3.86774568440103d+00]; end;
if firstCall, dbk1(3) =[1.80441072356289d+00]; end;
if firstCall, dbk1(4) =[5.78070764125328d-01]; end;
if firstCall, dbk1(5) =[1.63011468174708d-01]; end;
if firstCall, dbk1(6) =[3.92044409961855d-02]; end;
if firstCall, dbk1(7) =[7.90964210433812d-03]; end;
if firstCall, dbk1(8) =[1.50640863167338d-03]; end;
if firstCall, dbk1(9) =[2.56651976920042d-04]; end;
if firstCall, dbk1(10) =[3.93826605867715d-05]; end;
if firstCall, dbk1(11) =[5.81097771463818d-06]; end;
if firstCall, dbk1(12)=[7.86881233754659d-07]; end;
if firstCall, dbk1(13) =[9.93272957325739d-08]; end;
if firstCall, dbk1(14) =[1.21424205575107d-08]; end;
if firstCall, dbk1(15) =[1.38528332697707d-09]; end;
if firstCall, dbk1(16) =[1.50190067586758d-10]; end;
if firstCall, dbk1(17) =[1.58271945457594d-11]; end;
if firstCall, dbk1(18) =[1.57531847699042d-12]; end;
if firstCall, dbk1(19) =[1.50774055398181d-13]; end;
if firstCall, dbk1(20) =[1.40594335806564d-14]; end;
if firstCall, dbk1(21)=[1.24942698777218d-15]; end;
if firstCall, dbk2(1) =[5.49756809432471d-01]; end;
if firstCall, dbk2(2) =[9.13556983276901d-03]; end;
if firstCall, dbk2(3) =[-2.53635048605507d-03]; end;
if firstCall, dbk2(4) =[6.60423795342054d-04]; end;
if firstCall, dbk2(5) =[-1.55217243135416d-04]; end;
if firstCall, dbk2(6) =[3.00090325448633d-05]; end;
if firstCall, dbk2(7) =[-3.76454339467348d-06]; end;
if firstCall, dbk2(8) =[-1.33291331611616d-07]; end;
if firstCall, dbk2(9) =[2.42587371049013d-07]; end;
if firstCall, dbk2(10) =[-8.07861075240228d-08]; end;
if firstCall, dbk2(11) =[1.71092818861193d-08]; end;
if firstCall, dbk2(12)=[-2.41087357570599d-09]; end;
if firstCall, dbk2(13) =[1.53910848162371d-10]; end;
if firstCall, dbk2(14) =[2.56465373190630d-11]; end;
if firstCall, dbk2(15) =[-9.88581911653212d-12]; end;
if firstCall, dbk2(16) =[1.60877986412631d-12]; end;
if firstCall, dbk2(17) =[-1.20952524741739d-13]; end;
if firstCall, dbk2(18) =[-1.06978278410820d-14]; end;
if firstCall, dbk2(19) =[5.02478557067561d-15]; end;
if firstCall, dbk2(20)=[-8.68986130935886d-16]; end;
if firstCall, dbk3(1) =[5.60598509354302d-01]; end;
if firstCall, dbk3(2) =[-3.64870013248135d-03]; end;
if firstCall, dbk3(3) =[-5.98147152307417d-05]; end;
if firstCall, dbk3(4) =[-2.33611595253625d-06]; end;
if firstCall, dbk3(5) =[-1.64571516521436d-07]; end;
if firstCall, dbk3(6) =[-2.06333012920569d-08]; end;
if firstCall, dbk3(7) =[-4.27745431573110d-09]; end;
if firstCall, dbk3(8) =[-1.08494137799276d-09]; end;
if firstCall, dbk3(9) =[-2.37207188872763d-10]; end;
if firstCall, dbk3(10) =[-2.22132920864966d-11]; end;
if firstCall, dbk3(11) =[1.07238008032138d-11]; end;
if firstCall, dbk3(12)=[5.71954845245808d-12]; end;
if firstCall, dbk3(13) =[7.51102737777835d-13]; end;
if firstCall, dbk3(14) =[-3.81912369483793d-13]; end;
if firstCall, dbk3(15) =[-1.75870057119257d-13]; end;
if firstCall, dbk3(16) =[6.69641694419084d-15]; end;
if firstCall, dbk3(17) =[2.26866724792055d-14]; end;
if firstCall, dbk3(18) =[2.69898141356743d-15]; end;
if firstCall, dbk3(19) =[-2.67133612397359d-15]; end;
if firstCall, dbk3(20)=[-6.54121403165269d-16]; end;
if firstCall, dbk4(1) =[4.93072999188036d-01]; end;
if firstCall, dbk4(2) =[4.38335419803815d-03]; end;
if firstCall, dbk4(3) =[-8.37413882246205d-05]; end;
if firstCall, dbk4(4) =[3.20268810484632d-06]; end;
if firstCall, dbk4(5) =[-1.75661979548270d-07]; end;
if firstCall, dbk4(6) =[1.22269906524508d-08]; end;
if firstCall, dbk4(7) =[-1.01381314366052d-09]; end;
if firstCall, dbk4(8) =[9.63639784237475d-11]; end;
if firstCall, dbk4(9) =[-1.02344993379648d-11]; end;
if firstCall, dbk4(10) =[1.19264576554355d-12]; end;
if firstCall, dbk4(11) =[-1.50443899103287d-13]; end;
if firstCall, dbk4(12)=[2.03299052379349d-14]; end;
if firstCall, dbk4(13) =[-2.91890652008292d-15]; end;
if firstCall, dbk4(14)=[4.42322081975475d-16]; end;
if firstCall, dbjp(1) =[1.13140872390745d-01]; end;
if firstCall, dbjp(2) =[-2.08301511416328d-01]; end;
if firstCall, dbjp(3) =[1.69396341953138d-02]; end;
if firstCall, dbjp(4) =[2.90895212478621d-02]; end;
if firstCall, dbjp(5) =[-3.41467131311549d-03]; end;
if firstCall, dbjp(6) =[-1.46455339197417d-03]; end;
if firstCall, dbjp(7) =[1.63313272898517d-04]; end;
if firstCall, dbjp(8) =[3.91145328922162d-05]; end;
if firstCall, dbjp(9) =[-3.96757190808119d-06]; end;
if firstCall, dbjp(10) =[-6.51846913772395d-07]; end;
if firstCall, dbjp(11) =[5.98707495269280d-08]; end;
if firstCall, dbjp(12)=[7.44108654536549d-09]; end;
if firstCall, dbjp(13) =[-6.21241056522632d-10]; end;
if firstCall, dbjp(14) =[-6.18768017313526d-11]; end;
if firstCall, dbjp(15) =[4.72323484752324d-12]; end;
if firstCall, dbjp(16) =[3.91652459802532d-13]; end;
if firstCall, dbjp(17) =[-2.74985937845226d-14]; end;
if firstCall, dbjp(18) =[-1.95036497762750d-15]; end;
if firstCall, dbjp(19)=[1.26669643809444d-16]; end;
if firstCall, dbjn(1) =[-1.88091260068850d-02]; end;
if firstCall, dbjn(2) =[-1.47798180826140d-01]; end;
if firstCall, dbjn(3) =[5.46075900433171d-01]; end;
if firstCall, dbjn(4) =[1.52146932663116d-01]; end;
if firstCall, dbjn(5) =[-9.58260412266886d-02]; end;
if firstCall, dbjn(6) =[-1.63102731696130d-02]; end;
if firstCall, dbjn(7) =[5.75364806680105d-03]; end;
if firstCall, dbjn(8) =[7.12145408252655d-04]; end;
if firstCall, dbjn(9) =[-1.75452116846724d-04]; end;
if firstCall, dbjn(10) =[-1.71063171685128d-05]; end;
if firstCall, dbjn(11) =[3.24435580631680d-06]; end;
if firstCall, dbjn(12)=[2.61190663932884d-07]; end;
if firstCall, dbjn(13) =[-4.03026865912779d-08]; end;
if firstCall, dbjn(14) =[-2.76435165853895d-09]; end;
if firstCall, dbjn(15) =[3.59687929062312d-10]; end;
if firstCall, dbjn(16) =[2.14953308456051d-11]; end;
if firstCall, dbjn(17) =[-2.41849311903901d-12]; end;
if firstCall, dbjn(18) =[-1.28068004920751d-13]; end;
if firstCall, dbjn(19)=[1.26939834401773d-14]; end;
if firstCall, daa(1) =[2.77571356944231d-01]; end;
if firstCall, daa(2) =[-4.44212833419920d-03]; end;
if firstCall, daa(3) =[8.42328522190089d-05]; end;
if firstCall, daa(4) =[2.58040318418710d-06]; end;
if firstCall, daa(5) =[-3.42389720217621d-07]; end;
if firstCall, daa(6) =[6.24286894709776d-09]; end;
if firstCall, daa(7)=[2.36377836844577d-09]; end;
if firstCall, daa(8) =[-3.16991042656673d-10]; end;
if firstCall, daa(9) =[4.40995691658191d-12]; end;
if firstCall, daa(10) =[5.18674221093575d-12]; end;
if firstCall, daa(11) =[-9.64874015137022d-13]; end;
if firstCall, daa(12) =[4.90190576608710d-14]; end;
if firstCall, daa(13) =[1.77253430678112d-14]; end;
if firstCall, daa(14)=[-5.55950610442662d-15]; end;
if firstCall, dbb(1) =[4.91627321104601d-01]; end;
if firstCall, dbb(2) =[3.11164930427489d-03]; end;
if firstCall, dbb(3) =[8.23140762854081d-05]; end;
if firstCall, dbb(4) =[-4.61769776172142d-06]; end;
if firstCall, dbb(5) =[-6.13158880534626d-08]; end;
if firstCall, dbb(6) =[2.87295804656520d-08]; end;
if firstCall, dbb(7)=[-1.81959715372117d-09]; end;
if firstCall, dbb(8) =[-1.44752826642035d-10]; end;
if firstCall, dbb(9) =[4.53724043420422d-11]; end;
if firstCall, dbb(10) =[-3.99655065847223d-12]; end;
if firstCall, dbb(11) =[-3.24089119830323d-13]; end;
if firstCall, dbb(12) =[1.62098952568741d-13]; end;
if firstCall, dbb(13) =[-2.40765247974057d-14]; end;
if firstCall, dbb(14)=[1.69384811284491d-16]; end;
firstCall=0;
%***FIRST EXECUTABLE STATEMENT DYAIRY
ax = abs(x);
rx = sqrt(ax);
c = con1.*ax.*rx;
if( x<0.0d0 )
%
if( c>5.0d0 )
%
rtrx = sqrt(rx);
t = 10.0d0./c - 1.0d0;
tt = t + t;
j = fix(n3);
f1 = aa(j);
e1 = bb(j);
f2 = 0.0d0;
e2 = 0.0d0;
for i = 1 : m3;
j = fix(j - 1);
temp1 = f1;
temp2 = e1;
f1 = tt.*f1 - f2 + aa(j);
e1 = tt.*e1 - e2 + bb(j);
f2 = temp1;
e2 = temp2;
end; i = fix(m3+1);
temp1 = t.*f1 - f2 + aa(1);
temp2 = t.*e1 - e2 + bb(1);
cv = c - fpi12;
bi =(temp1.*cos(cv)+temp2.*sin(cv))./rtrx;
j = fix(n4d);
f1 = daa(j);
e1 = dbb(j);
f2 = 0.0d0;
e2 = 0.0d0;
for i = 1 : m4d;
j = fix(j - 1);
temp1 = f1;
temp2 = e1;
f1 = tt.*f1 - f2 + daa(j);
e1 = tt.*e1 - e2 + dbb(j);
f2 = temp1;
e2 = temp2;
end; i = fix(m4d+1);
temp1 = t.*f1 - f2 + daa(1);
temp2 = t.*e1 - e2 + dbb(1);
cv = c - spi12;
dbi =(temp1.*cos(cv)-temp2.*sin(cv)).*rtrx;
else;
t = 0.4d0.*c - 1.0d0;
tt = t + t;
j = fix(n2);
f1 = bjp(j);
e1 = bjn(j);
f2 = 0.0d0;
e2 = 0.0d0;
for i = 1 : m2;
j = fix(j - 1);
temp1 = f1;
temp2 = e1;
f1 = tt.*f1 - f2 + bjp(j);
e1 = tt.*e1 - e2 + bjn(j);
f2 = temp1;
e2 = temp2;
end; i = fix(m2+1);
bi =(t.*e1-e2+bjn(1)) - ax.*(t.*f1-f2+bjp(1));
j = fix(n3d);
f1 = dbjp(j);
e1 = dbjn(j);
f2 = 0.0d0;
e2 = 0.0d0;
for i = 1 : m3d;
j = fix(j - 1);
temp1 = f1;
temp2 = e1;
f1 = tt.*f1 - f2 + dbjp(j);
e1 = tt.*e1 - e2 + dbjn(j);
f2 = temp1;
e2 = temp2;
end; i = fix(m3d+1);
dbi = x.*x.*(t.*f1-f2+dbjp(1)) +(t.*e1-e2+dbjn(1));
return;
end;
elseif( c>8.0d0 ) ;
%
rtrx = sqrt(rx);
t = 16.0d0./c - 1.0d0;
tt = t + t;
j = fix(n1);
f1 = bk3(j);
f2 = 0.0d0;
for i = 1 : m1;
j = fix(j - 1);
temp1 = f1;
f1 = tt.*f1 - f2 + bk3(j);
f2 = temp1;
end; i = fix(m1+1);
s1 = t.*f1 - f2 + bk3(1);
j = fix(n2d);
f1 = dbk3(j);
f2 = 0.0d0;
for i = 1 : m2d;
j = fix(j - 1);
temp1 = f1;
f1 = tt.*f1 - f2 + dbk3(j);
f2 = temp1;
end; i = fix(m2d+1);
d1 = t.*f1 - f2 + dbk3(1);
tc = c + c;
ex = exp(c);
if( tc>35.0d0 )
bi = ex.*s1./rtrx;
dbi = ex.*rtrx.*d1;
return;
else;
t = 10.0d0./c - 1.0d0;
tt = t + t;
j = fix(n3);
f1 = bk4(j);
f2 = 0.0d0;
for i = 1 : m3;
j = fix(j - 1);
temp1 = f1;
f1 = tt.*f1 - f2 + bk4(j);
f2 = temp1;
end; i = fix(m3+1);
s2 = t.*f1 - f2 + bk4(1);
bi =(s1+exp(-tc).*s2)./rtrx;
bi = bi.*ex;
j = fix(n4d);
f1 = dbk4(j);
f2 = 0.0d0;
for i = 1 : m4d;
j = fix(j - 1);
temp1 = f1;
f1 = tt.*f1 - f2 + dbk4(j);
f2 = temp1;
end; i = fix(m4d+1);
d2 = t.*f1 - f2 + dbk4(1);
dbi = rtrx.*(d1+exp(-tc).*d2);
dbi = dbi.*ex;
return;
end;
elseif( x>2.5d0 ) ;
rtrx = sqrt(rx);
t =(x+x-con2).*con3;
tt = t + t;
j = fix(n1);
f1 = bk2(j);
f2 = 0.0d0;
for i = 1 : m1;
j = fix(j - 1);
temp1 = f1;
f1 = tt.*f1 - f2 + bk2(j);
f2 = temp1;
end; i = fix(m1+1);
bi =(t.*f1-f2+bk2(1))./rtrx;
ex = exp(c);
bi = bi.*ex;
j = fix(n2d);
f1 = dbk2(j);
f2 = 0.0d0;
for i = 1 : m2d;
j = fix(j - 1);
temp1 = f1;
f1 = tt.*f1 - f2 + dbk2(j);
f2 = temp1;
end; i = fix(m2d+1);
dbi =(t.*f1-f2+dbk2(1)).*rtrx;
dbi = dbi.*ex;
return;
else;
t =(x+x-2.5d0).*0.4d0;
tt = t + t;
j = fix(n1);
f1 = bk1(j);
f2 = 0.0d0;
for i = 1 : m1;
j = fix(j - 1);
temp1 = f1;
f1 = tt.*f1 - f2 + bk1(j);
f2 = temp1;
end; i = fix(m1+1);
bi = t.*f1 - f2 + bk1(1);
j = fix(n1d);
f1 = dbk1(j);
f2 = 0.0d0;
for i = 1 : m1d;
j = fix(j - 1);
temp1 = f1;
f1 = tt.*f1 - f2 + dbk1(j);
f2 = temp1;
end; i = fix(m1d+1);
dbi = t.*f1 - f2 + dbk1(1);
return;
end;
end
%DECK DY
|
|
