function [ x, y, z, w ] = ld2702 ( ) %*****************************************************************************80; % %% LD2702 computes the 2702 point Lebedev angular grid. % % Modified: % % 14 September 2010 % % Author: % % Dmitri Laikov % % Reference: % % Vyacheslav Lebedev, Dmitri Laikov, % A quadrature formula for the sphere of the 131st % algebraic order of accuracy, % Russian Academy of Sciences Doklady Mathematics, % Volume 59, Number 3, 1999, pages 477-481. % % Parameters: % % Output, real X(N), Y(N), Z(N), W(N), the coordinates % and weights of the points. % n = 0; x = zeros(2702,1); y = zeros(2702,1); z = zeros(2702,1); w = zeros(2702,1); a = 0.0; b = 0.0; v = 0.2998675149888161E-04; [ n, x, y, z, w ] = gen_oh ( 1, n, a, b, v, x, y, z, w ); v = 0.4077860529495355E-03; [ n, x, y, z, w ] = gen_oh ( 3, n, a, b, v, x, y, z, w ); a = 0.2065562538818703E-01; v = 0.1185349192520667E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.5250918173022379E-01; v = 0.1913408643425751E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.8993480082038376E-01; v = 0.2452886577209897E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.1306023924436019; v = 0.2862408183288702E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.1732060388531418; v = 0.3178032258257357E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.2168727084820249; v = 0.3422945667633690E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.2609528309173586; v = 0.3612790520235922E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.3049252927938952; v = 0.3758638229818521E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.3483484138084404; v = 0.3868711798859953E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.3908321549106406; v = 0.3949429933189938E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.4320210071894814; v = 0.4006068107541156E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.4715824795890053; v = 0.4043192149672723E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.5091984794078453; v = 0.4064947495808078E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.5445580145650803; v = 0.4075245619813152E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6072575796841768; v = 0.4076423540893566E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6339484505755803; v = 0.4074280862251555E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6570718257486958; v = 0.4074163756012244E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6762557330090709; v = 0.4077647795071246E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6911161696923790; v = 0.4084517552782530E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.7012841911659961; v = 0.4092468459224052E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.7064559272410020; v = 0.4097872687240906E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6123554989894765E-01; v = 0.1738986811745028E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.1533070348312393; v = 0.2659616045280191E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.2563902605244206; v = 0.3240596008171533E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.3629346991663361; v = 0.3621195964432943E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.4683949968987538; v = 0.3868838330760539E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.5694479240657952; v = 0.4018911532693111E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.6634465430993955; v = 0.4089929432983252E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.1033958573552305; b = 0.3034544009063584E-01; v = 0.2279907527706409E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.1473521412414395; b = 0.6618803044247135E-01; v = 0.2715205490578897E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.1924552158705967; b = 0.1054431128987715; v = 0.3057917896703976E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.2381094362890328; b = 0.1468263551238858; v = 0.3326913052452555E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.2838121707936760; b = 0.1894486108187886; v = 0.3537334711890037E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3291323133373415; b = 0.2326374238761579; v = 0.3700567500783129E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3736896978741460; b = 0.2758485808485768; v = 0.3825245372589122E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4171406040760013; b = 0.3186179331996921; v = 0.3918125171518296E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4591677985256915; b = 0.3605329796303794; v = 0.3984720419937579E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4994733831718418; b = 0.4012147253586509; v = 0.4029746003338211E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5377731830445096; b = 0.4403050025570692; v = 0.4057428632156627E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5737917830001331; b = 0.4774565904277483; v = 0.4071719274114857E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.2027323586271389; b = 0.3544122504976147E-01; v = 0.2990236950664119E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.2516942375187273; b = 0.7418304388646328E-01; v = 0.3262951734212878E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3000227995257181; b = 0.1150502745727186; v = 0.3482634608242413E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3474806691046342; b = 0.1571963371209364; v = 0.3656596681700892E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3938103180359209; b = 0.1999631877247100; v = 0.3791740467794218E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4387519590455703; b = 0.2428073457846535; v = 0.3894034450156905E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4820503960077787; b = 0.2852575132906155; v = 0.3968600245508371E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5234573778475101; b = 0.3268884208674639; v = 0.4019931351420050E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5627318647235282; b = 0.3673033321675939; v = 0.4052108801278599E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5996390607156954; b = 0.4061211551830290; v = 0.4068978613940934E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3084780753791947; b = 0.3860125523100059E-01; v = 0.3454275351319704E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3589988275920223; b = 0.7928938987104867E-01; v = 0.3629963537007920E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4078628415881973; b = 0.1212614643030087; v = 0.3770187233889873E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4549287258889735; b = 0.1638770827382693; v = 0.3878608613694378E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5000278512957279; b = 0.2065965798260176; v = 0.3959065270221274E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5429785044928199; b = 0.2489436378852235; v = 0.4015286975463570E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5835939850491711; b = 0.2904811368946891; v = 0.4050866785614717E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6216870353444856; b = 0.3307941957666609; v = 0.4069320185051913E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4151104662709091; b = 0.4064829146052554E-01; v = 0.3760120964062763E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4649804275009218; b = 0.8258424547294755E-01; v = 0.3870969564418064E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5124695757009662; b = 0.1251841962027289; v = 0.3955287790534055E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5574711100606224; b = 0.1679107505976331; v = 0.4015361911302668E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5998597333287227; b = 0.2102805057358715; v = 0.4053836986719548E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6395007148516600; b = 0.2518418087774107; v = 0.4073578673299117E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5188456224746252; b = 0.4194321676077518E-01; v = 0.3954628379231406E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5664190707942778; b = 0.8457661551921499E-01; v = 0.4017645508847530E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6110464353283153; b = 0.1273652932519396; v = 0.4059030348651293E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6526430302051563; b = 0.1698173239076354; v = 0.4080565809484880E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6167551880377548; b = 0.4266398851548864E-01; v = 0.4063018753664651E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6607195418355383; b = 0.8551925814238349E-01; v = 0.4087191292799671E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); return end