clear all; close all;
%
% *****Obr(2) **************** výdajové omezení ************************
%
p1=4;
p2=6;
M0=48;
x1=0:0.01:12;
xM=(M0-p1.*x1)./p2;
%
% *** Obr(2A) ********** lineární užitková funkce **********************
%
hold on;
 uti1=40;
 uti2=60;
 uti3=80;
 alf1=5;
 alf2=6;
%
 lif1=(uti1-alf1.*x1)./alf2;
 lif2=(uti2-alf1.*x1)./alf2;
 lif3=(uti3-alf1.*x1)./alf2;
%
subplot(2,3,1);
gr=plot(x1,lif1,'b',x1,lif2,'b',x1,lif3,'b',x1,xM,'r');
axis([0 12 0 12]);box on;
title('Lineární užitková funkce');

text(1.8,3,'\rmu^0=40','fontsize',10,'color','b','FontName','Tahoma');
text(1.8,7,'\rmu^1=60','fontsize',10,'color','b','FontName','Tahoma');
text(5.6,8,'\rmu^2=80','fontsize',10,'color','b','FontName','Tahoma');

xpos=5.8;
ypos=4.45;
text(xpos,ypos,'\rm4x_1+6x_2 =48','color','r','fontsize',10,'FontName','Tahoma');
set(gca,'XTickMode','Manual');set(gca,'XTick',[0 2 4 6 8 10 12]);
set(gca,'YTickMode','Manual');set(gca,'YTick',[0 2 4 6 8 10 12]);
hold off;
%
% *** Obr(2B) ************* kvadratická užitková funkce ******************
%
hold on;
 uti1=120;
 uti2=224;
 uti3=360;
 alf1=4;
 alf2=6;
%
 kvf1=sqrt(max((uti1-alf1.*x1.^2)./alf2,0));
 kvf2=sqrt(max((uti2-alf1.*x1.^2)./alf2,0));
 kvf3=sqrt(max((uti3-alf1.*x1.^2)./alf2,0));
%
subplot(2,3,2);
gr=plot(x1,kvf1,'b',x1,kvf2,'b',x1,kvf3,'b',x1,xM,'r');
axis([0 12 0 12]);box on;
title('Kvadratická užitková funkce');

text(4,3.5,'\rmu^0=120','fontsize',10,'FontName','Tahoma','color','b');
text(5,5.05,'\rmu^1=224','fontsize',10,'FontName','Tahoma','color','b');
text(6,6.5,'\rmu^2=360','fontsize',10,'FontName','Tahoma','color','b');

text(6,14,'Příklady indiferenčních křivek dvoukomoditních užitkových funkcí','fontweight','bold','fontsize',15,'HorizontalAlignment','center');
xpos=8.0;
ypos=3.0;
text(xpos,ypos,'\rm4x_1+6x_2=48','color','r','fontsize',10,'FontName','Tahoma');
set(gca,'XTickMode','Manual');set(gca,'XTick',[0 2 4 6 8 10 12]);
set(gca,'YTickMode','Manual');set(gca,'YTick',[0 2 4 6 8 10 12]);
hold off;

% *** Obr(2C)*************** odmocninná  užitková funkce ***************
%
hold on;
 uti1=18;
 uti2=22;
 uti3=26;
 alf1=4;
 alf2=6;
%
 odf1=((uti1-alf1.*sqrt(x1))/alf2).^2;
 odf2=((uti2-alf1.*sqrt(x1))/alf2).^2;
 odf3=((uti3-alf1.*sqrt(x1))/alf2).^2;
%
subplot(2,3,3);
gr=plot(x1,odf1,'b',x1,odf2,'b',x1,odf3,'b',x1,xM,'r');
axis([0 12 0 12]);box on;hold off;
title('odmocninná užitková funkce');

text(3.5,3.55,'\rmu^0=18','fontsize',10,'FontName','Tahoma','Color','b');
text(4.55,5.5,'\rmu^1=22','fontsize',10,'FontName','Tahoma','Color','b');
text(5.6,8.05,'\rmu^2=26','fontsize',10,'FontName','Tahoma','Color','b');

xpos=6.6;
ypos=3.95;
text(xpos,ypos,'\rm4x_1+6x_2=48','color','r','fontsize',10,'FontName','Tahoma');
set(gca,'XTickMode','Manual');set(gca,'XTick',[0 2 4 6 8 10 12]);
set(gca,'YTickMode','Manual');set(gca,'YTick',[0 2 4 6 8 10 12]);
hold off;
%
% *** Obr(2D) ************* Leontiefova užitková funkce *****************
%
hold on;
 uti1=12;
 uti2=22;
 uti3=32;
 alf1=5;
 alf2=6;
 lef1=uti1./alf1;
 lef2=uti2./alf1;
 lef3=uti3./alf1;
 arrx=(alf1.*x1)./alf2;
 subplot(2,3,4);
 gr=plot(x1,lef1,'b',x1,lef2,'b',x1,lef3,'b',x1,xM,'r',x1,arrx,'k:');

line([5.3 5.3],[4.417 12]);line([5.25 12],[4.417 4.417]);
line([3 3],[2.5 12]);line([3 12],[2.5 2.5]);
line([8 8],[6.667 12]);line([8 12],[6.667 6.667]);
%
text(2.8,10,'\rmu^0=12','fontsize',10,'FontName','Tahoma','Color','b','HorizontalAlignment','right');
text(5.8,6,'\rmu^1=22','fontsize',10,'FontName','Tahoma','Color','b');
text(8.3,10,'\rmu^2=32','fontsize',10,'FontName','Tahoma','Color','b');

xpos=9.5;
ypos=1.0;
text(xpos,ypos,'\rm4x_1+6x_2 =48','color','r','fontsize',10,'FontName','Tahoma','HorizontalAlignment','right');
set(gca,'XTickMode','Manual');set(gca,'XTick',[0 2 4 6 8 10 12]);
set(gca,'YTickMode','Manual');set(gca,'YTick',[0 2 4 6 8 10 12]);

 axis([0 12 0 12]);
 grid off;
 box on;hold off;
 title('Leontiefova užitková funkce');
%
% *** Obr(2E) *********** užitková funkce přímý ADDILOG *****************
%
 uti1=9;
 uti2=12;
 uti3=15;
 alf1=2;
 alf2=3;
 bet1=0.4;
 bet2=0.6;
%
 add1=((uti1-alf1.*x1.^bet1)./alf2).^(1./bet2);
 add2=((uti2-alf1.*x1.^bet1)./alf2).^(1./bet2);
 add3=((uti3-alf1.*x1.^bet1)./alf2).^(1./bet2);
%
 subplot(2,3,5);
 gr=plot(x1,add1,'b',x1,add2,'b',x1,add3,'b',x1,xM,'r');

xpos=7.45;
ypos=3.1;
text(xpos,ypos,'\rm4x_1+6x_2 =48','color','r','fontsize',10,'FontName','Tahoma');
set(gca,'XTickMode','Manual');set(gca,'XTick',[0 2 4 6 8 10 12]);
set(gca,'YTickMode','Manual');set(gca,'YTick',[0 2 4 6 8 10 12]);

text(2,4,'\rmu^0=9','fontsize',10,'FontName','Tahoma','Color','b');
text(4.5,6,'\rmu^1=12','fontsize',10,'FontName','Tahoma','Color','b');
text(7.8,8.6,'\rmu^2=15','fontsize',10,'FontName','Tahoma','Color','b');

 axis([0 12 0 12]);
 title('užitková funkce přímý ADDILOG');
%
% *** Obr(2F) ************* Stone-Gearyho užitková funkce ***************
%
hold on;
 z1=0.1:0.01:12;
 xM=(M0-p1.*z1)./p2;
 pro1=3;
 pro2=6;
 pro3=9;
 gama=2;
 alf1=2/5;
 alf2=3/5;
 del1=0.1;
 del2=3;
 sgf1=(pro1./(gama.*(z1-del1).^alf1)).^(1/alf2)+del2;
 sgf2=(pro2./(gama.*(z1-del1).^alf1)).^(1/alf2)+del2;
 sgf3=(pro3./(gama.*(z1-del1).^alf1)).^(1/alf2)+del2;
 subplot(2,3,6);
 gr=plot(z1,sgf1,'b',z1,sgf2,'b',z1,sgf3,'b',z1,xM,'r');

 axis([0 12 0 12]);
 box on;
 title('Stone-Gearyho užitková funkce');hold off;

xpos=9.5;
ypos=1.0;
text(xpos,ypos,'\rm4x_1+6x_2 =48','color','r','fontsize',10,'FontName','Tahoma','HorizontalAlignment','right');
set(gca,'XTickMode','Manual');set(gca,'XTick',[0 2 4 6 8 10 12]);
set(gca,'YTickMode','Manual');set(gca,'YTick',[0 2 4 6 8 10 12]);

text(2.0,4.7,'\rmu^0=3','fontsize',10,'FontName','Tahoma','Color','b');
text(4.5,5.8,'\rmu^1=6','fontsize',10,'FontName','Tahoma','Color','b');
text(7.1,6.8,'\rmu^2=9','fontsize',10,'FontName','Tahoma','Color','b');
%-------------------------------------------------------------------------
set(gcf,'Name','Příklady indiferenčních křivek dvoukomoditních užitkových funkcí');
 pause