static double Sum(double [] m, int n, double x, double y, double z)
{
double S = 0;
for (int i = p, j = 0; i < k; i++, j++)
{
S += m[j] * L_Pol.SF(y, z, i);
}
return(S);
}
public static double Int3(Function f, double[] b, double[] a, int n, int num) //полином * якобиан (собственные ф-ции и якобиан) * подинтегральная функция
{
double hx = (b[0] - a[0]) / n, hy = (b[1] - a[1]) / n, hz = (b[2] - a[2]) / n, S = 0;
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++)
{
for (int k = 0; k < n; k++)
{
double xi = a[0] + hx / 2 + i * hx;
double yj = a[1] + hy / 2 + j * hy;
double zk = a[2] + hz / 2 + k * hz;
S -= hx * hy * hz * xi * xi * Math.Sin(zk) * L_Pol.koef(yj, zk, num) * f(xi * Math.Sin(zk) * Math.Cos(yj), xi * Math.Sin(zk) * Math.Sin(yj), xi * Math.Cos(zk));
}
}
}
return(S);
}
public static double L(double x, double y, double z)
{
return(L_Pol.SF(y, z, 1));// * L_Pol.SF(y, z, 2);
}