function [ x, y, z, w ] = ld2354 ( ) %*****************************************************************************80; % %% LD2354 computes the 2354 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(2354,1); y = zeros(2354,1); z = zeros(2354,1); w = zeros(2354,1); a = 0.0; b = 0.0; v = 0.3922616270665292E-04; [ n, x, y, z, w ] = gen_oh ( 1, n, a, b, v, x, y, z, w ); v = 0.4703831750854424E-03; [ n, x, y, z, w ] = gen_oh ( 2, n, a, b, v, x, y, z, w ); v = 0.4678202801282136E-03; [ n, x, y, z, w ] = gen_oh ( 3, n, a, b, v, x, y, z, w ); a = 0.2290024646530589E-01; v = 0.1437832228979900E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.5779086652271284E-01; v = 0.2303572493577644E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.9863103576375984E-01; v = 0.2933110752447454E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.1428155792982185; v = 0.3402905998359838E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.1888978116601463; v = 0.3759138466870372E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.2359091682970210; v = 0.4030638447899798E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.2831228833706171; v = 0.4236591432242211E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.3299495857966693; v = 0.4390522656946746E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.3758840802660796; v = 0.4502523466626247E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.4204751831009480; v = 0.4580577727783541E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.4633068518751051; v = 0.4631391616615899E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.5039849474507313; v = 0.4660928953698676E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.5421265793440747; v = 0.4674751807936953E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6092660230557310; v = 0.4676414903932920E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6374654204984869; v = 0.4674086492347870E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6615136472609892; v = 0.4674928539483207E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6809487285958127; v = 0.4680748979686447E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6952980021665196; v = 0.4690449806389040E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.7041245497695400; v = 0.4699877075860818E-03; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.6744033088306065E-01; v = 0.2099942281069176E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.1678684485334166; v = 0.3172269150712804E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.2793559049539613; v = 0.3832051358546523E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.3935264218057639; v = 0.4252193818146985E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.5052629268232558; v = 0.4513807963755000E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.6107905315437531; v = 0.4657797469114178E-03; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); a = 0.1135081039843524; b = 0.3331954884662588E-01; v = 0.2733362800522836E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.1612866626099378; b = 0.7247167465436538E-01; v = 0.3235485368463559E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.2100786550168205; b = 0.1151539110849745; v = 0.3624908726013453E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.2592282009459942; b = 0.1599491097143677; v = 0.3925540070712828E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3081740561320203; b = 0.2058699956028027; v = 0.4156129781116235E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3564289781578164; b = 0.2521624953502911; v = 0.4330644984623263E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4035587288240703; b = 0.2982090785797674; v = 0.4459677725921312E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4491671196373903; b = 0.3434762087235733; v = 0.4551593004456795E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4928854782917489; b = 0.3874831357203437; v = 0.4613341462749918E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5343646791958988; b = 0.4297814821746926; v = 0.4651019618269806E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5732683216530990; b = 0.4699402260943537; v = 0.4670249536100625E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.2214131583218986; b = 0.3873602040643895E-01; v = 0.3549555576441708E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.2741796504750071; b = 0.8089496256902013E-01; v = 0.3856108245249010E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3259797439149485; b = 0.1251732177620872; v = 0.4098622845756882E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3765441148826891; b = 0.1706260286403185; v = 0.4286328604268950E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4255773574530558; b = 0.2165115147300408; v = 0.4427802198993945E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4727795117058430; b = 0.2622089812225259; v = 0.4530473511488561E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5178546895819012; b = 0.3071721431296201; v = 0.4600805475703138E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5605141192097460; b = 0.3508998998801138; v = 0.4644599059958017E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6004763319352512; b = 0.3929160876166931; v = 0.4667274455712508E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3352842634946949; b = 0.4202563457288019E-01; v = 0.4069360518020356E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.3891971629814670; b = 0.8614309758870850E-01; v = 0.4260442819919195E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4409875565542281; b = 0.1314500879380001; v = 0.4408678508029063E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4904893058592484; b = 0.1772189657383859; v = 0.4518748115548597E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5375056138769549; b = 0.2228277110050294; v = 0.4595564875375116E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5818255708669969; b = 0.2677179935014386; v = 0.4643988774315846E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6232334858144959; b = 0.3113675035544165; v = 0.4668827491646946E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.4489485354492058; b = 0.4409162378368174E-01; v = 0.4400541823741973E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5015136875933150; b = 0.8939009917748489E-01; v = 0.4514512890193797E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5511300550512623; b = 0.1351806029383365; v = 0.4596198627347549E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5976720409858000; b = 0.1808370355053196; v = 0.4648659016801781E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6409956378989354; b = 0.2257852192301602; v = 0.4675502017157673E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.5581222330827514; b = 0.4532173421637160E-01; v = 0.4598494476455523E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6074705984161695; b = 0.9117488031840314E-01; v = 0.4654916955152048E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6532272537379033; b = 0.1369294213140155; v = 0.4684709779505137E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); a = 0.6594761494500487; b = 0.4589901487275583E-01; v = 0.4691445539106986E-03; [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w ); return end