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. % % Parameters: % % 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