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

%*****************************************************************************80
%
%% ld1202() computes the 1202 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.
%
%  Output:
%
%    real X(N), Y(N), Z(N), W(N), the coordinates
%    and weights of the points.
%
  n = 0;
  x = zeros(1202,1);
  y = zeros(1202,1);
  z = zeros(1202,1);
  w = zeros(1202,1);
  a = 0.0;
  b = 0.0;
  v = 0.1105189233267572E-03;
  [ n, x, y, z, w ] = gen_oh ( 1, n, a, b, v, x, y, z, w );
  v = 0.9205232738090741E-03;
  [ n, x, y, z, w ] = gen_oh ( 2, n, a, b, v, x, y, z, w );
  v = 0.9133159786443561E-03;
  [ n, x, y, z, w ] = gen_oh ( 3, n, a, b, v, x, y, z, w );
  a = 0.3712636449657089E-01;
  v = 0.3690421898017899E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.9140060412262223E-01;
  v = 0.5603990928680660E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1531077852469906;
  v = 0.6865297629282609E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2180928891660612;
  v = 0.7720338551145630E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2839874532200175;
  v = 0.8301545958894795E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3491177600963764;
  v = 0.8686692550179628E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4121431461444309;
  v = 0.8927076285846890E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4718993627149127;
  v = 0.9060820238568219E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.5273145452842337;
  v = 0.9119777254940867E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6209475332444019;
  v = 0.9128720138604181E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6569722711857291;
  v = 0.9130714935691735E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6841788309070143;
  v = 0.9152873784554116E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.7012604330123631;
  v = 0.9187436274321654E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1072382215478166;
  v = 0.5176977312965694E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.2582068959496968;
  v = 0.7331143682101417E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.4172752955306717;
  v = 0.8463232836379928E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.5700366911792503;
  v = 0.9031122694253992E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.9827986018263947;
  b = 0.1771774022615325;
  v = 0.6485778453163257E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.9624249230326228;
  b = 0.2475716463426288;
  v = 0.7435030910982369E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.9402007994128811;
  b = 0.3354616289066489;
  v = 0.7998527891839054E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.9320822040143202;
  b = 0.3173615246611977;
  v = 0.8101731497468018E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.9043674199393299;
  b = 0.4090268427085357;
  v = 0.8483389574594331E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.8912407560074747;
  b = 0.3854291150669224;
  v = 0.8556299257311812E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.8676435628462708;
  b = 0.4932221184851285;
  v = 0.8803208679738260E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.8581979986041619;
  b = 0.4785320675922435;
  v = 0.8811048182425720E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.8396753624049856;
  b = 0.4507422593157064;
  v = 0.8850282341265444E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.8165288564022188;
  b = 0.5632123020762100;
  v = 0.9021342299040653E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.8015469370783529;
  b = 0.5434303569693900;
  v = 0.9010091677105086E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.7773563069070351;
  b = 0.5123518486419871;
  v = 0.9022692938426915E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.7661621213900394;
  b = 0.6394279634749102;
  v = 0.9158016174693465E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.7553584143533510;
  b = 0.6269805509024392;
  v = 0.9131578003189435E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.7344305757559503;
  b = 0.6031161693096310;
  v = 0.9107813579482705E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.7043837184021765;
  b = 0.5693702498468441;
  v = 0.9105760258970126E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  
  return
end
