clc
clear all close all
Dim = 50;  %максимальная размерность матриц
t_prod = 100 /   (2*4^3-4^2);  %время вычисления
t_send = t_prod/5;  %время передачи матрицы
for n = 1:Dim
%%%Parallel algorithm  
N2=2*n^3-n^2;
T2=0;%Т - время работы алгоритма с простаиванием.  Т2 - без простаивания  
T = t_send*(n^2)*2;  %Передача А1 и А2 в Т1  
T = T + N2*t_prod;
T = T + t_send*(n^2);  %Передача A3 в Т1  
T = T + N2*t_prod;
T = T + t_send*(n^2);  %Передача Т2 в Т1
T2=T+t_send*((n^2)/2);%Передача Т1/2 в Т2
T2 = T2 +   (N2*t_prod)/2;  %Умножение Т1/2 на Т2  
T2 = T2 + 2*t_send*((n^2)/2);  %вывод матриц  
T = T + N2*t_prod;  %умножение Т1 на Т2  
T = T + t_send*(n^2);  %вывод
%Ts - время работы последовательного алгоритма  
Ts = t_send*(n^2)*7;  
Ts = Ts + t_prod*N2*5;  
e(n)  = Ts/T;  
f(n)=Ts/T2;
end;
plot(1:Dim,e,'r','LineWidth',2)  
hold on
plot(1:Dim,f,'b','LineWidth',2);
grid on;  
hold off;  
% pause;
figure;
for n = 1:Dim %аналогичный процесс для трех транспьюетров.
N2=2*n^3-n^2; T2=0;
T = t_send*(n^2)*2;
T = T + N2*t_prod;
T = T + t_send*(n^2) ;
T3=T+t_send*((n^2)/2) ;
T3 - T3 +   (N2*t_prod)/2;
T3 = T3 + 2*t_send*((n^2)/2);
T = T + N2*t_prod;
T = T + t_send*(n^2) ;
T4=T3+2*t_send*((n^2)/2) ;
T4 = T4 +   (N2*t_prod)/2;
T4 = T4 + 2*t_send*((n^2)/2);
T2=T+3*t_send*((n^2)/3);
T2 = T2 +   (N2*t_prod)/3;
T2 = T2 + 3*t_send*((n^2)/3);
T3=T3+2*t_send*((n^2)/3)+ t_send*(n^2);
T3=T3 +   (N2*t_prod)/3;
T3 = T3 + 3*t_send*((n^2)/3);
T = T + N2*t_prod;
T = T + t_send*(n^2);
%%%Series algorithm
Ts = t_send*(n^2)*7;
Ts = Ts + t_prod*N2*5;
j(n)  = Ts/T;
l(n)=Ts/T2;
m(n)=Ts/T3;
k(n)=Ts/T4;
end;
plot(1:Dim,j,'r','LineWidth',2)  
hold on
plot(1:Dim,m,'b','LineWidth',2)  
plot(1:Dim,k,'y','LineWidth',2)  
plot(1:Dim,1,'g','LineWidth',2);
grid on;