function [Mz,ALFA,BETA,L6,Fi,Zi,Zk,Fi3,BETA24] = MZaALFAproPCHM2(L1x,L1y,L2x,L2y,L3,L4,L5x)

%PREDDEFINOVANE ROZMERY (mm)
L5y=0;%cep PCHM na zakladu (E)
L6y=0;%pocatecni poloha druheho konce paky (F)
%PCHM
Lp=505;%min delka PCHM
Z=340;%zdvih PCHM
d=50;%prumer pistnice
S=(pi*d^2)/4;%plocha prurezu pistnice (mm2)
p=12; %tlak hydrauliky (MPa)
N0=p*S;%sila v pistnici (N)

%dopoctene rozmery
L1=sqrt(L1x^2+L1y^2); %vzdalenost cepu na zakladu
L2=sqrt(L2x^2+L2y^2); %vzdalenost cepu na lopate
alfa1=acos(L1x/L1); %uhel cepu na zakladu vuci vodorovine

%konfigurace v pocatecni poloze
alfa2=acos(L2x/L2);
l5=sqrt((L1x-L2x)^2+(L1y-L2y)^2);
alfa15=acos((L1^2+l5^2-L2^2)/(2*L1*l5));
alfa35=acos((L3^2+l5^2-L4^2)/(2*L3*l5));
alfa3=alfa35+alfa15-(pi-alfa1);
fi3=pi/2-alfa3;%rad pocatecni odkloneni delsiho ramene paky od svislice
Fi3=fi3*180/pi;%deg

%dopoctene rozmery druheho ramena paky
L6x=L5x-Lp;%druhy rozmer ramene aby byl PCHM v rovine
L6=round(sqrt((L6x-L1x)^2+(L6y-L1y)^2)-0.5);%delka druheho ramene paky (zaokrouhlena)
Fi=round(((fi3+atan((L6x-L1x)/(L1y-L6y)))*180/pi)-0.5);%deg uhel mezi rameny paky
fi=Fi*pi/180;%rad
L6x=L1x+L6*sin(fi-fi3);
L6y=L1y-L6*cos(fi-fi3);

%VYPOCET UHLU ALFA A BETA

%vypocet uhlu pohybu paky (beta)
LP=sqrt((L6x-L5x)^2+(L6y-L5y)^2);%delka neuplne zasunute pistnice
L15=sqrt((L1x-L5x)^2+(L1y-L5y)^2);%vzdalenost mezi body B a E
psiZ=acos((L15^2+L6^2-(LP)^2)/(2*L15*L6));
psiK=acos((L15^2+L6^2-(Lp+Z)^2)/(2*L15*L6));
beta=psiK-psiZ; %rad
BETA=beta*180/pi; %deg

%vypocet uhlu pohybu lopaty (alfa)
beta3=alfa3-beta;
beta13=beta3+pi-alfa1;
l6=sqrt(-cos(beta13)*(2*L1*L3)+L1^2+L3^2);
beta26=acos((L2^2+l6^2-L4^2)/(2*L2*l6));
beta16=acos((L1^2+l6^2-L3^2)/(2*L1*l6));
beta2=alfa1-beta16-beta26;
beta24=acos((L2^2+L4^2-l6^2)/(2*L2*L4));
BETA24=beta24*180/pi; %deg
alfa=alfa2-beta2; %rad
ALFA=alfa*180/pi; %deg

%VYPOCET POHYBU A CHARAKTERISTIK

%nepohyblive body (zaklad/ram)
xA=0;
yA=0;
xB=L1x;
yB=L1y;
xE=L5x;
yE=L5y;

%pocatecni podminky cyklu
n=340; %pocet poloh v grafu
delta3=alfa3; %pocatecni odklon delsiho ramene paky od horizontaly
zi=LP-Lp; %pocatecni zdvih pistnice
ni=(Lp+Z-LP)/(n-1); %krok zdvihu (rovnomerny)

%cyklus
for i=1:n
    %moment na pace
    Mo(i)=(N0*L6*sin(acos(((Lp+zi)^2+L6^2-(L15)^2)/(2*(Lp+zi)*L6))))/1000;
    %souradnice bodu C a F
    xC(i)=xB+L3*cos(delta3);
    yC(i)=yB+L3*sin(delta3);
    xF(i)=xB-L6*cos(delta3+fi);
    yF(i)=yB-L6*sin(delta3+fi);
    %pomocne rozmery a uhly
    l7=sqrt((xC(i))^2+(yC(i))^2);
    delta17(i)=acos((L1^2+l7^2-L3^2)/(2*L1*l7));
    delta27(i)=acos((L2^2+l7^2-L4^2)/(2*L2*l7));
    delta2(i)=alfa1-delta17(i)-delta27(i);
    %souradnice bodu D
    xD(i)=L2*cos(delta2(i));
    yD(i)=L2*sin(delta2(i));
    %vypocet vyklopneho momentu
    delta24(i)=acos(((xC(i)-xD(i))*(-xD(i))+(yC(i)-yD(i))*(-yD(i)))/(L2*L4));
    delta34(i)=acos(((xB-xC(i))*(xD(i)-xC(i))+(yB-yC(i))*(yD(i)-yC(i)))/(L3*L4));
    N4(i)=(Mo(i)/(L3*sin(delta34(i))))*1000;%sila v tahle (vzper)
    Mz(i)=(N4(i)*L2*sin(delta24(i)))/1000;%vyklopny moment
    %dopocet velicin do dalsiho cyklu
    zi=zi+ni;
    delta3=alfa3-(acos((L15^2+L6^2-(Lp+zi)^2)/(2*L15*L6))-psiZ);
end
Zk=zi-ni;
Zi=LP-Lp; %pocatecni zdvih pistnice
end