function [ x, y, z, w ] = ld3074 ( )

%*****************************************************************************80;
%
%% LD3074 computes the 3074 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(3074,1);
  y = zeros(3074,1);
  z = zeros(3074,1);
  w = zeros(3074,1);
  a = 0.0;
  b = 0.0;
  v = 0.2599095953754734E-04;
  [ n, x, y, z, w ] = gen_oh ( 1, n, a, b, v, x, y, z, w );
  v = 0.3603134089687541E-03;
  [ n, x, y, z, w ] = gen_oh ( 2, n, a, b, v, x, y, z, w );
  v = 0.3586067974412447E-03;
  [ n, x, y, z, w ] = gen_oh ( 3, n, a, b, v, x, y, z, w );
  a = 0.1886108518723392E-01;
  v = 0.9831528474385880E-04;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4800217244625303E-01;
  v = 0.1605023107954450E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.8244922058397242E-01;
  v = 0.2072200131464099E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1200408362484023;
  v = 0.2431297618814187E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1595773530809965;
  v = 0.2711819064496707E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2002635973434064;
  v = 0.2932762038321116E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2415127590139982;
  v = 0.3107032514197368E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2828584158458477;
  v = 0.3243808058921213E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3239091015338138;
  v = 0.3349899091374030E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3643225097962194;
  v = 0.3430580688505218E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4037897083691802;
  v = 0.3490124109290343E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4420247515194127;
  v = 0.3532148948561955E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4787572538464938;
  v = 0.3559862669062833E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.5137265251275234;
  v = 0.3576224317551411E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.5466764056654611;
  v = 0.3584050533086076E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6054859420813535;
  v = 0.3584903581373224E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6308106701764562;
  v = 0.3582991879040586E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6530369230179584;
  v = 0.3582371187963125E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6718609524611158;
  v = 0.3584353631122350E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6869676499894013;
  v = 0.3589120166517785E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6980467077240748;
  v = 0.3595445704531601E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.7048241721250522;
  v = 0.3600943557111074E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.5591105222058232E-01;
  v = 0.1456447096742039E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.1407384078513916;
  v = 0.2252370188283782E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.2364035438976309;
  v = 0.2766135443474897E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.3360602737818170;
  v = 0.3110729491500851E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.4356292630054665;
  v = 0.3342506712303391E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.5321569415256174;
  v = 0.3491981834026860E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.6232956305040554;
  v = 0.3576003604348932E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.9469870086838469E-01;
  b = 0.2778748387309470E-01;
  v = 0.1921921305788564E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1353170300568141;
  b = 0.6076569878628364E-01;
  v = 0.2301458216495632E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1771679481726077;
  b = 0.9703072762711040E-01;
  v = 0.2604248549522893E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2197066664231751;
  b = 0.1354112458524762;
  v = 0.2845275425870697E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2624783557374927;
  b = 0.1750996479744100;
  v = 0.3036870897974840E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3050969521214442;
  b = 0.2154896907449802;
  v = 0.3188414832298066E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3472252637196021;
  b = 0.2560954625740152;
  v = 0.3307046414722089E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3885610219026360;
  b = 0.2965070050624096;
  v = 0.3398330969031360E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4288273776062765;
  b = 0.3363641488734497;
  v = 0.3466757899705373E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4677662471302948;
  b = 0.3753400029836788;
  v = 0.3516095923230054E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5051333589553359;
  b = 0.4131297522144286;
  v = 0.3549645184048486E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5406942145810492;
  b = 0.4494423776081795;
  v = 0.3570415969441392E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5742204122576457;
  b = 0.4839938958841502;
  v = 0.3581251798496118E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1865407027225188;
  b = 0.3259144851070796E-01;
  v = 0.2543491329913348E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2321186453689432;
  b = 0.6835679505297343E-01;
  v = 0.2786711051330776E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2773159142523882;
  b = 0.1062284864451989;
  v = 0.2985552361083679E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3219200192237254;
  b = 0.1454404409323047;
  v = 0.3145867929154039E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3657032593944029;
  b = 0.1854018282582510;
  v = 0.3273290662067609E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4084376778363622;
  b = 0.2256297412014750;
  v = 0.3372705511943501E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4499004945751427;
  b = 0.2657104425000896;
  v = 0.3448274437851510E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4898758141326335;
  b = 0.3052755487631557;
  v = 0.3503592783048583E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5281547442266309;
  b = 0.3439863920645423;
  v = 0.3541854792663162E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5645346989813992;
  b = 0.3815229456121914;
  v = 0.3565995517909428E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5988181252159848;
  b = 0.4175752420966734;
  v = 0.3578802078302898E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2850425424471603;
  b = 0.3562149509862536E-01;
  v = 0.2958644592860982E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3324619433027876;
  b = 0.7330318886871096E-01;
  v = 0.3119548129116835E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3785848333076282;
  b = 0.1123226296008472;
  v = 0.3250745225005984E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4232891028562115;
  b = 0.1521084193337708;
  v = 0.3355153415935208E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4664287050829722;
  b = 0.1921844459223610;
  v = 0.3435847568549328E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5078458493735726;
  b = 0.2321360989678303;
  v = 0.3495786831622488E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5473779816204180;
  b = 0.2715886486360520;
  v = 0.3537767805534621E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5848617133811376;
  b = 0.3101924707571355;
  v = 0.3564459815421428E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6201348281584888;
  b = 0.3476121052890973;
  v = 0.3578464061225468E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3852191185387871;
  b = 0.3763224880035108E-01;
  v = 0.3239748762836212E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4325025061073423;
  b = 0.7659581935637135E-01;
  v = 0.3345491784174287E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4778486229734490;
  b = 0.1163381306083900;
  v = 0.3429126177301782E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5211663693009000;
  b = 0.1563890598752899;
  v = 0.3492420343097421E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5623469504853703;
  b = 0.1963320810149200;
  v = 0.3537399050235257E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6012718188659246;
  b = 0.2357847407258738;
  v = 0.3566209152659172E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6378179206390117;
  b = 0.2743846121244060;
  v = 0.3581084321919782E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4836936460214534;
  b = 0.3895902610739024E-01;
  v = 0.3426522117591512E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5293792562683797;
  b = 0.7871246819312640E-01;
  v = 0.3491848770121379E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5726281253100033;
  b = 0.1187963808202981;
  v = 0.3539318235231476E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6133658776169068;
  b = 0.1587914708061787;
  v = 0.3570231438458694E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6515085491865307;
  b = 0.1983058575227646;
  v = 0.3586207335051714E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5778692716064976;
  b = 0.3977209689791542E-01;
  v = 0.3541196205164025E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6207904288086192;
  b = 0.7990157592981152E-01;
  v = 0.3574296911573953E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6608688171046802;
  b = 0.1199671308754309;
  v = 0.3591993279818963E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6656263089489130;
  b = 0.4015955957805969E-01;
  v = 0.3595855034661997E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  
  return
end
