<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#ffffff" text="#000000">
    Hi Séverine,<br>
    <br>
    The string you can pass to scilab has a limited size.<br>
    <br>
    A workaround:<br>
    i) put your big line in a file: expression.txt<br>
    ii) get your line with : l=mgetl('expression.txt');<br>
    iii) J(i,k)=evstr(l);<br>
    <br>
    Best regards,<br>
    <br>
    Calixte<br>
    <br>
    <br>
    On 03/10/2011 15:53, Séverine Paul wrote:
    <blockquote
cite="mid:CAHvBYcN+sNjsOFOzQBLbfQ4z4WXx+AvZ0VV+gYb=e4exxjyPAw@mail.gmail.com"
      type="cite">Hi again,<br>
      <br>
      this is the answer of scilab:<br>
      <br>
      Command is too long (more than 512 characters long): could not
      send it to Scilab<br>
      <br>
      and below is the command. It is actually quite long, but I didn't
      know that this could be a problem. In Maple (from which this
      comes), there is no problem<br>
      <br>
      Thanks for your help.<br>
      Séverine.<br>
      <br>
      J(i,k) =<br>
      -exp(-b*(x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      (-0.144e3 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i)
      ^ 2)) * sqrt(b) + 0.144e3 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k
      + 1) + x(i + 1) ^ 2)) * sqrt(b) - 0.56e2 * exp(b * (x(k + 1) ^ 2 +
      0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(k) ^ 2 * x(i) * x(i + 1)
      * b ^ (0.5e1 / 0.2e1) - 0.66e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i)
      * x(k + 1) + x(i + 1) ^ 2)) * x(i) ^ 2 * b ^ (0.3e1 / 0.2e1) -
      0.420e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i) ^ 4 * x(i +
      1) * sqrt(%pi) * b ^ 3 - 0.168e3 * b ^ (0.7e1 / 0.2e1) * exp(b *
      (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i) * x(k
      + 1) ^ 4 * x(k) + 0.32e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) *
      x(k) + x(i + 1) ^ 2)) * x(i) ^ 5 * x(k) * b ^ (0.7e1 / 0.2e1) +
      0.32e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^
      2)) * x(i) ^ 4 * x(k) ^ 2 * b ^ (0.7e1 / 0.2e1) + 0.112e3 * exp(b
      * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i) ^ 5
      * x(i + 1) * b ^ (0.7e1 / 0.2e1) + 0.112e3 * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i) * x(i + 1) ^ 5 * b
      ^ (0.7e1 / 0.2e1) + 0.32e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1)
      * x(k + 1) + x(i) ^ 2)) * x(i + 1) ^ 5 * x(k + 1) * b ^ (0.7e1 /
      0.2e1) - 0.126e3 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1)
      + x(i) ^ 2)) * x(i + 1) * x(i) * b ^ (0.3e1 / 0.2e1) + 0.48e2 *
      exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(k
      + 1) * x(i) ^ 3 * b ^ (0.5e1 / 0.2e1) + 0.168e3 * b ^ 4 * exp(b *
      (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i +
      1) * sqrt(b) - x(k) * sqrt(b)) * x(i + 1) ^ 5 * x(i) * sqrt(%pi) *
      x(k + 1) + 0.196e3 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(i) ^ 2 *
      x(k) - 0.56e2 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 +
      0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) * x(k) ^ 2 * x(k +
      1) + 0.56e2 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) * x(k + 1) ^ 2 * x(k)
      + 0.280e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i +
      1) * sqrt(%pi) * x(i) * x(k + 1) ^ 4 * x(k) - 0.112e3 * b ^ 4 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i) ^ 6 * sqrt(%pi) *
      x(k + 1) + 0.40e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) +
      x(i + 1) ^ 2)) * x(k + 1) ^ 2 * x(i) ^ 2 * b ^ (0.5e1 / 0.2e1) -
      0.14e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^
      2)) * x(k + 1) * x(i + 1) * b ^ (0.3e1 / 0.2e1) - 0.105e3 * exp(b
      * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i) *
      sqrt(%pi) * erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b)) * b +
      0.212e3 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^
      2)) * x(k) ^ 4 * b ^ (0.5e1 / 0.2e1) - 0.212e3 * exp(b * (x(k + 1)
      ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(k) ^ 4 * b ^ (0.5e1
      / 0.2e1) - 0.48e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) +
      x(i + 1) ^ 2)) * x(k) ^ 3 * x(i) * b ^ (0.5e1 / 0.2e1) + 0.364e3 *
      exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) *
      x(i) * x(i + 1) ^ 3 * b ^ (0.5e1 / 0.2e1) - 0.140e3 * b ^ 3 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * sqrt(%pi) * x(i) ^ 4 *
      x(k + 1) - 0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k +
      1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b))
      * x(i + 1) * sqrt(%pi) * x(k + 1) ^ 5 * x(k) + 0.140e3 * exp(b *
      (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(k) ^ 4 *
      erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i) * sqrt(%pi) * b ^
      3 + 0.56e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i +
      1) ^ 2)) * x(k) * x(i + 1) * x(i) ^ 2 * b ^ (0.5e1 / 0.2e1) -
      0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) *
      x(k + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(i) * x(k + 1) ^ 3 * x(k)
      - 0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(i) ^ 5 *
      x(i + 1) * sqrt(%pi) * x(k + 1) - 0.168e3 * b ^ (0.5e1 / 0.2e1) *
      exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i
      + 1) * x(i) * x(k) * x(k + 1) - 0.56e2 * b ^ (0.7e1 / 0.2e1) *
      exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) *
      x(i) ^ 3 * x(k) ^ 2 * x(k + 1) - 0.196e3 * b ^ (0.5e1 / 0.2e1) *
      exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) *
      x(i) * x(k + 1) ^ 2 * x(k) - 0.210e3 * b ^ 2 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i + 1) * x(i) *
      sqrt(%pi) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * x(k) +
      0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i +
      1) * x(k + 1) + x(i) ^ 2)) * x(i) * x(i + 1) ^ 3 * x(k + 1) * x(k)
      + 0.144e3 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i)
      ^ 2)) * sqrt(b) - 0.252e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k +
      1) ^ 2 + x(i + 1) ^ 2)) * x(k) ^ 5 * erf(x(i + 1) * sqrt(b) - x(k)
      * sqrt(b)) * sqrt(%pi) * b ^ 3 + 0.280e3 * b ^ 3 * exp(b * (x(k) ^
      2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) *
      sqrt(b) - x(k) * sqrt(b)) * x(i + 1) ^ 3 * x(i) * sqrt(%pi) * x(k)
      + 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i
      + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) * x(i) * x(k + 1) ^ 3 *
      x(k) + 0.56e2 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 +
      0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i) ^ 2 * x(k) * x(k + 1)
      - 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i)
      * x(k + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(i) ^ 2 * x(k + 1) ^ 2 *
      x(k) + 0.210e3 * b ^ 2 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^
      2 + x(i + 1) ^ 2)) * x(i + 1) * x(i) * sqrt(%pi) * erf(x(i + 1) *
      sqrt(b) - x(k + 1) * sqrt(b)) * x(k) - 0.112e3 * b ^ (0.7e1 /
      0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^
      2)) * x(i) * x(i + 1) ^ 2 * x(k) ^ 2 * x(k + 1) - 0.40e2 * exp(b *
      (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(k + 1) ^
      2 * x(i + 1) ^ 2 * b ^ (0.5e1 / 0.2e1) - 0.364e3 * exp(b * (x(k) ^
      2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(i) ^ 3
      * b ^ (0.5e1 / 0.2e1) - 0.140e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 +
      x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(k) ^ 4 * erf(x(i) * sqrt(b) -
      x(k) * sqrt(b)) * x(i + 1) * sqrt(%pi) * b ^ 3 + 0.140e3 * exp(b *
      (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(k) ^ 4 *
      erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i + 1) * sqrt(%pi) *
      b ^ 3 + 0.56e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) +
      x(i) ^ 2)) * x(k + 1) * x(i) * x(i + 1) ^ 2 * b ^ (0.5e1 / 0.2e1)
      - 0.56e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1)
      ^ 2)) * x(i) ^ 4 * x(i + 1) * x(k) * b ^ (0.7e1 / 0.2e1) + 0.112e3
      * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i + 1) ^ 6 * x(i)
      * sqrt(%pi) * b ^ 4 - 0.112e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^
      2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k) *
      sqrt(b)) * sqrt(%pi) * x(k) ^ 6 * x(k + 1) + 0.224e3 * b ^ (0.3e1
      / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1)
      ^ 2)) * x(i + 1) * x(k + 1) + 0.14e2 * b ^ (0.3e1 / 0.2e1) * exp(b
      * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i) *
      x(k + 1) + 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 +
      0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(k) ^ 5 * x(k + 1) +
      0.56e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^
      2)) * x(i + 1) ^ 3 * x(i) * x(k) ^ 2 * b ^ (0.7e1 / 0.2e1) +
      0.66e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^
      2)) * x(i) ^ 2 * b ^ (0.3e1 / 0.2e1) - 0.56e2 * exp(b * (x(k + 1)
      ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i + 1) * x(i) *
      x(k) ^ 4 * b ^ (0.7e1 / 0.2e1) + 0.364e3 * exp(b * (x(k + 1) ^ 2 +
      0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i + 1) * x(i) ^ 3 * b ^
      (0.5e1 / 0.2e1) - 0.92e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i +
      1) * x(k) + x(i) ^ 2)) * x(k) * x(i + 1) * b ^ (0.3e1 / 0.2e1) +
      0.126e3 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^
      2)) * x(i + 1) * x(i) * b ^ (0.3e1 / 0.2e1) - 0.56e2 * exp(b *
      (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i) * x(i
      + 1) ^ 3 * x(k + 1) ^ 2 * b ^ (0.7e1 / 0.2e1) + 0.14e2 * exp(b *
      (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(k + 1) *
      x(i) * b ^ (0.3e1 / 0.2e1) - 0.80e2 * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) ^ 6 * b ^ (0.7e1 /
      0.2e1) + 0.80e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) +
      x(i + 1) ^ 2)) * x(i) ^ 6 * b ^ (0.7e1 / 0.2e1) + 0.105e3 * b *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      sqrt(%pi) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(k + 1) +
      0.126e3 * b ^ (0.3e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) *
      x(k + 1) + x(i + 1) ^ 2)) * x(k + 1) * x(k)+0.84e2 * b ^ (0.5e1 /
      0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^
      2)) * x(i + 1) ^ 3 * x(k) + 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b
      * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i) ^ 5
      * x(k + 1) - 0.126e3 * b ^ (0.3e1 / 0.2e1) * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(k + 1) * x(k) -
      0.364e3 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 *
      x(i + 1) * x(k) + x(i) ^ 2)) * x(k) ^ 3 * x(k + 1) + 0.364e3 * b ^
      (0.5e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) +
      x(i + 1) ^ 2)) * x(k) ^ 3 * x(k + 1) + 0.84e2 * b ^ (0.5e1 /
      0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^
      2)) * x(i) ^ 3 * x(k + 1) - 0.364e3 * b ^ (0.5e1 / 0.2e1) * exp(b
      * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(k + 1)
      ^ 3 * x(k) + 0.224e3 * b ^ (0.3e1 / 0.2e1) * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i) * x(k) - 0.32e2 *
      exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) *
      x(i) ^ 5 * x(k + 1) * b ^ (0.7e1 / 0.2e1) + 0.56e2 * exp(b * (x(k
      + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) * x(i)
      * x(k) ^ 4 * b ^ (0.7e1 / 0.2e1) - 0.420e3 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b)
      - x(k) * sqrt(b)) * x(i + 1) ^ 4 * x(i) * sqrt(%pi) * b ^ 3 -
      0.80e2 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(i) ^ 7 * sqrt(%pi)
      * b ^ 4 + 0.280e3 * x(k + 1) * x(k) * b ^ 4 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b)
      - x(k + 1) * sqrt(b)) * x(i) * sqrt(%pi) * x(i + 1) ^ 4 - 0.168e3
      * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(i) ^ 5 * x(i + 1) *
      sqrt(%pi) * x(k) - 0.168e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k +
      1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i + 1) * x(k) ^
      4 * x(k + 1) + 0.56e2 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i) ^ 3 * x(k + 1) ^ 2
      * x(k) + 0.105e3 * b * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2
      + x(i + 1) ^ 2)) * sqrt(%pi) * erf(x(i + 1) * sqrt(b) - x(k + 1) *
      sqrt(b)) * x(k) + 0.56e2 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k + 1)
      ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 3 * x(k) ^
      2 * x(k + 1) + 0.280e3 * b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 +
      x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) *
      sqrt(b)) * x(i + 1) * sqrt(%pi) * x(i) ^ 3 * x(k + 1) + 0.140e3 *
      b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(i) ^ 4 * sqrt(%pi)
      * x(k) + 0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1)
      ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k + 1) *
      sqrt(b)) * x(i) * sqrt(%pi) * x(k + 1) ^ 5 * x(k) + 0.196e3 * b ^
      (0.5e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) +
      x(i + 1) ^ 2)) * x(i + 1) * x(k + 1) ^ 2 * x(k) - 0.56e2 * b ^
      (0.7e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) +
      x(i + 1) ^ 2)) * x(i) ^ 4 * x(k) * x(k + 1) + 0.196e3 * b ^ (0.5e1
      / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i)
      ^ 2)) * x(i) * x(k) ^ 2 * x(k + 1) + 0.32e2 * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) ^ 2 * x(k + 1)
      ^ 4 * b ^ (0.7e1 / 0.2e1) + 0.280e3 * b ^ 3 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) -
      x(k) * sqrt(b)) * x(i) * sqrt(%pi) * x(k) ^ 3 * x(k + 1) - 0.280e3
      * b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * x(i + 1) ^ 3 * x(i) * sqrt(%pi) * erf(x(i + 1) * sqrt(b) -
      x(k + 1) * sqrt(b)) * x(k + 1) - 0.168e3 * b ^ 4 * exp(b * (x(k) ^
      2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) *
      sqrt(b) - x(k + 1) * sqrt(b)) * x(i + 1) ^ 5 * x(i) * sqrt(%pi) *
      x(k) + 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) ^ 2 * x(i) * x(k + 1)
      ^ 2 * x(k) - 0.280e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k
      + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k + 1) *
      sqrt(b)) * x(i + 1) * x(i) * sqrt(%pi) * x(k + 1) ^ 4 * x(k) -
      0.112e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i +
      1) ^ 6 * sqrt(%pi) * x(k + 1) + 0.105e3 * exp(b * (x(k) ^ 2 + x(i)
      ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i + 1) * sqrt(%pi) *
      erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * b - 0.280e3 * b ^ 3 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      x(i) ^ 3 * x(i + 1) * sqrt(%pi) * erf(x(i) * sqrt(b) - x(k) *
      sqrt(b)) * x(k) + 0.56e2 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k + 1)
      ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 4 * x(k) *
      x(k + 1) + 0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k +
      1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b))
      * x(i) ^ 5 * x(i + 1) * sqrt(%pi) * x(k) + 0.210e3 * b ^ 2 * exp(b
      * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i + 1)
      * x(i) * sqrt(%pi) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(k +
      1) + 0.168e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(k + 1) ^ 4 * x(k)
      - 0.84e2 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 *
      x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 3 * x(k + 1) - 0.280e3 *
      b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i + 1) *
      sqrt(%pi) * x(k) ^ 3 * x(k + 1) + 0.420e3 * b ^ 3 * exp(b * (x(k)
      ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) *
      sqrt(b) - x(k) * sqrt(b)) * sqrt(%pi) * x(k) ^ 4 * x(k + 1) -
      0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i) *
      sqrt(%pi) * x(k) ^ 5 * x(k + 1) - 0.280e3 * b ^ 3 * exp(b * (x(k)
      ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) *
      sqrt(b) - x(k + 1) * sqrt(b)) * x(i + 1) * sqrt(%pi) * x(k + 1) ^
      3 * x(k) - 0.210e3 * b ^ 2 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k +
      1) ^ 2 + x(i + 1) ^ 2)) * x(i + 1) * x(i) * sqrt(%pi) * erf(x(i) *
      sqrt(b) - x(k + 1) * sqrt(b)) * x(k + 1) + 0.168e3 * exp(b * (x(k)
      ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) *
      sqrt(b) - x(k) * sqrt(b)) * x(i + 1) * x(i) * sqrt(%pi) * x(k) ^ 5
      * b ^ 4 - 0.252e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(i) ^ 5 *
      sqrt(%pi) * b ^ 3 - 0.140e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k +
      1) ^ 2 + x(i + 1) ^ 2)) * x(k + 1) ^ 4 * erf(x(i + 1) * sqrt(b) -
      x(k + 1) * sqrt(b)) * x(i + 1) * sqrt(%pi) * b ^ 3 - 0.80e2 *
      exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(k
      + 1) ^ 6 * b ^ (0.7e1 / 0.2e1) + 0.56e2 * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(k + 1) ^ 2 * x(i) *
      x(i + 1) * b ^ (0.5e1 / 0.2e1) + 0.56e2 * b ^ (0.5e1 / 0.2e1) *
      exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i
      + 1) ^ 2 * x(k + 1) * x(k) - 0.210e3 * b ^ 2 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i + 1) * x(i) *
      sqrt(%pi) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(k + 1) +
      0.280e3 * b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i +
      1) * sqrt(%pi) * x(i) ^ 3 * x(k) - 0.420e3 * b ^ 3 * exp(b * (x(k)
      ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) *
      sqrt(b) - x(k + 1) * sqrt(b)) * sqrt(%pi) * x(k + 1) ^ 4 * x(k) -
      0.80e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^
      2)) * x(i) ^ 6 * b ^ (0.7e1 / 0.2e1) + 0.56e2 * b ^ (0.7e1 /
      0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^
      2)) * x(i) * x(k + 1) ^ 4 * x(k) + 0.140e3 * b ^ 3 * exp(b * (x(k)
      ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) *
      sqrt(b) - x(k + 1) * sqrt(b)) * x(i + 1) ^ 4 * sqrt(%pi) * x(k +
      1) + 0.168e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1
      * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 4 * x(i) * x(k + 1) +
      0.140e3 * b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(i) ^ 4 *
      sqrt(%pi) * x(k + 1) - 0.168e3 * b ^ (0.7e1 / 0.2e1) * exp(b *
      (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i) * x(i
      + 1) ^ 4 * x(k) + 0.112e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 +
      x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k + 1)
      * sqrt(b)) * x(i + 1) ^ 6 * sqrt(%pi) * x(k) - 0.196e3 * b ^
      (0.5e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) +
      x(i + 1) ^ 2)) * x(i + 1) * x(k) ^ 2 * x(k + 1) + 0.212e3 * exp(b
      * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(k + 1)
      ^ 4 * b ^ (0.5e1 / 0.2e1) - 0.66e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1
      * x(i + 1) * x(k) + x(i) ^ 2)) * x(k) ^ 2 * b ^ (0.3e1 / 0.2e1) -
      0.105e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * x(i + 1) * sqrt(%pi) * erf(x(i + 1) * sqrt(b) - x(k + 1) *
      sqrt(b)) * b - 0.14e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) *
      x(k) + x(i) ^ 2)) * x(k) * x(i) * b ^ (0.3e1 / 0.2e1) - 0.40e2 *
      exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) *
      x(k) ^ 2 * x(i) ^ 2 * b ^ (0.5e1 / 0.2e1) + 0.56e2 * exp(b * (x(k
      + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 4 *
      x(i) * x(k) * b ^ (0.7e1 / 0.2e1) - 0.56e2 * b ^ (0.7e1 / 0.2e1) *
      exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i
      + 1) ^ 4 * x(k + 1) * x(k) - 0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b)
      - x(k + 1) * sqrt(b)) * x(i + 1) ^ 5 * x(i) * sqrt(%pi) * x(k + 1)
      - 0.14e2 * b ^ (0.3e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i)
      * x(k + 1) + x(i + 1) ^ 2)) * x(i) * x(k) - 0.14e2 * b ^ (0.3e1 /
      0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^
      2)) * x(i + 1) * x(k + 1) + 0.126e3 * b ^ (0.3e1 / 0.2e1) * exp(b
      * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(k) *
      x(k + 1) + 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) ^ 5 * x(k) +
      0.364e3 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i +
      1) * x(k + 1) + x(i) ^ 2)) * x(k + 1) ^ 3 * x(k) + 0.140e3 * exp(b
      * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(k + 1)
      ^ 4 * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i + 1) *
      sqrt(%pi) * b ^ 3 + 0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2
      + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k +
      1) * sqrt(b)) * sqrt(%pi) * x(i + 1) * x(k + 1) ^ 5 * x(k) -
      0.105e3 * b * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i +
      1) ^ 2)) * sqrt(%pi) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) *
      x(k) - 0.105e3 * b * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * sqrt(%pi) * erf(x(i) * sqrt(b) - x(k + 1) *
      sqrt(b)) * x(k + 1) - 0.210e3 * x(k + 1) * x(k) * b ^ 2 * exp(b *
      (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i) *
      sqrt(%pi) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) + 0.105e3 * exp(b
      * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i) *
      sqrt(%pi) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * b + 0.92e2
      * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) *
      x(k + 1) * x(i + 1) * b ^ (0.3e1 / 0.2e1) + 0.84e2 * exp(b * (x(k
      + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(k) ^ 3 * x(i +
      1) * b ^ (0.5e1 / 0.2e1) - 0.92e2 * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(k + 1) * x(i) * b ^ (0.3e1 /
      0.2e1) + 0.420e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b)) *
      x(i) * sqrt(%pi) * x(i + 1) ^ 4 * b ^ 3 + 0.66e2 * exp(b * (x(k) ^
      2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(k + 1) ^ 2 * b ^
      (0.3e1 / 0.2e1) + 0.66e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) *
      x(k) + x(i + 1) ^ 2)) * x(k) ^ 2 * b ^ (0.3e1 / 0.2e1) + 0.252e3 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(k
      + 1) ^ 5 * erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b)) *
      sqrt(%pi) * b ^ 3 + 0.210e3 * b ^ 2 * exp(b * (x(k) ^ 2 + x(i) ^ 2
      + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i + 1) * x(i) * sqrt(%pi) *
      erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b)) * x(k + 1) + 0.112e3
      * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) *
      x(k) + x(i + 1) ^ 2)) * x(i + 1) * x(i) ^ 3 * x(k) * x(k + 1) -
      0.56e2 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 *
      x(i) * x(k) + x(i + 1) ^ 2)) * x(i) * x(k) ^ 4 * x(k + 1) +
      0.168e3 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) *
      x(k + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(i) * x(k + 1) * x(k) +
      0.112e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) *
      sqrt(%pi) * x(k) ^ 6 * x(k + 1) - 0.420e3 * b ^ 3 * exp(b * (x(k)
      ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) *
      sqrt(b) - x(k) * sqrt(b)) * sqrt(%pi) * x(k) ^ 4 * x(k + 1) -
      0.66e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^
      2)) * x(k + 1) ^ 2 * b ^ (0.3e1 / 0.2e1) + 0.66e2 * exp(b * (x(k)
      ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) ^ 2 * b
      ^ (0.3e1 / 0.2e1) + 0.212e3 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) *
      x(k + 1) + x(i + 1) ^ 2)) * x(i) ^ 4 * b ^ (0.5e1 / 0.2e1) -
      0.252e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * x(k + 1) ^ 5 * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) *
      sqrt(%pi) * b ^ 3 - 0.105e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k +
      1) ^ 2 + x(i + 1) ^ 2)) * x(i) * sqrt(%pi) * erf(x(i) * sqrt(b) -
      x(k) * sqrt(b)) * b + 0.32e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i +
      1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) * x(k + 1) ^ 5 * b ^ (0.7e1
      / 0.2e1) - 0.112e3 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) +
      x(i + 1) ^ 2)) * x(i + 1) * x(k + 1) ^ 5 * b ^ (0.7e1 / 0.2e1) -
      0.32e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^
      2)) * x(i + 1) ^ 5 * x(k) * b ^ (0.7e1 / 0.2e1) + 0.32e2 * exp(b *
      (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i) * x(k)
      ^ 5 * b ^ (0.7e1 / 0.2e1) - 0.32e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1
      * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) * x(k) ^ 5 * b ^ (0.7e1
      / 0.2e1) - 0.48e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k +
      1) + x(i) ^ 2)) * x(k + 1) ^ 3 * x(i + 1) * b ^ (0.5e1 / 0.2e1) +
      0.92e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^
      2)) * x(k) * x(i) * b ^ (0.3e1 / 0.2e1) + 0.126e3 * exp(b * (x(k +
      1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) * x(i) *
      b ^ (0.3e1 / 0.2e1) - 0.112e3 * exp(b * (x(k + 1) ^ 2 + 0.2e1 *
      x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 5 * x(i) * b ^ (0.7e1 /
      0.2e1) - 0.32e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k)
      + x(i) ^ 2)) * x(i + 1) ^ 3 * x(k) ^ 3 * b ^ (0.7e1 / 0.2e1) +
      0.112e3 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^
      2)) * x(i + 1) * x(k) ^ 5 * b ^ (0.7e1 / 0.2e1) + 0.80e2 * exp(b *
      (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(k + 1) ^
      6 * b ^ (0.7e1 / 0.2e1) + 0.56e2 * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(i) ^ 3 * x(k + 1)
      ^ 2 * b ^ (0.7e1 / 0.2e1) - 0.66e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1
      * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 2 * b ^ (0.3e1 /
      0.2e1) + 0.252e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i +
      1) ^ 5 * sqrt(%pi) * b ^ 3 - 0.80e2 * exp(b * (x(k + 1) ^ 2 +
      0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(k) ^ 6 * b ^ (0.7e1 /
      0.2e1) - 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 +
      0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i) * x(i + 1) * x(k) ^ 3
      * x(k + 1) + 0.112e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k
      + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) *
      x(i) ^ 6 * sqrt(%pi) * x(k) + 0.112e3 * b ^ 4 * exp(b * (x(k) ^ 2
      + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) -
      x(k) * sqrt(b)) * x(i) ^ 6 * sqrt(%pi) * x(k + 1) - 0.56e2 * b ^
      (0.5e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) +
      x(i + 1) ^ 2)) * x(i) ^ 2 * x(k + 1) * x(k) - 0.84e2 * exp(b *
      (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(k) ^ 3 *
      x(i) * b ^ (0.5e1 / 0.2e1) - 0.126e3 * exp(b * (x(k + 1) ^ 2 +
      0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i + 1) * x(i) * b ^
      (0.3e1 / 0.2e1) + 0.14e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) *
      x(k) + x(i + 1) ^ 2)) * x(k) * x(i + 1) * b ^ (0.3e1 / 0.2e1) +
      0.210e3 * x(k + 1) * x(k) * b ^ 2 * exp(b * (x(k) ^ 2 + x(i) ^ 2 +
      x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i + 1) * sqrt(%pi) * erf(x(i +
      1) * sqrt(b) - x(k) * sqrt(b)) + 0.168e3 * x(k + 1) * x(k) * b ^ 4
      * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * sqrt(%pi) * x(i) ^ 5 +
      0.168e3 * x(k + 1) * x(k) * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 +
      x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) *
      sqrt(b)) * x(i + 1) ^ 5 * sqrt(%pi) - 0.210e3 * x(k + 1) * x(k) *
      b ^ 2 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * x(i + 1) * sqrt(%pi) * erf(x(i) * sqrt(b) - x(k) * sqrt(b))
      + 0.210e3 * x(k + 1) * x(k) * b ^ 2 * exp(b * (x(k) ^ 2 + x(i) ^ 2
      + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i) * sqrt(%pi) * erf(x(i + 1)
      * sqrt(b) - x(k) * sqrt(b)) - 0.168e3 * x(k + 1) * x(k) * b ^ 4 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i + 1) ^ 5 *
      sqrt(%pi) + 0.280e3 * x(k + 1) * x(k) * b ^ 4 * exp(b * (x(k) ^ 2
      + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) -
      x(k) * sqrt(b)) * x(i) ^ 4 * x(i + 1) * sqrt(%pi) + 0.140e3 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(k
      + 1) ^ 4 * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * sqrt(%pi) *
      x(i) * b ^ 3 + 0.252e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^
      2 + x(i + 1) ^ 2)) * x(k) ^ 5 * erf(x(i) * sqrt(b) - x(k) *
      sqrt(b)) * sqrt(%pi) * b ^ 3 - 0.56e2 * exp(b * (x(k) ^ 2 + 0.2e1
      * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(k + 1) * x(i + 1) * x(i) ^
      2 * b ^ (0.5e1 / 0.2e1) - 0.212e3 * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(k + 1) ^ 4 * b ^ (0.5e1 /
      0.2e1) - 0.168e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i +
      1) * sqrt(%pi) * x(i) * x(k + 1) ^ 5 * b ^ 4 + 0.40e2 * exp(b *
      (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(k) ^ 2 *
      x(i + 1) ^ 2 * b ^ (0.5e1 / 0.2e1) - 0.84e2 * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(k + 1) ^ 3 * x(i + 1)
      * b ^ (0.5e1 / 0.2e1) + 0.48e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 *
      x(i + 1) * x(k) + x(i) ^ 2)) * x(k) ^ 3 * x(i + 1) * b ^ (0.5e1 /
      0.2e1) + 0.112e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i +
      1) * sqrt(%pi) * x(k) ^ 6 * b ^ 4 - 0.112e3 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) -
      x(k) * sqrt(b)) * x(i + 1) * sqrt(%pi) * x(k) ^ 6 * b ^ 4 - 0.56e2
      * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) *
      x(k) * x(i) * x(i + 1) ^ 2 * b ^ (0.5e1 / 0.2e1) + 0.252e3 * exp(b
      * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i)
      * sqrt(b) - x(k + 1) * sqrt(b)) * sqrt(%pi) * x(i) ^ 5 * b ^ 3 -
      0.105e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * x(i + 1) * sqrt(%pi) * erf(x(i) * sqrt(b) - x(k) * sqrt(b))
      * b - 0.212e3 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i
      + 1) ^ 2)) * x(i) ^ 4 * b ^ (0.5e1 / 0.2e1) + 0.212e3 * exp(b *
      (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^
      4 * b ^ (0.5e1 / 0.2e1) + 0.14e2 * b ^ (0.3e1 / 0.2e1) * exp(b *
      (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) *
      x(k) - 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i) ^ 5 * x(k) - 0.280e3 * b ^
      3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * sqrt(%pi) * x(i) * x(k
      + 1) ^ 3 * x(k) + 0.420e3 * b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 +
      x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) *
      sqrt(b)) * sqrt(%pi) * x(k + 1) ^ 4 * x(k) - 0.112e3 * b ^ 4 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i) ^ 6 * sqrt(%pi) *
      x(k) + 0.168e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i) ^ 4 * x(i + 1) * x(k) -
      0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i) *
      sqrt(%pi) * x(k + 1) ^ 5 * x(k) - 0.144e3 * exp(b * (x(k + 1) ^ 2
      + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * sqrt(b) + 0.80e2 * exp(b
      * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(k) ^ 6
      * b ^ (0.7e1 / 0.2e1) - 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b *
      (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i) * x(i
      + 1) ^ 3 * x(k) * x(k + 1) - 0.105e3 * b * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * sqrt(%pi) * erf(x(i) *
      sqrt(b) - x(k + 1) * sqrt(b)) * x(k) + 0.105e3 * b * exp(b * (x(k)
      ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * sqrt(%pi) *
      erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(k) - 0.280e3 * x(k + 1) *
      x(k) * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i +
      1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i + 1) ^ 4
      * x(i) * sqrt(%pi) + 0.210e3 * x(k + 1) * x(k) * b ^ 2 * exp(b *
      (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i + 1) *
      sqrt(%pi) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) + 0.56e2 *
      exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) *
      x(k) ^ 2 * x(i + 1) * x(i) * b ^ (0.5e1 / 0.2e1) + 0.80e2 * exp(b
      * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i +
      1) * sqrt(b) - x(k + 1) * sqrt(b)) * sqrt(%pi) * x(k + 1) ^ 7 * b
      ^ 4 + 0.280e3 * b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^
      2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) *
      x(i + 1) ^ 3 * x(i) * sqrt(%pi) * x(k + 1) - 0.210e3 * b ^ 2 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i
      + 1) * x(i) * sqrt(%pi) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b))
      * x(k) - 0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1)
      ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) *
      x(i + 1) * sqrt(%pi) * x(k) ^ 5 * x(k + 1) + 0.105e3 * b * exp(b *
      (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * sqrt(%pi) *
      erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b)) * x(k + 1) - 0.56e2 *
      b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k +
      1) + x(i) ^ 2)) * x(i + 1) ^ 3 * x(k + 1) ^ 2 * x(k) + 0.210e3 * b
      ^ 2 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2))
      * x(i + 1) * x(i) * sqrt(%pi) * erf(x(i) * sqrt(b) - x(k) *
      sqrt(b)) * x(k) - 0.280e3 * x(k + 1) * x(k) * b ^ 4 * exp(b *
      (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) *
      sqrt(b) - x(k + 1) * sqrt(b)) * x(i) ^ 4 * x(i + 1) * sqrt(%pi) -
      0.280e3 * b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i) *
      sqrt(%pi) * x(k) ^ 3 * x(k + 1) + 0.280e3 * b ^ 3 * exp(b * (x(k)
      ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) *
      sqrt(b) - x(k) * sqrt(b)) * x(i + 1) * sqrt(%pi) * x(k) ^ 3 * x(k
      + 1) + 0.280e3 * b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^
      2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b))
      * x(i + 1) * sqrt(%pi) * x(k + 1) ^ 3 * x(k) - 0.56e2 * b ^ (0.7e1
      / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i)
      ^ 2)) * x(i + 1) * x(k + 1) ^ 4 * x(k) - 0.56e2 * b ^ (0.7e1 /
      0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^
      2)) * x(i) ^ 2 * x(k) ^ 3 * x(k + 1) - 0.224e3 * b ^ (0.3e1 /
      0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^
      2)) * x(i) * x(k + 1) - 0.168e3 * x(k + 1) * x(k) * b ^ 4 * exp(b
      * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i)
      * sqrt(b) - x(k) * sqrt(b)) * x(i) ^ 5 * sqrt(%pi) - 0.210e3 * x(k
      + 1) * x(k) * b ^ 2 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2
      + x(i + 1) ^ 2)) * x(i + 1) * sqrt(%pi) * erf(x(i + 1) * sqrt(b) -
      x(k + 1) * sqrt(b)) + 0.210e3 * x(k + 1) * x(k) * b ^ 2 * exp(b *
      (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i) *
      sqrt(%pi) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) - 0.56e2 *
      exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i
      + 1) * x(i) ^ 3 * x(k) ^ 2 * b ^ (0.7e1 / 0.2e1) - 0.212e3 * exp(b
      * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1)
      ^ 4 * b ^ (0.5e1 / 0.2e1) + 0.105e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2
      + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i + 1) * sqrt(%pi) * erf(x(i +
      1) * sqrt(b) - x(k) * sqrt(b)) * b - 0.364e3 * exp(b * (x(k + 1) ^
      2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i) * x(i + 1) ^ 3 * b
      ^ (0.5e1 / 0.2e1) + 0.48e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i +
      1) * x(k) + x(i) ^ 2)) * x(k) * x(i + 1) ^ 3 * b ^ (0.5e1 / 0.2e1)
      - 0.48e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i)
      ^ 2)) * x(k + 1) * x(i + 1) ^ 3 * b ^ (0.5e1 / 0.2e1) - 0.105e3 *
      b * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      sqrt(%pi) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(k + 1) +
      0.196e3 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 *
      x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 2 * x(i) * x(k + 1) +
      0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i) ^
      5 * x(i + 1) * sqrt(%pi) * x(k + 1) - 0.280e3 * b ^ 3 * exp(b *
      (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i + 1) ^
      3 * x(i) * sqrt(%pi) * erf(x(i + 1) * sqrt(b) - x(k + 1) *
      sqrt(b)) * x(k) + 0.56e2 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2
      + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i) ^ 4 * x(k + 1) *
      x(k) - 0.56e2 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) ^ 2 * x(k + 1) ^ 3 *
      x(k) + 0.112e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i
      + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i) *
      sqrt(%pi) * x(k) ^ 6 * b ^ 4 - 0.112e3 * exp(b * (x(k) ^ 2 + x(i)
      ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) -
      x(k) * sqrt(b)) * x(i + 1) ^ 6 * x(i) * sqrt(%pi) * b ^ 4 -
      0.140e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * x(k) ^ 4 * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(i) *
      sqrt(%pi) * b ^ 3 - 0.80e2 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k +
      1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b))
      * sqrt(%pi) * x(k + 1) ^ 7 * b ^ 4 - 0.252e3 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b)
      - x(k + 1) * sqrt(b)) * x(i + 1) ^ 5 * sqrt(%pi) * b ^ 3 + 0.420e3
      * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(i) ^ 4 * x(i + 1) *
      sqrt(%pi) * b ^ 3 + 0.168e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k +
      1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k + 1) *
      sqrt(b)) * x(i + 1) * x(i) * sqrt(%pi) * x(k + 1) ^ 5 * b ^ 4 -
      0.112e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i) *
      sqrt(%pi) * x(k + 1) ^ 6 * b ^ 4 - 0.80e2 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b)
      - x(k + 1) * sqrt(b)) * x(i + 1) ^ 7 * sqrt(%pi) * b ^ 4 - 0.140e3
      * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      x(k + 1) ^ 4 * erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i)
      * sqrt(%pi) * b ^ 3 - 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k
      + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 5 *
      x(k + 1) + 0.112e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i +
      1) * sqrt(%pi) * x(k + 1) ^ 6 * b ^ 4 + 0.105e3 * exp(b * (x(k) ^
      2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i) * sqrt(%pi) *
      erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * b - 0.56e2 * exp(b *
      (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i + 1) *
      x(i) ^ 2 * x(k) ^ 3 * b ^ (0.7e1 / 0.2e1) - 0.84e2 * b ^ (0.5e1 /
      0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^
      2)) * x(i) ^ 3 * x(k) - 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b *
      (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(k + 1) ^
      5 * x(k) - 0.224e3 * b ^ (0.3e1 / 0.2e1) * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(k) +
      0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i +
      1) * x(k + 1) + x(i) ^ 2)) * x(k + 1) ^ 5 * x(k) - 0.112e3 * b ^
      (0.7e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k)
      + x(i) ^ 2)) * x(k) ^ 5 * x(k + 1) - 0.126e3 * b ^ (0.3e1 / 0.2e1)
      * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) *
      x(k) * x(k + 1) - 0.112e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k +
      1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b))
      * x(i) ^ 6 * x(i + 1) * sqrt(%pi) * b ^ 4 - 0.32e2 * exp(b * (x(k
      + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 4 *
      x(k) ^ 2 * b ^ (0.7e1 / 0.2e1) - 0.32e2 * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i) ^ 4 * x(k + 1) ^ 2
      * b ^ (0.7e1 / 0.2e1) + 0.32e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i +
      1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) ^ 3 * x(k + 1) ^ 3 * b ^
      (0.7e1 / 0.2e1) - 0.48e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) *
      x(k) + x(i + 1) ^ 2)) * x(k) * x(i) ^ 3 * b ^ (0.5e1 / 0.2e1) +
      0.84e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^
      2)) * x(k + 1) ^ 3 * x(i) * b ^ (0.5e1 / 0.2e1) + 0.48e2 * exp(b *
      (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(k + 1) ^
      3 * x(i) * b ^ (0.5e1 / 0.2e1) + 0.80e2 * exp(b * (x(k) ^ 2 + x(i)
      ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k +
      1) * sqrt(b)) * x(i) ^ 7 * sqrt(%pi) * b ^ 4 - 0.56e2 * exp(b *
      (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) ^
      4 * x(i) * x(k + 1) * b ^ (0.7e1 / 0.2e1) + 0.56e2 * exp(b * (x(k)
      ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(i) *
      x(k + 1) ^ 4 * b ^ (0.7e1 / 0.2e1) - 0.56e2 * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i) * x(i + 1) ^ 2 *
      x(k + 1) ^ 3 * b ^ (0.7e1 / 0.2e1) - 0.112e3 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b)
      - x(k + 1) * sqrt(b)) * sqrt(%pi) * x(i + 1) * x(k + 1) ^ 6 * b ^
      4 + 0.56e2 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i) ^ 2 * x(k + 1) ^ 3 * x(k)
      - 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i)
      * x(k + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(i) ^ 3 * x(k + 1) *
      x(k) - 0.112e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^
      2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b))
      * sqrt(%pi) * x(k + 1) ^ 6 * x(k) - 0.140e3 * b ^ 3 * exp(b *
      (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i +
      1) * sqrt(b) - x(k) * sqrt(b)) * x(i + 1) ^ 4 * sqrt(%pi) * x(k +
      1) + 0.56e2 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1
      * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 2 * x(k) ^ 3 * x(k +
      1) + 0.56e2 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1
      * x(i) * x(k) + x(i + 1) ^ 2)) * x(i) * x(k) ^ 2 * x(k + 1) -
      0.196e3 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 *
      x(i) * x(k) + x(i + 1) ^ 2)) * x(i + 1) * x(i) ^ 2 * x(k + 1) -
      0.112e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 +
      x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i +
      1) ^ 6 * sqrt(%pi) * x(k) - 0.56e2 * exp(b * (x(k) ^ 2 + 0.2e1 *
      x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(k + 1) ^ 2 * x(i + 1) * x(i)
      * b ^ (0.5e1 / 0.2e1) + 0.112e3 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i
      + 1) * x(k + 1) + x(i) ^ 2)) * x(i) * x(k + 1) ^ 5 * b ^ (0.7e1 /
      0.2e1) + 0.112e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1)
      ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) *
      sqrt(%pi) * x(k + 1) ^ 6 * x(k) + 0.56e2 * b ^ (0.7e1 / 0.2e1) *
      exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i
      + 1) * x(k) ^ 4 * x(k + 1) + 0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) -
      x(k) * sqrt(b)) * x(i + 1) * sqrt(%pi) * x(k) ^ 5 * x(k + 1) +
      0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1 *
      x(i) * x(k) + x(i + 1) ^ 2)) * x(i + 1) * x(i) * x(k) ^ 3 * x(k +
      1) - 0.210e3 * x(k + 1) * x(k) * b ^ 2 * exp(b * (x(k) ^ 2 + x(i)
      ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * x(i) * sqrt(%pi) * erf(x(i +
      1) * sqrt(b) - x(k + 1) * sqrt(b)) - 0.168e3 * exp(b * (x(k) ^ 2 +
      x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b)
      - x(k) * sqrt(b)) * x(i + 1) * x(i) * sqrt(%pi) * x(k) ^ 5 * b ^ 4
      + 0.32e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1)
      ^ 2)) * x(i) ^ 3 * x(k) ^ 3 * b ^ (0.7e1 / 0.2e1) - 0.56e2 * b ^
      (0.5e1 / 0.2e1) * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) +
      x(i + 1) ^ 2)) * x(i) * x(k + 1) ^ 2 * x(k) - 0.280e3 * b ^ 4 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(i + 1) * x(i) * sqrt(%pi)
      * x(k) ^ 4 * x(k + 1) - 0.168e3 * b ^ (0.7e1 / 0.2e1) * exp(b *
      (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i) ^ 4 *
      x(i + 1) * x(k + 1) - 0.112e3 * exp(b * (x(k + 1) ^ 2 + 0.2e1 *
      x(i + 1) * x(k) + x(i) ^ 2)) * x(i) * x(k) ^ 5 * b ^ (0.7e1 /
      0.2e1) - 0.112e3 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) +
      x(i + 1) ^ 2)) * x(i + 1) * x(i) ^ 5 * b ^ (0.7e1 / 0.2e1) -
      0.32e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k + 1) + x(i + 1) ^
      2)) * x(i) * x(k + 1) ^ 5 * b ^ (0.7e1 / 0.2e1) - 0.32e2 * exp(b *
      (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^
      2 * x(k) ^ 4 * b ^ (0.7e1 / 0.2e1) + 0.56e2 * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(i) ^ 2 *
      x(k + 1) ^ 3 * b ^ (0.7e1 / 0.2e1) + 0.80e2 * exp(b * (x(k + 1) ^
      2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 6 * b ^
      (0.7e1 / 0.2e1) + 0.56e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i) * x(k
      + 1) + x(i + 1) ^ 2)) * x(i + 1) * x(i) ^ 4 * x(k + 1) * b ^
      (0.7e1 / 0.2e1) + 0.280e3 * b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 +
      x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k + 1)
      * sqrt(b)) * x(i) * sqrt(%pi) * x(k + 1) ^ 3 * x(k) - 0.80e2 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * sqrt(%pi) * x(k) ^ 7 *
      b ^ 4 + 0.168e3 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 +
      0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) * x(i) * x(k) *
      x(k + 1) - 0.168e3 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) * x(i) * x(k +
      1) * x(k) + 0.168e3 * b ^ (0.7e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2
      + 0.2e1 * x(i + 1) * x(k) + x(i) ^ 2)) * x(i) * x(k) ^ 4 * x(k +
      1) - 0.56e2 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k + 1) ^ 2 + 0.2e1
      * x(i + 1) * x(k) + x(i) ^ 2)) * x(i + 1) ^ 2 * x(k) * x(k + 1) -
      0.112e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * sqrt(%pi) * x(i) *
      x(k) ^ 6 * b ^ 4 - 0.196e3 * b ^ (0.5e1 / 0.2e1) * exp(b * (x(k) ^
      2 + 0.2e1 * x(i + 1) * x(k + 1) + x(i) ^ 2)) * x(i + 1) ^ 2 * x(i)
      * x(k) + 0.140e3 * b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1)
      ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k + 1) *
      sqrt(b)) * x(i + 1) ^ 4 * sqrt(%pi) * x(k) + 0.168e3 * b ^ 4 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b)) * x(i + 1) ^ 5 * x(i) *
      sqrt(%pi) * x(k) + 0.168e3 * b ^ 4 * exp(b * (x(k) ^ 2 + x(i) ^ 2
      + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k) *
      sqrt(b)) * x(i) * sqrt(%pi) * x(k) ^ 5 * x(k + 1) - 0.140e3 * b ^
      3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b)) * sqrt(%pi) * x(i) ^ 4 *
      x(k) + 0.80e2 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i
      + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * sqrt(%pi) *
      x(k) ^ 7 * b ^ 4 + 0.80e2 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k +
      1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) * sqrt(b))
      * x(i + 1) ^ 7 * sqrt(%pi) * b ^ 4 - 0.32e2 * exp(b * (x(k) ^ 2 +
      0.2e1 * x(i) * x(k + 1) + x(i + 1) ^ 2)) * x(i) ^ 2 * x(k + 1) ^ 4
      * b ^ (0.7e1 / 0.2e1) - 0.32e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i)
      * x(k + 1) + x(i + 1) ^ 2)) * x(i) ^ 3 * x(k + 1) ^ 3 * b ^ (0.7e1
      / 0.2e1) + 0.32e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1) * x(k +
      1) + x(i) ^ 2)) * x(i + 1) ^ 4 * x(k + 1) ^ 2 * b ^ (0.7e1 /
      0.2e1) + 0.32e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) +
      x(i + 1) ^ 2)) * x(i) ^ 2 * x(k) ^ 4 * b ^ (0.7e1 / 0.2e1) +
      0.56e2 * exp(b * (x(k + 1) ^ 2 + 0.2e1 * x(i + 1) * x(k) + x(i) ^
      2)) * x(i + 1) ^ 2 * x(i) * x(k) ^ 3 * b ^ (0.7e1 / 0.2e1) +
      0.112e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^
      2)) * erf(x(i) * sqrt(b) - x(k) * sqrt(b)) * x(i) ^ 6 * x(i + 1) *
      sqrt(%pi) * b ^ 4 - 0.56e2 * exp(b * (x(k) ^ 2 + 0.2e1 * x(i + 1)
      * x(k + 1) + x(i) ^ 2)) * x(i + 1) * x(i) * x(k + 1) ^ 4 * b ^
      (0.7e1 / 0.2e1) + 0.112e3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k +
      1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i) * sqrt(b) - x(k + 1) * sqrt(b))
      * sqrt(%pi) * x(i) * x(k + 1) ^ 6 * b ^ 4 + 0.112e3 * b ^ 4 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      erf(x(i + 1) * sqrt(b) - x(k + 1) * sqrt(b)) * x(i + 1) ^ 6 *
      sqrt(%pi) * x(k + 1) + 0.112e3 * b ^ (0.7e1 / 0.2e1) * exp(b *
      (x(k + 1) ^ 2 + 0.2e1 * x(i) * x(k) + x(i + 1) ^ 2)) * x(i + 1) *
      x(i) ^ 2 * x(k) ^ 2 * x(k + 1) + 0.280e3 * b ^ 4 * exp(b * (x(k) ^
      2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) *
      sqrt(b) - x(k) * sqrt(b)) * x(i + 1) * x(i) * sqrt(%pi) * x(k) ^ 4
      * x(k + 1) - 0.140e3 * b ^ 3 * exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k
      + 1) ^ 2 + x(i + 1) ^ 2)) * erf(x(i + 1) * sqrt(b) - x(k) *
      sqrt(b)) * x(i + 1) ^ 4 * sqrt(%pi) * x(k) - 0.280e3 * b ^ 3 *
      exp(b * (x(k) ^ 2 + x(i) ^ 2 + x(k + 1) ^ 2 + x(i + 1) ^ 2)) *
      x(i) ^ 3 * x(i + 1) * sqrt(%pi) * erf(x(i) * sqrt(b) - x(k) *
      sqrt(b)) * x(k + 1)) * b ^ (-0.9e1 / 0.2e1) / (x(i + 1) ^ 2 * x(k
      + 1) ^ 2 - 0.2e1 * x(i + 1) ^ 2 * x(k + 1) * x(k) + x(i + 1) ^ 2 *
      x(k) ^ 2 - 0.2e1 * x(i + 1) * x(k + 1) ^ 2 * x(i) + 0.4e1 * x(i +
      1) * x(i) * x(k + 1) * x(k) - 0.2e1 * x(i + 1) * x(i) * x(k) ^ 2 +
      x(i) ^ 2 * x(k + 1) ^ 2 - 0.2e1 * x(i) ^ 2 * x(k + 1) * x(k) +
      x(i) ^ 2 * x(k) ^ 2) / 0.6720e4;<br>
      <br>
      <div class="gmail_quote">2011/10/3 Mike Page <span dir="ltr"><<a
            moz-do-not-send="true"
            href="mailto:Mike@page-one.waitrose.com">Mike@page-one.waitrose.com</a>></span><br>
        <blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
          0.8ex; border-left: 1px solid rgb(204, 204, 204);
          padding-left: 1ex;">
          <div>
            <div><span><font color="#0000ff" face="Arial" size="2">Can
                  you post some code that shows the problem?</font></span></div>
            <div><span></span> </div>
            <div><span><font color="#0000ff" face="Arial" size="2">Sounds
                  like maybe you are creating a string instead of a
                  numeric matrix.</font></span></div>
            <div><span></span> </div>
            <div><span><font color="#0000ff" face="Arial" size="2">Regards,</font></span></div>
            <div><span><font color="#0000ff" face="Arial" size="2">Mike.</font></span></div>
            <font color="#888888">
              <div><span></span> </div>
            </font>
            <blockquote>
              <div>
                <div class="h5">
                  <div dir="ltr" align="left"><font face="Tahoma"
                      size="2">-----Original Message-----<br>
                      <b>From:</b> <a moz-do-not-send="true"
                        href="http://severine.pl" target="_blank">severine.pl</a>
                      [mailto:<a moz-do-not-send="true"
                        href="mailto:severine.pl@gmail.com"
                        target="_blank">severine.pl@gmail.com</a>]<br>
                      <b>Sent:</b> 03 October 2011 14:10<br>
                      <b>To:</b> <a moz-do-not-send="true"
                        href="mailto:users@lists.scilab.org"
                        target="_blank">users@lists.scilab.org</a><br>
                      <b>Subject:</b> [scilab-Users] too large string<br>
                      <br>
                    </font></div>
                  Hi!<br>
                  <br>
                  I woulf like to find an answer to my problem. <br>
                  <br>
                  I'm doing very big calculus in scilab, and calculating
                  very big matrix.<br>
                  And when i'm trying to executing the programm, Scilab
                  says me:<br>
                  <br>
                  "Too large string"<br>
                  <br>
                  What must I do?<br>
                  <br>
                  Séverine Paul<br>
                  <br>
                </div>
              </div>
              <div class="im">
                <hr align="left" width="300"> View this message in
                context: <a moz-do-not-send="true"
href="http://mailinglists.scilab.org/too-large-string-tp3389713p3389713.html"
                  target="_blank">too large string</a><br>
                Sent from the <a moz-do-not-send="true"
href="http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html"
                  target="_blank">Scilab users - Mailing Lists Archives
                  mailing list archive</a> at Nabble.com.<br>
              </div>
            </blockquote>
          </div>
        </blockquote>
      </div>
      <br>
    </blockquote>
    <br>
  </body>
</html>