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