Elastik qovushqoq plastinkaning kuchlanish-deformatsion holatini matematik modellashtirish

nuqta[3].Y = Mkk - (float)(MK[MN[i].i - 1].y * koef);

graphics.DrawPolygon(myPen, nuqta);
}
for (i = 1; i <= n; i++)
{
x = (float)((MK[i - 1].x) * koef) + xr;
y = Mkk - (float)((MK[i - 1].y) * koef);
graphics.DrawString(i.ToString(), font, myBrush, x, y);
}*/

}
double Li(int e,int Ai )
{
double xi, xj, xk, x, yi, yj, yk, y, s1, s2, s3, a1, a2, a3, st;
xi = MK[MN[e].i - 1].x; yi = MK[MN[e].i - 1].y;
xj = MK[MN[e].j - 1].x; yj = MK[MN[e].j - 1].y;
xk = MK[MN[e].k - 1].x; yk = MK[MN[e].k - 1].y;
x = (xi + xj + xk) / 3; y = (yi + yj + yk) / 3;
st = Math.Abs(xi * yj + xj * yk + xk * yi - xi * yk - xj * yi - xk * yj) / 2;
s1 = Math.Abs(x * yj + xj * yk + xk * y - x * yk - xj * y - xk * yj) / 2;
s2 = Math.Abs(xi * y + x * yk + xk * yi - xi * yk - x * yi - xk * y) / 2;
s3 = Math.Abs(xi * yj + xj * y + x * yi - xi * y - xj * yi - x * yj) / 2;
a1 = s1 / st; a2 = s2 / st; a3 = s3 / st;
if (Ai == 1) return a1;
else if(Ai == 2) return a2;
else if(Ai == 3) return a3;

return 1;
}
public void k_e(int e, bool pech)
{
double xx, yy,Ni,Nj,Nk,r,Ri,Rj, V,rt,zt;
double[] b = new double[3], c = new double[3], fe = new double[6], fh = new double[6], a=new double[3];
a[0] = MK[MN[e].j - 1].x * MK[MN[e].k - 1].y - MK[MN[e].k - 1].x * MK[MN[e].j - 1].y;
a[1] = MK[MN[e].k - 1].x * MK[MN[e].i - 1].y - MK[MN[e].i - 1].x * MK[MN[e].k - 1].y;
a[2] = MK[MN[e].i - 1].x * MK[MN[e].j - 1].y - MK[MN[e].j - 1].x * MK[MN[e].i - 1].y;
b[0] = MK[MN[e].j - 1].y - MK[MN[e].k - 1].y;
b[1] = MK[MN[e].k - 1].y - MK[MN[e].i - 1].y;
b[2] = MK[MN[e].i - 1].y - MK[MN[e].j - 1].y;