:5:100];% інтенсивність рекламної компанії
aa2 = a2/50; % Ступінь спілкування потенційних покупців
N0 = 6 + a1/aa2; % Число потенційних покупців (тис.чол)
m = a1/aa2;
[t, N] = ode45 (@ mma2, [0 T], m, [], N0, a1, aa2);
P = P * (N-(a1/aa2));
S = S1 * a1 * t;
SumD = 0;
for i = 1: length (t)
D = P-S;
SumD = SumD + D (i);
end
SumD;
P_m = 0.25 * p * aa2 * ((a1/aa2 + N0) ^ 2) ;% Формула максимального прибутку
plot (aa2, P_m, 'rp');
plot (aa2, D, 'gh');
xlabel ('a2');
ylabel ('P m, D, SumD ');
legend ('P m ',' Dohod ')
plot (aa2, SumD, 'c *');
grid on;
end
% modp.m
clear, clf, clc, hold on, grid on;
global a1 a2;
T_max = 10;
N0 = 6;% число потенційних покупців (тис.чол)
p = 20;% Прибуток від одного продажу
s = 5;% Витрати на одну рекламну акцію
[t, N] = ode45 (@ p, [0 T_max], 0, [], N0);
P = P * N;
S = S * (sin (2 * t) +0.5 * t);% 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;
% f.m
function dN = f (t, N, N0, a1)
dN = A1 * N0;
% f1.m
function dN = f1 (t, N, N0, a1, a2)
dN = (A1 + a2 * N) * (N0-N);
% f41.m
function dN = f41 (t, N, N0, a1, a2)
global a1 a2;
a1 = 1/(t +1);% інтенсивність рекламної компанії
a2 = Sqrt (t); % Ступінь спілкування потенційних покупців
dN = (A1 + a2 * N) * (N0-N);
% mm.m
function dN = mm (t, N, N0, a1, a2)
dN = A2 * N * (N0-N);
% mma1.m
function dN = mma1 (t, N, N0, aa1, a2)
dN = A2 * N * (N0-N);
% mma2.m
function dN = mma2 (t, N, N0, a1, aa2)
dN = Aa2 * N * (N0-N);
% p.m
function dN = p (t, N, N0, a1, a2)
global a1 a2;
a1 = Cos (2 * t) +0.5; % Інтенсивність рекламної компанії
a2 = Sqrt (t); % Ступінь спілкування потенційних покупців
dN = (a1 + a2 * N) * (N0-N);
/
