%nezname q1 a q2
syms q1 q2
%nákladové funkce obou firem
c1 = input('Zadejte náklady firmy A na jeden kus produktu (Kč):')
c2 = input('Zadejte náklady firmy B na jeden kus produktu (Kč):')
%určení odhadovaných cifer tržní ceny
a = 10^(numel(num2str(c1)))
b = 10^(numel(num2str(c2)))
if a>=b
%trzni cena vyrobku
p = 10^(numel(num2str(c1)))-(q1+q2)
else
p = 10^(numel(num2str(c2)))-(q1+q2)
end
%zisk firem
z1 = p*q1-c1*q1
z2 = p*q2-c2*q2
%derivace
firmaA = diff(z1,q1)
firmaB = diff(z2,q2)
%soustava rovnic pro nalezení velikosti produkce obou firem
eqns = [solve(firmaA,q1) == q1, solve(firmaB,q2) == q2];
soustavaRovnic = solve(eqns,[q1,q2])
%hodnota produkce
produkceA = vpa(soustavaRovnic.q1,4)
produkceB = vpa(soustavaRovnic.q2,4)
%tržní cena produktu
if a>=b
trzniCena = vpa(10^(numel(num2str(c1)))-(produkceA+produkceB),4)
else
trzniCena = vpa(10^(numel(num2str(c2)))-(produkceA+produkceB),4)
end
%zisk firem za stanovenou produkci
ziskA = vpa(trzniCena*produkceA-(c1*produkceA),6)
ziskB = vpa(trzniCena*produkceB-(c2*produkceB),6)
ziskZaJedenKusA = vpa(ziskA/produkceA,4)
ziskZaJedenKusB = vpa(ziskB/produkceB,4)
%vysledky setreni
TrzniCena = ['Rovnovazna cena produktu je:' int2str(trzniCena) ' Kč']
ZiskFirem = ['Zisk firmy A z jednoho kusu produktu bude:' int2str(ziskZaJedenKusA) ' Kč' ' a zisk firmy B z jednoho kusu produktu bude:' int2str(ziskZaJedenKusB) ' Kč']