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

%*****************************************************************************80
%
%% LD0974 computes the 974 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(974,1);
  y = zeros(974,1);
  z = zeros(974,1);
  w = zeros(974,1);
  a = 0.0;
  b = 0.0;
  v = 0.1438294190527431E-03;
  [ n, x, y, z, w ] = gen_oh ( 1, n, a, b, v, x, y, z, w );
  v = 0.1125772288287004E-02;
  [ n, x, y, z, w ] = gen_oh ( 3, n, a, b, v, x, y, z, w );
  a = 0.4292963545341347E-01;
  v = 0.4948029341949241E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1051426854086404;
  v = 0.7357990109125470E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1750024867623087;
  v = 0.8889132771304384E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2477653379650257;
  v = 0.9888347838921435E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3206567123955957;
  v = 0.1053299681709471E-02;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3916520749849983;
  v = 0.1092778807014578E-02;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4590825874187624;
  v = 0.1114389394063227E-02;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.5214563888415861;
  v = 0.1123724788051555E-02;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6253170244654199;
  v = 0.1125239325243814E-02;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6637926744523170;
  v = 0.1126153271815905E-02;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6910410398498301;
  v = 0.1130286931123841E-02;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.7052907007457760;
  v = 0.1134986534363955E-02;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1236686762657990;
  v = 0.6823367927109931E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.2940777114468387;
  v = 0.9454158160447096E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.4697753849207649;
  v = 0.1074429975385679E-02;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.6334563241139567;
  v = 0.1129300086569132E-02;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.5974048614181342E-01;
  b = 0.2029128752777523;
  v = 0.8436884500901954E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1375760408473636;
  b = 0.4602621942484054;
  v = 0.1075255720448885E-02;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3391016526336286;
  b = 0.5030673999662036;
  v = 0.1108577236864462E-02;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1271675191439820;
  b = 0.2817606422442134;
  v = 0.9566475323783357E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2693120740413512;
  b = 0.4331561291720157;
  v = 0.1080663250717391E-02;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1419786452601918;
  b = 0.6256167358580814;
  v = 0.1126797131196295E-02;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6709284600738255E-01;
  b = 0.3798395216859157;
  v = 0.1022568715358061E-02;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.7057738183256172E-01;
  b = 0.5517505421423520;
  v = 0.1108960267713108E-02;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2783888477882155;
  b = 0.6029619156159187;
  v = 0.1122790653435766E-02;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1979578938917407;
  b = 0.3589606329589096;
  v = 0.1032401847117460E-02;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2087307061103274;
  b = 0.5348666438135476;
  v = 0.1107249382283854E-02;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4055122137872836;
  b = 0.5674997546074373;
  v = 0.1121780048519972E-02;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  
  return
end
