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

%*****************************************************************************80;
%
%% LD3470 computes the 3470 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(3470,1);
  y = zeros(3470,1);
  z = zeros(3470,1);
  w = zeros(3470,1);
  a = 0.0;
  b = 0.0;
  v = 0.2040382730826330E-04;
  [ n, x, y, z, w ] = gen_oh ( 1, n, a, b, v, x, y, z, w );
  v = 0.3178149703889544E-03;
  [ n, x, y, z, w ] = gen_oh ( 3, n, a, b, v, x, y, z, w );
  a = 0.1721420832906233E-01;
  v = 0.8288115128076110E-04;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4408875374981770E-01;
  v = 0.1360883192522954E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.7594680813878681E-01;
  v = 0.1766854454542662E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1108335359204799;
  v = 0.2083153161230153E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1476517054388567;
  v = 0.2333279544657158E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1856731870860615;
  v = 0.2532809539930247E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2243634099428821;
  v = 0.2692472184211158E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2633006881662727;
  v = 0.2819949946811885E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3021340904916283;
  v = 0.2920953593973030E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3405594048030089;
  v = 0.2999889782948352E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3783044434007372;
  v = 0.3060292120496902E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4151194767407910;
  v = 0.3105109167522192E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4507705766443257;
  v = 0.3136902387550312E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4850346056573187;
  v = 0.3157984652454632E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.5176950817792470;
  v = 0.3170516518425422E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.5485384240820989;
  v = 0.3176568425633755E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6039117238943308;
  v = 0.3177198411207062E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6279956655573113;
  v = 0.3175519492394733E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6493636169568952;
  v = 0.3174654952634756E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6677644117704504;
  v = 0.3175676415467654E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6829368572115624;
  v = 0.3178923417835410E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6946195818184121;
  v = 0.3183788287531909E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.7025711542057026;
  v = 0.3188755151918807E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.7066004767140119;
  v = 0.3191916889313849E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.5132537689946062E-01;
  v = 0.1231779611744508E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.1297994661331225;
  v = 0.1924661373839880E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.2188852049401307;
  v = 0.2380881867403424E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.3123174824903457;
  v = 0.2693100663037885E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.4064037620738195;
  v = 0.2908673382834366E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.4984958396944782;
  v = 0.3053914619381535E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.5864975046021365;
  v = 0.3143916684147777E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.6686711634580175;
  v = 0.3187042244055363E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.8715738780835950E-01;
  b = 0.2557175233367578E-01;
  v = 0.1635219535869790E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1248383123134007;
  b = 0.5604823383376681E-01;
  v = 0.1968109917696070E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1638062693383378;
  b = 0.8968568601900765E-01;
  v = 0.2236754342249974E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2035586203373176;
  b = 0.1254086651976279;
  v = 0.2453186687017181E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2436798975293774;
  b = 0.1624780150162012;
  v = 0.2627551791580541E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2838207507773806;
  b = 0.2003422342683208;
  v = 0.2767654860152220E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3236787502217692;
  b = 0.2385628026255263;
  v = 0.2879467027765895E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3629849554840691;
  b = 0.2767731148783578;
  v = 0.2967639918918702E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4014948081992087;
  b = 0.3146542308245309;
  v = 0.3035900684660351E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4389818379260225;
  b = 0.3519196415895088;
  v = 0.3087338237298308E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4752331143674377;
  b = 0.3883050984023654;
  v = 0.3124608838860167E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5100457318374018;
  b = 0.4235613423908649;
  v = 0.3150084294226743E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5432238388954868;
  b = 0.4574484717196220;
  v = 0.3165958398598402E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5745758685072442;
  b = 0.4897311639255524;
  v = 0.3174320440957372E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1723981437592809;
  b = 0.3010630597881105E-01;
  v = 0.2182188909812599E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2149553257844597;
  b = 0.6326031554204694E-01;
  v = 0.2399727933921445E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2573256081247422;
  b = 0.9848566980258631E-01;
  v = 0.2579796133514652E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2993163751238106;
  b = 0.1350835952384266;
  v = 0.2727114052623535E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3407238005148000;
  b = 0.1725184055442181;
  v = 0.2846327656281355E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3813454978483264;
  b = 0.2103559279730725;
  v = 0.2941491102051334E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4209848104423343;
  b = 0.2482278774554860;
  v = 0.3016049492136107E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4594519699996300;
  b = 0.2858099509982883;
  v = 0.3072949726175648E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4965640166185930;
  b = 0.3228075659915428;
  v = 0.3114768142886460E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5321441655571562;
  b = 0.3589459907204151;
  v = 0.3143823673666223E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5660208438582166;
  b = 0.3939630088864310;
  v = 0.3162269764661535E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5980264315964364;
  b = 0.4276029922949089;
  v = 0.3172164663759821E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2644215852350733;
  b = 0.3300939429072552E-01;
  v = 0.2554575398967435E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3090113743443063;
  b = 0.6803887650078501E-01;
  v = 0.2701704069135677E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3525871079197808;
  b = 0.1044326136206709;
  v = 0.2823693413468940E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3950418005354029;
  b = 0.1416751597517679;
  v = 0.2922898463214289E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4362475663430163;
  b = 0.1793408610504821;
  v = 0.3001829062162428E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4760661812145854;
  b = 0.2170630750175722;
  v = 0.3062890864542953E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5143551042512103;
  b = 0.2545145157815807;
  v = 0.3108328279264746E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5509709026935597;
  b = 0.2913940101706601;
  v = 0.3140243146201245E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5857711030329428;
  b = 0.3274169910910705;
  v = 0.3160638030977130E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6186149917404392;
  b = 0.3623081329317265;
  v = 0.3171462882206275E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3586894569557064;
  b = 0.3497354386450040E-01;
  v = 0.2812388416031796E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4035266610019441;
  b = 0.7129736739757095E-01;
  v = 0.2912137500288045E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4467775312332510;
  b = 0.1084758620193165;
  v = 0.2993241256502206E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4883638346608543;
  b = 0.1460915689241772;
  v = 0.3057101738983822E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5281908348434601;
  b = 0.1837790832369980;
  v = 0.3105319326251432E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5661542687149311;
  b = 0.2212075390874021;
  v = 0.3139565514428167E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6021450102031452;
  b = 0.2580682841160985;
  v = 0.3161543006806366E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6360520783610050;
  b = 0.2940656362094121;
  v = 0.3172985960613294E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4521611065087196;
  b = 0.3631055365867002E-01;
  v = 0.2989400336901431E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4959365651560963;
  b = 0.7348318468484350E-01;
  v = 0.3054555883947677E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5376815804038283;
  b = 0.1111087643812648;
  v = 0.3104764960807702E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5773314480243768;
  b = 0.1488226085145408;
  v = 0.3141015825977616E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6148113245575056;
  b = 0.1862892274135151;
  v = 0.3164520621159896E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6500407462842380;
  b = 0.2231909701714456;
  v = 0.3176652305912204E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5425151448707213;
  b = 0.3718201306118944E-01;
  v = 0.3105097161023939E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5841860556907931;
  b = 0.7483616335067346E-01;
  v = 0.3143014117890550E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6234632186851500;
  b = 0.1125990834266120;
  v = 0.3168172866287200E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6602934551848843;
  b = 0.1501303813157619;
  v = 0.3181401865570968E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6278573968375105;
  b = 0.3767559930245720E-01;
  v = 0.3170663659156037E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6665611711264577;
  b = 0.7548443301360158E-01;
  v = 0.3185447944625510E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );

  return
end
