ZADACHA_UPROSHEN_URAVNENIAY

clear
clc
close all
% M1=[24  -4  -10;-4  29  -4;-10  -4  42];
% r1=[6;-18;-12];
%
% Dis1 = M1\r1

E1=6;   E2=10;  E3=14;  J=2;    r1=10;  r2=5;   r3=8;   r4=20;  r5=16;  r6=4;   r7=10;  r8=22;  r9=2;   r10=16;     r11=12;     EJ=8;

Mout=[];

% строки составленных уравнений
str1='  r7*i11+r6*(i11-i22)+r1*(i11-i33)-E1                 ;';
str2='  (r2+r5)*i22+r6*(i22-i33)+r6*(i22-i11)+E2+EJ         ;';
str3='  r1*(i33-i11)+r6*(i33-i22)+(r3+r4)*i33+E1+E3-EJ      ;';

Mout(1,:)= Mupr(str1,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ);
Mout(2,:)= Mupr(str2,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ);
Mout(3,:)= Mupr(str3,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ);

Mout % матрица
Mout(1:3,1:3)\Mout(1:3,4:4) %  решение i11 i22 i33




function[val] = M(str,i11,i22,i33,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ)
% str='r7*i11+r6*(i11-i22)+r1*(i11-i33)-E1;';
val=eval(str);
end

function[val] = Mupr(str,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ)
% str='r7*i11+r6*(i11-i22)+r1*(i11-i33)-E1;';
i11=0;
i22=0;
i33=0;
Chislo=M(str,i11,i22,i33,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ);
i11=1;
i22=0;
i33=0;
i11get=M(str,i11,i22,i33,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ)-Chislo;
i11=0;
i22=1;
i33=0;
i22get=M(str,i11,i22,i33,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ)-Chislo;
i11=0;
i22=0;
i33=1;
i33get=M(str,i11,i22,i33,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ)-Chislo;
Chislo=-Chislo;
Arr=[i11get i22get i33get Chislo];
val=Arr;
end



% Mout(2,:)=
% Mupr(str(2,1:end),r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ);
% Mout
% i11=0;
% i22=0;
% i33=0;
% Chislo=M(i11,i22,i33,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ);
% i11=1;
% i22=0;
% i33=0;
% i11get=M(i11,i22,i33,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ)-Chislo;
% i11=0;
% i22=1;
% i33=0;
% i22get=M(i11,i22,i33,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ)-Chislo;
% i11=0;
% i22=0;
% i33=1;
% i33get=M(i11,i22,i33,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,E1,E2,E3,EJ)-Chislo;
% Chislo=-Chislo;
% Arr=[i11get i22get i33get Chislo]