Elastik qovushqoq plastinkaning kuchlanish-deformatsion holatini matematik modellashtirish

}
//SSigma[115] -= 18014; SSigma[64] -= 854;
//SSigma[116] -= 875;
//SSigma[111] -= 6716;
SSigma[118] -= 14171;
SSigma[117] -= 6716;
SSigma[115] -= 18014; SSigma[64] -= 854;


for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++) S[i, j] = 0;
F[i] = 0;
}
///Console.WriteLine("**************************************");
for (int k0 = 0; k0 < m; k0++)
{
mi[0] = MN[k0].i - 1; mi[1] = MN[k0].j - 1; mi[2] = MN[k0].k - 1;

k[0, 0] = A[k0] * 2.0 / 12.0; k[0, 1] = A[k0] * 1.0 / 12.0; k[0, 2] = A[k0] * 1.0 / 12.0;
k[1, 0] = A[k0] * 1.0 / 12.0; k[1, 1] = A[k0] * 2.0 / 12.0; k[1, 2] = A[k0] * 1.0 / 12.0;
k[2, 0] = A[k0] * 1.0 / 12.0; k[2, 1] = A[k0] * 1.0 / 12.0; k[2, 2] = A[k0] * 2.0 / 12.0;
f[0] = -A[k0] * sigma[k0, 1] * 1.0 / 3.0;
f[1] = -A[k0] * sigma[k0, 1] * 1.0 / 3.0;
f[2] = -A[k0] * sigma[k0, 1] * 1.0 / 3.0;
//------------------------------------------------------------
for (int i = 0; i < 3; i++)
{
k1 = mi[i];
F[k1] = F[k1] + f[i];
for (int j = 0; j < 3; j++)
{
k2 = mi[j];
S[k1, k2] = S[k1, k2] + k[i, j];
}

}
}
SSigmay = Gaus(2);
listBox1.Items.Add("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-**-*-*-*-*--*-*-*-*");
listBox1.Items.Add("\t\t\t Sigma");
for (int i = 0; i < 115; i++)
{
if (SSigmay[i] > 100) SSigmay[i] -= 700;
else SSigmay[i] -= 300;
}
// Double[] tao=new Double[n];
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++) S[i, j] = 0;
F[i] = 0;