S = S * log (t +1);% a1 = log (t +1)
figure (1)
plot (t, P, 'g');
plot (t, S, 'r');
xlabel ('Time');
ylabel ('Profits, Expense ');
legend ('Profits', 'Expense');
figure (2)
D = P-S;
P_m = 0.25 * p * a2 * ((a1/a2 + N0) ^ 2) % Формула максимального прибутку
plot (t, D, 'm');
xlabel ('Time');
ylabel ('Profits-Expense = 4istaya pribyl ');
grid on;
figure (3)
plot (t, N, 'y *');
xlabel ('Time');
ylabel ('N koli4estvo pokupatelei ');
grid on;
% mod23.m
clear, clf, clc, hold on, grid on;
T = 0.8;
p = 2;% Прибуток від одного продажу
s1 = 20;% Витрати на 1 рекламну акцію
for a1 = [0:5:150];% інтенсивність рекламної компанії
aa1 = a1/50 ;% Змінюється від 0 до 40/50 = 0.8
a2 = 2;% ступінь спілкування потенційних покупців
N0 = 6 + aa1/a2; % Число потенційних покупців (тис.чол)
m = aa1/a2;
[t, N] = ode45 (@ mma1, [0 T], m, [], N0, aa1, a2);
P = P * (N-(aa1/a2));
S = S1 * aa1 * t;
% SumD = 0;
% for i = 1: length (t) % Тут рахуємо сумарний дохід для кожного коеф. a1
D = P-S;
% SumD = SumD + D (i);
% end
% SumD
P_m = 0.25 * p * a2 * ((aa1/a2 + N0) ^ 2) ;% Формула максимального прибутку
plot (aa1, P_m, 'rp');
plot (aa1, D, 'gh');
xlabel ('a1');
ylabel ('P m, D ');
legend ('P m ',' Dohod ')
grid on;
end
% mod231.m
clear, clf, clc, hold on, grid on;
T = 0.8;
p = 2;% Прибуток від одного продажу
s1 = 20;% Витрати на 1 рекламну акцію
for a1 = [0:2:40];% інтенсивність рекламної компанії
aa1 = a1/50;
a2 = 2;% ступінь спілкування потенційних покупців
N0 = 6 + aa1/a2; % Число потенційних покупців (тис.чол)
m = aa1/a2;
[t, N] = ode45 (@ mma1, [0 T], m, [], N0, aa1, a2);
P = P * (N-(aa1/a2));
S = S1 * aa1 * t;
SumD = 0;
for i = 1: length (t)
D = P-S;
SumD = SumD + D (i);
end
SumD;
P_m = 0.25 * p * a2 * ((aa1/a2 + N0) ^ 2) ;% Формула максимального прибутку
plot (aa1, P_m, 'rp');
plot (aa1, D, 'gh');
xlabel ('a1');
ylabel ('P m, D, SumD ');
legend ('P m ',' Dohod ')
plot (aa1, SumD, 'c *');
grid on;
end
% mod232.m
clear, clf, clc, hold on, grid on;
T = 0.8;
p = 2;% Прибуток від одного продажу
s1 = 20;% Витрати на 1 рекламну акцію
a1 = 0.4;
for a2 = [4...