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

%*****************************************************************************80
%
%% LD2030 computes the 2030 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(2030,1);
  y = zeros(2030,1);
  z = zeros(2030,1);
  w = zeros(2030,1);
  a = 0.0;
  b = 0.0;
  v = 0.4656031899197431E-04;
  [ n, x, y, z, w ] = gen_oh ( 1, n, a, b, v, x, y, z, w );
  v = 0.5421549195295507E-03;
  [ n, x, y, z, w ] = gen_oh ( 3, n, a, b, v, x, y, z, w );
  a = 0.2540835336814348E-01;
  v = 0.1778522133346553E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6399322800504915E-01;
  v = 0.2811325405682796E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1088269469804125;
  v = 0.3548896312631459E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1570670798818287;
  v = 0.4090310897173364E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2071163932282514;
  v = 0.4493286134169965E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2578914044450844;
  v = 0.4793728447962723E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3085687558169623;
  v = 0.5015415319164265E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3584719706267024;
  v = 0.5175127372677937E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4070135594428709;
  v = 0.5285522262081019E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4536618626222638;
  v = 0.5356832703713962E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4979195686463577;
  v = 0.5397914736175170E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.5393075111126999;
  v = 0.5416899441599930E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6115617676843916;
  v = 0.5419308476889938E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6414308435160159;
  v = 0.5416936902030596E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6664099412721607;
  v = 0.5419544338703164E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6859161771214913;
  v = 0.5428983656630975E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6993625593503890;
  v = 0.5442286500098193E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.7062393387719380;
  v = 0.5452250345057301E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.7479028168349763E-01;
  v = 0.2568002497728530E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.1848951153969366;
  v = 0.3827211700292145E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.3059529066581305;
  v = 0.4579491561917824E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.4285556101021362;
  v = 0.5042003969083574E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.5468758653496526;
  v = 0.5312708889976025E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.6565821978343439;
  v = 0.5438401790747117E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.1253901572367117;
  b = 0.3681917226439641E-01;
  v = 0.3316041873197344E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1775721510383941;
  b = 0.7982487607213301E-01;
  v = 0.3899113567153771E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2305693358216114;
  b = 0.1264640966592335;
  v = 0.4343343327201309E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2836502845992063;
  b = 0.1751585683418957;
  v = 0.4679415262318919E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3361794746232590;
  b = 0.2247995907632670;
  v = 0.4930847981631031E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3875979172264824;
  b = 0.2745299257422246;
  v = 0.5115031867540091E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4374019316999074;
  b = 0.3236373482441118;
  v = 0.5245217148457367E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4851275843340022;
  b = 0.3714967859436741;
  v = 0.5332041499895321E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5303391803806868;
  b = 0.4175353646321745;
  v = 0.5384583126021542E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5726197380596287;
  b = 0.4612084406355461;
  v = 0.5411067210798852E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2431520732564863;
  b = 0.4258040133043952E-01;
  v = 0.4259797391468714E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3002096800895869;
  b = 0.8869424306722721E-01;
  v = 0.4604931368460021E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3558554457457432;
  b = 0.1368811706510655;
  v = 0.4871814878255202E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4097782537048887;
  b = 0.1860739985015033;
  v = 0.5072242910074885E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4616337666067458;
  b = 0.2354235077395853;
  v = 0.5217069845235350E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5110707008417874;
  b = 0.2842074921347011;
  v = 0.5315785966280310E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5577415286163795;
  b = 0.3317784414984102;
  v = 0.5376833708758905E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6013060431366950;
  b = 0.3775299002040700;
  v = 0.5408032092069521E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3661596767261781;
  b = 0.4599367887164592E-01;
  v = 0.4842744917904866E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4237633153506581;
  b = 0.9404893773654421E-01;
  v = 0.5048926076188130E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4786328454658452;
  b = 0.1431377109091971;
  v = 0.5202607980478373E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5305702076789774;
  b = 0.1924186388843570;
  v = 0.5309932388325743E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5793436224231788;
  b = 0.2411590944775190;
  v = 0.5377419770895208E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6247069017094747;
  b = 0.2886871491583605;
  v = 0.5411696331677717E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4874315552535204;
  b = 0.4804978774953206E-01;
  v = 0.5197996293282420E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5427337322059053;
  b = 0.9716857199366665E-01;
  v = 0.5311120836622945E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5943493747246700;
  b = 0.1465205839795055;
  v = 0.5384309319956951E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6421314033564943;
  b = 0.1953579449803574;
  v = 0.5421859504051886E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6020628374713980;
  b = 0.4916375015738108E-01;
  v = 0.5390948355046314E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6529222529856881;
  b = 0.9861621540127005E-01;
  v = 0.5433312705027845E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  
  return
end
