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 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;">
<u></u>
<div>
<div><span><font color="#0000ff" face="Arial" size="2">Can
you post some code that shows the problem?</font></span></div>
<div><span><font color="#0000ff" face="Arial" size="2"></font></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><font color="#0000ff" face="Arial" size="2"></font></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><font color="#0000ff" face="Arial" size="2"></font></span> </div>
</font><blockquote><div><div></div><div class="h5">
<div dir="ltr" align="left"><font face="Tahoma" size="2">-----Original Message-----<br><b>From:</b> <a href="http://severine.pl" target="_blank">severine.pl</a>
[mailto:<a 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 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 width="300" align="left">
View this message in context: <a href="http://mailinglists.scilab.org/too-large-string-tp3389713p3389713.html" target="_blank">too
large string</a><br>Sent from the <a 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>