/* TEMPO Version 1.193, 02/03/98, Copyright (c) IRISA/INRIA-Universite de Rennes */


/* 
call signature of entry point romberg:
binding times of formals: r D, a D, b D 
 */
extern float f(float);
extern void romberg_8(float (*)[10], float, float);


extern void romberg_8/*0*/(float (*r)[10], float a, float b)
/* 
binding times of formals: r D, a D, b D 
evaluation time of body: D 
evaluation times of pairings: D, D, D 
 */
  {
    float h;
    float s;

    h = b - a;
    (r[0])[0] = (double)((f(a) + f(b)) * h) / 2.0;
    h = (double)h / 2.0;
    s = 0.0F + f(a + (float)1.000000 * h);
    (r[1])[0] = (double)(r[0])[0] / 2.0 + (double)(h * s);
    (r[1])[1] = (r[1])[0] + (float)0.333333 * ((r[1])[0] - (r[0])[0]);
    h = (double)h / 2.0;
    s = 0.0F + f(a + (float)1.000000 * h) + f(a + (float)3.000000 * h);
    (r[2])[0] = (double)(r[1])[0] / 2.0 + (double)(h * s);
    (r[2])[1] = (r[2])[0] + (float)0.333333 * ((r[2])[0] - (r[1])[0]);
    (r[2])[2] = (r[2])[1] + (float)0.066667 * ((r[2])[1] - (r[1])[1]);
    h = (double)h / 2.0;
    s = 0.0F + f(a + (float)1.000000 * h) + f(a + (float)3.000000 * h) + f(a + (float)5.000000 * h) + f(a + (float)7.000000 * h);
    (r[3])[0] = (double)(r[2])[0] / 2.0 + (double)(h * s);
    (r[3])[1] = (r[3])[0] + (float)0.333333 * ((r[3])[0] - (r[2])[0]);
    (r[3])[2] = (r[3])[1] + (float)0.066667 * ((r[3])[1] - (r[2])[1]);
    (r[3])[3] = (r[3])[2] + (float)0.015873 * ((r[3])[2] - (r[2])[2]);
    h = (double)h / 2.0;
    s = 0.0F + f(a + (float)1.000000 * h) + f(a + (float)3.000000 * h) + f(a + (float)5.000000 * h) + f(a + (float)7.000000 * h) + f(a + (float)9.000000 * h) + f(a + (float)11.000000 * h) + f(a + (float)13.000000 * h) + f(a + (float)15.000000 * h);
    (r[4])[0] = (double)(r[3])[0] / 2.0 + (double)(h * s);
    (r[4])[1] = (r[4])[0] + (float)0.333333 * ((r[4])[0] - (r[3])[0]);
    (r[4])[2] = (r[4])[1] + (float)0.066667 * ((r[4])[1] - (r[3])[1]);
    (r[4])[3] = (r[4])[2] + (float)0.015873 * ((r[4])[2] - (r[3])[2]);
    (r[4])[4] = (r[4])[3] + (float)0.003922 * ((r[4])[3] - (r[3])[3]);
    h = (double)h / 2.0;
    s = 0.0F + f(a + (float)1.000000 * h) + f(a + (float)3.000000 * h) + f(a + (float)5.000000 * h) + f(a + (float)7.000000 * h) + f(a + (float)9.000000 * h) + f(a + (float)11.000000 * h) + f(a + (float)13.000000 * h) + f(a + (float)15.000000 * h) + f(a + (float)17.000000 * h) + f(a + (float)19.000000 * h) + f(a + (float)21.000000 * h) + f(a + (float)23.000000 * h) + f(a + (float)25.000000 * h) + f(a + (float)27.000000 * h) + f(a + (float)29.000000 * h) + f(a + (float)31.000000 * h);
    (r[5])[0] = (double)(r[4])[0] / 2.0 + (double)(h * s);
    (r[5])[1] = (r[5])[0] + (float)0.333333 * ((r[5])[0] - (r[4])[0]);
    (r[5])[2] = (r[5])[1] + (float)0.066667 * ((r[5])[1] - (r[4])[1]);
    (r[5])[3] = (r[5])[2] + (float)0.015873 * ((r[5])[2] - (r[4])[2]);
    (r[5])[4] = (r[5])[3] + (float)0.003922 * ((r[5])[3] - (r[4])[3]);
    (r[5])[5] = (r[5])[4] + (float)0.000978 * ((r[5])[4] - (r[4])[4]);
    h = (double)h / 2.0;
    s = 0.0F + f(a + (float)1.000000 * h) + f(a + (float)3.000000 * h) + f(a + (float)5.000000 * h) + f(a + (float)7.000000 * h) + f(a + (float)9.000000 * h) + f(a + (float)11.000000 * h) + f(a + (float)13.000000 * h) + f(a + (float)15.000000 * h) + f(a + (float)17.000000 * h) + f(a + (float)19.000000 * h) + f(a + (float)21.000000 * h) + f(a + (float)23.000000 * h) + f(a + (float)25.000000 * h) + f(a + (float)27.000000 * h) + f(a + (float)29.000000 * h) + f(a + (float)31.000000 * h) + f(a + (float)33.000000 * h) + f(a + (float)35.000000 * h) + f(a + (float)37.000000 * h) + f(a + (float)39.000000 * h) + f(a + (float)41.000000 * h) + f(a + (float)43.000000 * h) + f(a + (float)45.000000 * h) + f(a + (float)47.000000 * h) + f(a + (float)49.000000 * h) + f(a + (float)51.000000 * h) + f(a + (float)53.000000 * h) + f(a + (float)55.000000 * h) + f(a + (float)57.000000 * h) + f(a + (float)59.000000 * h) + f(a + (float)61.000000 * h) + f(a + (float)63.000000 * h);
    (r[6])[0] = (double)(r[5])[0] / 2.0 + (double)(h * s);
    (r[6])[1] = (r[6])[0] + (float)0.333333 * ((r[6])[0] - (r[5])[0]);
    (r[6])[2] = (r[6])[1] + (float)0.066667 * ((r[6])[1] - (r[5])[1]);
    (r[6])[3] = (r[6])[2] + (float)0.015873 * ((r[6])[2] - (r[5])[2]);
    (r[6])[4] = (r[6])[3] + (float)0.003922 * ((r[6])[3] - (r[5])[3]);
    (r[6])[5] = (r[6])[4] + (float)0.000978 * ((r[6])[4] - (r[5])[4]);
    (r[6])[6] = (r[6])[5] + (float)0.000244 * ((r[6])[5] - (r[5])[5]);
    h = (double)h / 2.0;
    s = 0.0F + f(a + (float)1.000000 * h) + f(a + (float)3.000000 * h) + f(a + (float)5.000000 * h) + f(a + (float)7.000000 * h) + f(a + (float)9.000000 * h) + f(a + (float)11.000000 * h) + f(a + (float)13.000000 * h) + f(a + (float)15.000000 * h) + f(a + (float)17.000000 * h) + f(a + (float)19.000000 * h) + f(a + (float)21.000000 * h) + f(a + (float)23.000000 * h) + f(a + (float)25.000000 * h) + f(a + (float)27.000000 * h) + f(a + (float)29.000000 * h) + f(a + (float)31.000000 * h) + f(a + (float)33.000000 * h) + f(a + (float)35.000000 * h) + f(a + (float)37.000000 * h) + f(a + (float)39.000000 * h) + f(a + (float)41.000000 * h) + f(a + (float)43.000000 * h) + f(a + (float)45.000000 * h) + f(a + (float)47.000000 * h) + f(a + (float)49.000000 * h) + f(a + (float)51.000000 * h) + f(a + (float)53.000000 * h) + f(a + (float)55.000000 * h) + f(a + (float)57.000000 * h) + f(a + (float)59.000000 * h) + f(a + (float)61.000000 * h) + f(a + (float)63.000000 * h) + f(a + (float)65.000000 * h) + f(a + (float)67.000000 * h) + f(a + (float)69.000000 * h) + f(a + (float)71.000000 * h) + f(a + (float)73.000000 * h) + f(a + (float)75.000000 * h) + f(a + (float)77.000000 * h) + f(a + (float)79.000000 * h) + f(a + (float)81.000000 * h) + f(a + (float)83.000000 * h) + f(a + (float)85.000000 * h) + f(a + (float)87.000000 * h) + f(a + (float)89.000000 * h) + f(a + (float)91.000000 * h) + f(a + (float)93.000000 * h) + f(a + (float)95.000000 * h) + f(a + (float)97.000000 * h) + f(a + (float)99.000000 * h) + f(a + (float)101.000000 * h) + f(a + (float)103.000000 * h) + f(a + (float)105.000000 * h) + f(a + (float)107.000000 * h) + f(a + (float)109.000000 * h) + f(a + (float)111.000000 * h) + f(a + (float)113.000000 * h) + f(a + (float)115.000000 * h) + f(a + (float)117.000000 * h) + f(a + (float)119.000000 * h) + f(a + (float)121.000000 * h) + f(a + (float)123.000000 * h) + f(a + (float)125.000000 * h) + f(a + (float)127.000000 * h);
    (r[7])[0] = (double)(r[6])[0] / 2.0 + (double)(h * s);
    (r[7])[1] = (r[7])[0] + (float)0.333333 * ((r[7])[0] - (r[6])[0]);
    (r[7])[2] = (r[7])[1] + (float)0.066667 * ((r[7])[1] - (r[6])[1]);
    (r[7])[3] = (r[7])[2] + (float)0.015873 * ((r[7])[2] - (r[6])[2]);
    (r[7])[4] = (r[7])[3] + (float)0.003922 * ((r[7])[3] - (r[6])[3]);
    (r[7])[5] = (r[7])[4] + (float)0.000978 * ((r[7])[4] - (r[6])[4]);
    (r[7])[6] = (r[7])[5] + (float)0.000244 * ((r[7])[5] - (r[6])[5]);
    (r[7])[7] = (r[7])[6] + (float)0.000061 * ((r[7])[6] - (r[6])[6]);
    h = (double)h / 2.0;
    s = 0.0F + f(a + (float)1.000000 * h) + f(a + (float)3.000000 * h) + f(a + (float)5.000000 * h) + f(a + (float)7.000000 * h) + f(a + (float)9.000000 * h) + f(a + (float)11.000000 * h) + f(a + (float)13.000000 * h) + f(a + (float)15.000000 * h) + f(a + (float)17.000000 * h) + f(a + (float)19.000000 * h) + f(a + (float)21.000000 * h) + f(a + (float)23.000000 * h) + f(a + (float)25.000000 * h) + f(a + (float)27.000000 * h) + f(a + (float)29.000000 * h) + f(a + (float)31.000000 * h) + f(a + (float)33.000000 * h) + f(a + (float)35.000000 * h) + f(a + (float)37.000000 * h) + f(a + (float)39.000000 * h) + f(a + (float)41.000000 * h) + f(a + (float)43.000000 * h) + f(a + (float)45.000000 * h) + f(a + (float)47.000000 * h) + f(a + (float)49.000000 * h) + f(a + (float)51.000000 * h) + f(a + (float)53.000000 * h) + f(a + (float)55.000000 * h) + f(a + (float)57.000000 * h) + f(a + (float)59.000000 * h) + f(a + (float)61.000000 * h) + f(a + (float)63.000000 * h) + f(a + (float)65.000000 * h) + f(a + (float)67.000000 * h) + f(a + (float)69.000000 * h) + f(a + (float)71.000000 * h) + f(a + (float)73.000000 * h) + f(a + (float)75.000000 * h) + f(a + (float)77.000000 * h) + f(a + (float)79.000000 * h) + f(a + (float)81.000000 * h) + f(a + (float)83.000000 * h) + f(a + (float)85.000000 * h) + f(a + (float)87.000000 * h) + f(a + (float)89.000000 * h) + f(a + (float)91.000000 * h) + f(a + (float)93.000000 * h) + f(a + (float)95.000000 * h) + f(a + (float)97.000000 * h) + f(a + (float)99.000000 * h) + f(a + (float)101.000000 * h) + f(a + (float)103.000000 * h) + f(a + (float)105.000000 * h) + f(a + (float)107.000000 * h) + f(a + (float)109.000000 * h) + f(a + (float)111.000000 * h) + f(a + (float)113.000000 * h) + f(a + (float)115.000000 * h) + f(a + (float)117.000000 * h) + f(a + (float)119.000000 * h) + f(a + (float)121.000000 * h) + f(a + (float)123.000000 * h) + f(a + (float)125.000000 * h) + f(a + (float)127.000000 * h) + f(a + (float)129.000000 * h) + f(a + (float)131.000000 * h) + f(a + (float)133.000000 * h) + f(a + (float)135.000000 * h) + f(a + (float)137.000000 * h) + f(a + (float)139.000000 * h) + f(a + (float)141.000000 * h) + f(a + (float)143.000000 * h) + f(a + (float)145.000000 * h) + f(a + (float)147.000000 * h) + f(a + (float)149.000000 * h) + f(a + (float)151.000000 * h) + f(a + (float)153.000000 * h) + f(a + (float)155.000000 * h) + f(a + (float)157.000000 * h) + f(a + (float)159.000000 * h) + f(a + (float)161.000000 * h) + f(a + (float)163.000000 * h) + f(a + (float)165.000000 * h) + f(a + (float)167.000000 * h) + f(a + (float)169.000000 * h) + f(a + (float)171.000000 * h) + f(a + (float)173.000000 * h) + f(a + (float)175.000000 * h) + f(a + (float)177.000000 * h) + f(a + (float)179.000000 * h) + f(a + (float)181.000000 * h) + f(a + (float)183.000000 * h) + f(a + (float)185.000000 * h) + f(a + (float)187.000000 * h) + f(a + (float)189.000000 * h) + f(a + (float)191.000000 * h) + f(a + (float)193.000000 * h) + f(a + (float)195.000000 * h) + f(a + (float)197.000000 * h) + f(a + (float)199.000000 * h) + f(a + (float)201.000000 * h) + f(a + (float)203.000000 * h) + f(a + (float)205.000000 * h) + f(a + (float)207.000000 * h) + f(a + (float)209.000000 * h) + f(a + (float)211.000000 * h) + f(a + (float)213.000000 * h) + f(a + (float)215.000000 * h) + f(a + (float)217.000000 * h) + f(a + (float)219.000000 * h) + f(a + (float)221.000000 * h) + f(a + (float)223.000000 * h) + f(a + (float)225.000000 * h) + f(a + (float)227.000000 * h) + f(a + (float)229.000000 * h) + f(a + (float)231.000000 * h) + f(a + (float)233.000000 * h) + f(a + (float)235.000000 * h) + f(a + (float)237.000000 * h) + f(a + (float)239.000000 * h) + f(a + (float)241.000000 * h) + f(a + (float)243.000000 * h) + f(a + (float)245.000000 * h) + f(a + (float)247.000000 * h) + f(a + (float)249.000000 * h) + f(a + (float)251.000000 * h) + f(a + (float)253.000000 * h) + f(a + (float)255.000000 * h);
    (r[8])[0] = (double)(r[7])[0] / 2.0 + (double)(h * s);
    (r[8])[1] = (r[8])[0] + (float)0.333333 * ((r[8])[0] - (r[7])[0]);
    (r[8])[2] = (r[8])[1] + (float)0.066667 * ((r[8])[1] - (r[7])[1]);
    (r[8])[3] = (r[8])[2] + (float)0.015873 * ((r[8])[2] - (r[7])[2]);
    (r[8])[4] = (r[8])[3] + (float)0.003922 * ((r[8])[3] - (r[7])[3]);
    (r[8])[5] = (r[8])[4] + (float)0.000978 * ((r[8])[4] - (r[7])[4]);
    (r[8])[6] = (r[8])[5] + (float)0.000244 * ((r[8])[5] - (r[7])[5]);
    (r[8])[7] = (r[8])[6] + (float)0.000061 * ((r[8])[6] - (r[7])[6]);
    (r[8])[8] = (r[8])[7] + (float)0.000015 * ((r[8])[7] - (r[7])[7]);
    return;
  }