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

%*****************************************************************************80;
%
%% LD4802 computes the 4802 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(4802,1);
  y = zeros(4802,1);
  z = zeros(4802,1);
  w = zeros(4802,1);
  a = 0.0;
  b = 0.0;
  v = 0.9687521879420705E-04;
  [ n, x, y, z, w ] = gen_oh ( 1, n, a, b, v, x, y, z, w );
  v = 0.2307897895367918E-03;
  [ n, x, y, z, w ] = gen_oh ( 2, n, a, b, v, x, y, z, w );
  v = 0.2297310852498558E-03;
  [ n, x, y, z, w ] = gen_oh ( 3, n, a, b, v, x, y, z, w );
  a = 0.2335728608887064E-01;
  v = 0.7386265944001919E-04;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4352987836550653E-01;
  v = 0.8257977698542210E-04;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6439200521088801E-01;
  v = 0.9706044762057630E-04;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.9003943631993181E-01;
  v = 0.1302393847117003E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1196706615548473;
  v = 0.1541957004600968E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1511715412838134;
  v = 0.1704459770092199E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.1835982828503801;
  v = 0.1827374890942906E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2165081259155405;
  v = 0.1926360817436107E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2496208720417563;
  v = 0.2008010239494833E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.2827200673567900;
  v = 0.2075635983209175E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3156190823994346;
  v = 0.2131306638690909E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3481476793749115;
  v = 0.2176562329937335E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.3801466086947226;
  v = 0.2212682262991018E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4114652119634011;
  v = 0.2240799515668565E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4419598786519751;
  v = 0.2261959816187525E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4714925949329543;
  v = 0.2277156368808855E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4999293972879466;
  v = 0.2287351772128336E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.5271387221431248;
  v = 0.2293490814084085E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.5529896780837761;
  v = 0.2296505312376273E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6000856099481712;
  v = 0.2296793832318756E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6210562192785175;
  v = 0.2295785443842974E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6401165879934240;
  v = 0.2295017931529102E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6571144029244334;
  v = 0.2295059638184868E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6718910821718863;
  v = 0.2296232343237362E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6842845591099010;
  v = 0.2298530178740771E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.6941353476269816;
  v = 0.2301579790280501E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.7012965242212991;
  v = 0.2304690404996513E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.7056471428242644;
  v = 0.2307027995907102E-03;
  [ n, x, y, z, w ] = gen_oh ( 4, n, a, b, v, x, y, z, w );
  a = 0.4595557643585895E-01;
  v = 0.9312274696671092E-04;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.1049316742435023;
  v = 0.1199919385876926E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.1773548879549274;
  v = 0.1598039138877690E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.2559071411236127;
  v = 0.1822253763574900E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.3358156837985898;
  v = 0.1988579593655040E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.4155835743763893;
  v = 0.2112620102533307E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.4937894296167472;
  v = 0.2201594887699007E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.5691569694793316;
  v = 0.2261622590895036E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.6405840854894251;
  v = 0.2296458453435705E-03;
  [ n, x, y, z, w ] = gen_oh ( 5, n, a, b, v, x, y, z, w );
  a = 0.7345133894143348E-01;
  b = 0.2177844081486067E-01;
  v = 0.1006006990267000E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1009859834044931;
  b = 0.4590362185775188E-01;
  v = 0.1227676689635876E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1324289619748758;
  b = 0.7255063095690877E-01;
  v = 0.1467864280270117E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1654272109607127;
  b = 0.1017825451960684;
  v = 0.1644178912101232E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1990767186776461;
  b = 0.1325652320980364;
  v = 0.1777664890718961E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2330125945523278;
  b = 0.1642765374496765;
  v = 0.1884825664516690E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2670080611108287;
  b = 0.1965360374337889;
  v = 0.1973269246453848E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3008753376294316;
  b = 0.2290726770542238;
  v = 0.2046767775855328E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3344475596167860;
  b = 0.2616645495370823;
  v = 0.2107600125918040E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3675709724070786;
  b = 0.2941150728843141;
  v = 0.2157416362266829E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4001000887587812;
  b = 0.3262440400919066;
  v = 0.2197557816920721E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4318956350436028;
  b = 0.3578835350611916;
  v = 0.2229192611835437E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4628239056795531;
  b = 0.3888751854043678;
  v = 0.2253385110212775E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4927563229773636;
  b = 0.4190678003222840;
  v = 0.2271137107548774E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5215687136707969;
  b = 0.4483151836883852;
  v = 0.2283414092917525E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5491402346984905;
  b = 0.4764740676087880;
  v = 0.2291161673130077E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5753520160126075;
  b = 0.5034021310998277;
  v = 0.2295313908576598E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1388326356417754;
  b = 0.2435436510372806E-01;
  v = 0.1438204721359031E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.1743686900537244;
  b = 0.5118897057342652E-01;
  v = 0.1607738025495257E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2099737037950268;
  b = 0.8014695048539634E-01;
  v = 0.1741483853528379E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2454492590908548;
  b = 0.1105117874155699;
  v = 0.1851918467519151E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2807219257864278;
  b = 0.1417950531570966;
  v = 0.1944628638070613E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3156842271975842;
  b = 0.1736604945719597;
  v = 0.2022495446275152E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3502090945177752;
  b = 0.2058466324693981;
  v = 0.2087462382438514E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3841684849519686;
  b = 0.2381284261195919;
  v = 0.2141074754818308E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4174372367906016;
  b = 0.2703031270422569;
  v = 0.2184640913748162E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4498926465011892;
  b = 0.3021845683091309;
  v = 0.2219309165220329E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4814146229807701;
  b = 0.3335993355165720;
  v = 0.2246123118340624E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5118863625734701;
  b = 0.3643833735518232;
  v = 0.2266062766915125E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5411947455119144;
  b = 0.3943789541958179;
  v = 0.2280072952230796E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5692301500357246;
  b = 0.4234320144403542;
  v = 0.2289082025202583E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5958857204139576;
  b = 0.4513897947419260;
  v = 0.2294012695120025E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2156270284785766;
  b = 0.2681225755444491E-01;
  v = 0.1722434488736947E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2532385054909710;
  b = 0.5557495747805614E-01;
  v = 0.1830237421455091E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2902564617771537;
  b = 0.8569368062950249E-01;
  v = 0.1923855349997633E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3266979823143256;
  b = 0.1167367450324135;
  v = 0.2004067861936271E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3625039627493614;
  b = 0.1483861994003304;
  v = 0.2071817297354263E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3975838937548699;
  b = 0.1803821503011405;
  v = 0.2128250834102103E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4318396099009774;
  b = 0.2124962965666424;
  v = 0.2174513719440102E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4651706555732742;
  b = 0.2445221837805913;
  v = 0.2211661839150214E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4974752649620969;
  b = 0.2762701224322987;
  v = 0.2240665257813102E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5286517579627517;
  b = 0.3075627775211328;
  v = 0.2262439516632620E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5586001195731895;
  b = 0.3382311089826877;
  v = 0.2277874557231869E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5872229902021319;
  b = 0.3681108834741399;
  v = 0.2287854314454994E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6144258616235123;
  b = 0.3970397446872839;
  v = 0.2293268499615575E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.2951676508064861;
  b = 0.2867499538750441E-01;
  v = 0.1912628201529828E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3335085485472725;
  b = 0.5867879341903510E-01;
  v = 0.1992499672238701E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3709561760636381;
  b = 0.8961099205022284E-01;
  v = 0.2061275533454027E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4074722861667498;
  b = 0.1211627927626297;
  v = 0.2119318215968572E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4429923648839117;
  b = 0.1530748903554898;
  v = 0.2167416581882652E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4774428052721736;
  b = 0.1851176436721877;
  v = 0.2206430730516600E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5107446539535904;
  b = 0.2170829107658179;
  v = 0.2237186938699523E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5428151370542935;
  b = 0.2487786689026271;
  v = 0.2260480075032884E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5735699292556964;
  b = 0.2800239952795016;
  v = 0.2277098884558542E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6029253794562866;
  b = 0.3106445702878119;
  v = 0.2287845715109671E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6307998987073145;
  b = 0.3404689500841194;
  v = 0.2293547268236294E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.3752652273692719;
  b = 0.2997145098184479E-01;
  v = 0.2056073839852528E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4135383879344028;
  b = 0.6086725898678011E-01;
  v = 0.2114235865831876E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4506113885153907;
  b = 0.9238849548435643E-01;
  v = 0.2163175629770551E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4864401554606072;
  b = 0.1242786603851851;
  v = 0.2203392158111650E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5209708076611709;
  b = 0.1563086731483386;
  v = 0.2235473176847839E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5541422135830122;
  b = 0.1882696509388506;
  v = 0.2260024141501235E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5858880915113817;
  b = 0.2199672979126059;
  v = 0.2277675929329182E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6161399390603444;
  b = 0.2512165482924867;
  v = 0.2289102112284834E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6448296482255090;
  b = 0.2818368701871888;
  v = 0.2295027954625118E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4544796274917948;
  b = 0.3088970405060312E-01;
  v = 0.2161281589879992E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.4919389072146628;
  b = 0.6240947677636835E-01;
  v = 0.2201980477395102E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5279313026985183;
  b = 0.9430706144280313E-01;
  v = 0.2234952066593166E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5624169925571135;
  b = 0.1263547818770374;
  v = 0.2260540098520838E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5953484627093287;
  b = 0.1583430788822594;
  v = 0.2279157981899988E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6266730715339185;
  b = 0.1900748462555988;
  v = 0.2291296918565571E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6563363204278871;
  b = 0.2213599519592567;
  v = 0.2297533752536649E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5314574716585696;
  b = 0.3152508811515374E-01;
  v = 0.2234927356465995E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.5674614932298185;
  b = 0.6343865291465561E-01;
  v = 0.2261288012985219E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6017706004970264;
  b = 0.9551503504223951E-01;
  v = 0.2280818160923688E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6343471270264178;
  b = 0.1275440099801196;
  v = 0.2293773295180159E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6651494599127802;
  b = 0.1593252037671960;
  v = 0.2300528767338634E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6050184986005704;
  b = 0.3192538338496105E-01;
  v = 0.2281893855065666E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6390163550880400;
  b = 0.6402824353962306E-01;
  v = 0.2295720444840727E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6711199107088448;
  b = 0.9609805077002909E-01;
  v = 0.2303227649026753E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  a = 0.6741354429572275;
  b = 0.3211853196273233E-01;
  v = 0.2304831913227114E-03;
  [ n, x, y, z, w ] = gen_oh ( 6, n, a, b, v, x, y, z, w );
  
  return
end
