<?php
//ini_set('display_errors',1);
//error_reporting(E_ALL);
{
$iter = 100000; //количество итераций
$Xmax = array(3); //верхние границы
$Xmin = array(1); //нижние границы
$F = 'cos($x[0])'; //подынтегральное выражение
}
$U = count($Xmin); //кратность интеграла
$Xint = array();
$x = array();
$SumH = 0;
for ($i = 0; $i <= $U; $i++)
{
$Xint[$i] = $Xmax[$i] - $Xmin[$i];
}
for ($j = 1; $j <= $iter; $j++)
{
for ($i = 0; $i <= $U; $i++)
{
$x[$i] = $Xmin[$i] + lcg_value() * $Xint[$i];
}
eval ("\$FA = $F;");
$SumH += $FA;
}
$rez = 1;
for ($i = 0; $i <= $U-1; $i++)
{
$rez = $rez*$Xint[$i];
}
$rez = ($SumH*$rez)/ $iter;
echo ("Answer: $rez");
?>