00001
00002
00003
00004
00005
00006 #ifdef __cplusplus
00007 extern "C" {
00008 #endif
00009 #include "f2c.h"
00010
00011 int d113re_(integer *wtleng, doublereal *w, doublereal *g,
00012 doublereal *errcof, doublereal *rulpts)
00013 {
00014
00015
00016 static doublereal dim3g[14] = { .19,.5,.75,.8,.9949999999999999,
00017 .99873449983514,.7793703685672423,.9999698993088767,
00018 .7902637224771788,.4403396687650737,.4378478459006862,
00019 .9549373822794593,.9661093133630748,.4577105877763134 };
00020 static doublereal dim3w[65] = { .007923078151105734,
00021 .0679717739278808,.001086986538805825,.1838633662212829,
00022 .03362119777829031,.01013751123334062,.001687648683985235,
00023 .1346468564512807,.001750145884600386,.07752336383837454,
00024 .2461864902770251,.06797944868483039,.01419962823300713,
00025 1.715006248224684,-.3755893815889209,.1488632145140549,
00026 -.2497046640620823,.1792501419135204,.00344612675897389,
00027 -.005140483185555825,.006536017839876425,-6.5134549392297e-4,
00028 -.006304672433547204,.01266959399788263,-.005454241018647931,
00029 .004826995274768427,1.936014978949526,-.3673449403754268,
00030 .02929778657898176,-.1151883520260315,.05086658220872218,
00031 .04453911087786469,-.022878282571259,.02908926216345833,
00032 -.002898884350669207,-.02805963413307495,.05638741361145884,
00033 -.02427469611942451,.02148307034182882,.517082819560576,
00034 .01445269144914044,-.3601489663995932,.3628307003418485,
00035 .007148802650872729,-.09222852896022966,.01719339732471725,
00036 -.102141653746035,-.007504397861080493,.01648362537726711,
00037 .05234610158469334,.01445432331613066,.003019236275367777,
00038 2.05440450381852,.0137775998849012,-.576806291790441,
00039 .03726835047700328,.006814878939777219,.05723169733851849,
00040 -.04493018743811285,.02729236573866348,3.54747395055699e-4,
00041 .01571366799739551,.04990099219278567,.0137791555266677,
00042 .002878206423099872 };
00043
00044 static integer i__, j;
00045
00046
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
00099
00100
00101 --rulpts;
00102 g -= 4;
00103 w -= 6;
00104 --errcof;
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114
00115
00116
00117 for (i__ = 1; i__ <= 13; ++i__) {
00118 for (j = 1; j <= 5; ++j) {
00119 w[j + i__ * 5] = dim3w[i__ + j * 13 - 14];
00120
00121 }
00122 }
00123
00124
00125
00126 for (i__ = 1; i__ <= 3; ++i__) {
00127 for (j = 1; j <= 13; ++j) {
00128 g[i__ + j * 3] = 0.;
00129
00130 }
00131 }
00132 g[7] = dim3g[0];
00133 g[10] = dim3g[1];
00134 g[13] = dim3g[2];
00135 g[16] = dim3g[3];
00136 g[19] = dim3g[4];
00137 g[22] = dim3g[5];
00138 g[23] = g[22];
00139 g[25] = dim3g[6];
00140 g[26] = g[25];
00141 g[28] = dim3g[7];
00142 g[29] = g[28];
00143 g[30] = g[28];
00144 g[31] = dim3g[8];
00145 g[32] = g[31];
00146 g[33] = g[31];
00147 g[34] = dim3g[9];
00148 g[35] = g[34];
00149 g[36] = g[34];
00150 g[37] = dim3g[11];
00151 g[38] = dim3g[10];
00152 g[39] = g[38];
00153 g[40] = dim3g[12];
00154 g[41] = g[40];
00155 g[42] = dim3g[13];
00156
00157
00158
00159 rulpts[1] = 1.;
00160 rulpts[2] = 6.;
00161 rulpts[3] = 6.;
00162 rulpts[4] = 6.;
00163 rulpts[5] = 6.;
00164 rulpts[6] = 6.;
00165 rulpts[7] = 12.;
00166 rulpts[8] = 12.;
00167 rulpts[9] = 8.;
00168 rulpts[10] = 8.;
00169 rulpts[11] = 8.;
00170 rulpts[12] = 24.;
00171 rulpts[13] = 24.;
00172
00173
00174
00175 errcof[1] = 4.;
00176 errcof[2] = (float)4.;
00177 errcof[3] = (float).5;
00178 errcof[4] = (float)3.;
00179 errcof[5] = (float).5;
00180 errcof[6] = (float).25;
00181
00182
00183
00184 return 0;
00185 }
00186
00187 #ifdef __cplusplus
00188 }
00189 #endif