function [ x, y, z, w ] = ld0074 ( ) %*****************************************************************************80 % %% LD0074 computes the 74 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(74,1); y = zeros(74,1); z = zeros(74,1); w = zeros(74,1); a = 0.0; b = 0.0; v = 0.5130671797338464E-03; [ n, x, y, z, w ] = gen_oh ( 1, n, a, b, v, x, y, z, w ); v = 0.1660406956574204E-01; [ n, x, y, z, w ] = gen_oh ( 2, n, a, b, v, x, y, z, w ); v = -0.2958603896103896E-01; [ n, x, y, z, w ] = gen_oh ( 3, n, a, b, v, x, y, z, w ); a = 0.4803844614152614; v = 0.2657620708215946E-01; [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w ); a = 0.3207726489807764; v = 0.1652217099371571E-01; [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w ); return end