Elastik qovushqoq plastinkaning kuchlanish-deformatsion holatini matematik modellashtirish

//s2[2] = SSigmay[10];
//s2[3] = SSigmay[15];
//s2[4] = SSigmay[20];
//s2[5] = SSigmay[25];
//s2[6] = SSigmay[30];
//s2[7] = SSigmay[35];
//s2[8] = SSigmay[40];
s2[0] = tao[0];
s2[1] = tao[5];
s2[2] = tao[10];
s2[3] = tao[15];
s2[4] = tao[20];
s2[5] = tao[25];
s2[6] = tao[30];
s2[7] = tao[35];
s2[8] = tao[40];
s1[0,0] = (2 + tao[0] / 10000) * Math.Cos(fi*Math.PI/180);
s1[0, 1] = (2 + tao[0]/10000) * Math.Sin(fi * Math.PI / 180); fi -= fk;
s1[1, 0] = (2+ tao[5] / 10000) *Math.Cos(fi * Math.PI / 180);
s1[1, 1] = (2 + tao[5] / 10000) *Math.Sin(fi * Math.PI / 180); fi -=fk;
s1[2, 0] = (2+ tao[10] / 10000) * Math.Cos(fi * Math.PI / 180);
s1[2, 1] = (2 + tao[10] / 10000) * Math.Sin(fi * Math.PI / 180); fi -= fk; ;
s1[3, 0] = (2+ tao[15] / 10000) * Math.Cos(fi * Math.PI / 180);
s1[3, 1] = (2 + tao[15] / 10000) * Math.Sin(fi * Math.PI / 180); fi -= fk;
s1[4, 0] = (2 + tao[20] / 10000)*Math.Cos(fi * Math.PI / 180);
s1[4, 1] = (2 + tao[20] / 10000) * Math.Sin(fi * Math.PI / 180); fi -= fk;
s1[5, 0] = (2 + tao[25] / 10000) * Math.Cos(fi * Math.PI / 180);
s1[5, 1] = (2 + tao[25] / 10000)*Math.Sin(fi * Math.PI / 180); fi -= fk;
s1[6, 0] = (2 + tao[30] / 10000) * Math.Cos(fi * Math.PI / 180);
s1[6, 1] = (2 + tao[30] / 10000)* Math.Sin(fi * Math.PI / 180); fi -= fk;
s1[7, 0] = (2 + tao[35] / 10000) * Math.Cos(fi * Math.PI / 180);
s1[7, 1] = (2 + tao[35] / 10000)* Math.Sin(fi * Math.PI / 180); fi -= fk;
s1[8, 0] = (2 + tao[40] / 10000) * Math.Cos(fi * Math.PI / 180);
s1[8, 1] = (2 + tao[40] / 10000)*Math.Sin(fi * Math.PI / 180); fi -= fk ;

for (i = 0; i < 9; i++)