00001
00002
00003
00004
00005
00006 #ifdef __cplusplus
00007 extern "C" {
00008 #endif
00009 #include "f2c.h"
00010
00011 int d132re_(integer *wtleng, doublereal *w, doublereal *g,
00012 doublereal *errcof, doublereal *rulpts)
00013 {
00014
00015
00016 static doublereal dim2g[16] = { .2517129343453109,.7013933644534266,
00017 .9590960631619962,.9956010478552127,.5,.1594544658297559,
00018 .3808991135940188,.6582769255267192,.8761473165029315,
00019 .998243184053198,.9790222658168462,.6492284325645389,
00020 .8727421201131239,.3582614645881228,.5666666666666666,
00021 .2077777777777778 };
00022 static doublereal dim2w[70] = { .0337969236013446,
00023 .09508589607597761,.1176006468056962,.0265777458632695,
00024 .0170144177020064,0.,.0162659309863741,.1344892658526199,
00025 .1328032165460149,.0563747476999187,.0039082790813105,
00026 .0301279877743215,.1030873234689166,.0625,.3213775489050763,
00027 -.1767341636743844,.07347600537466072,-.03638022004364754,
00028 .02125297922098712,.1460984204026913,.01747613286152099,
00029 .1444954045641582,1.307687976001325e-4,5.380992313941161e-4,
00030 1.042259576889814e-4,-.001401152865045733,.008041788181514763,
00031 -.1420416552759383,.3372900883288987,-.1644903060344491,
00032 .07707849911634622,-.0380447835850631,.02223559940380806,
00033 .1480693879765931,4.467143702185814e-6,.150894476707413,
00034 3.647200107516215e-5,5.77719899901388e-4,1.041757313688177e-4,
00035 -.001452822267047819,.008338339968783705,-.147279632923196,
00036 -.8264123822525677,.306583861409436,.002389292538329435,
00037 -.1343024157997222,.088333668405339,0.,9.786283074168292e-4,
00038 -.1319227889147519,.00799001220015063,.003391747079760626,
00039 .002294915718283264,-.01358584986119197,.04025866859057809,
00040 .003760268580063992,.6539094339575232,-.2041614154424632,
00041 -.174698151579499,.03937939671417803,.006974520545933992,0.,
00042 .006667702171778258,.05512960621544304,.05443846381278607,
00043 .02310903863953934,.01506937747477189,-.0605702164890189,
00044 .04225737654686337,.02561989142123099 };
00045
00046 static integer i__, j;
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098 --rulpts;
00099 g -= 3;
00100 w -= 6;
00101 --errcof;
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114 for (i__ = 1; i__ <= 14; ++i__) {
00115 for (j = 1; j <= 5; ++j) {
00116 w[j + i__ * 5] = dim2w[i__ + j * 14 - 15];
00117
00118 }
00119 }
00120
00121
00122
00123 for (i__ = 1; i__ <= 2; ++i__) {
00124 for (j = 1; j <= 14; ++j) {
00125 g[i__ + (j << 1)] = 0.;
00126
00127 }
00128 }
00129 g[5] = dim2g[0];
00130 g[7] = dim2g[1];
00131 g[9] = dim2g[2];
00132 g[11] = dim2g[3];
00133 g[13] = dim2g[4];
00134 g[15] = dim2g[5];
00135 g[16] = g[15];
00136 g[17] = dim2g[6];
00137 g[18] = g[17];
00138 g[19] = dim2g[7];
00139 g[20] = g[19];
00140 g[21] = dim2g[8];
00141 g[22] = g[21];
00142 g[23] = dim2g[9];
00143 g[24] = g[23];
00144 g[25] = dim2g[10];
00145 g[26] = dim2g[11];
00146 g[27] = dim2g[12];
00147 g[28] = dim2g[13];
00148 g[29] = dim2g[14];
00149 g[30] = dim2g[15];
00150
00151
00152
00153 rulpts[1] = 1.;
00154 for (i__ = 2; i__ <= 11; ++i__) {
00155 rulpts[i__] = 4.;
00156
00157 }
00158 rulpts[12] = 8.;
00159 rulpts[13] = 8.;
00160 rulpts[14] = 8.;
00161
00162
00163
00164 errcof[1] = 10.;
00165 errcof[2] = 10.;
00166 errcof[3] = (float)1.;
00167 errcof[4] = (float)5.;
00168 errcof[5] = (float).5;
00169 errcof[6] = (float).25;
00170
00171
00172
00173 return 0;
00174 }
00175
00176 #ifdef __cplusplus
00177 }
00178 #endif