
clear all 
close all
clc

%síla od prvni packy
%zname veliciny
%rozmery paky mm
a=30;
b=76.2;
c=52;
d=47;
alpha=0.7853981634;%zadat v rad
ksi=147*pi()/180;
S1=464.15; %zadana sila soukup_opraveno
%dalsi uhly vyjadrime
gama=asin(d*sin(ksi)/b);
beta=gama+alpha;

%vysledna sila N
S=(S1*b)/(c*sin(beta))


%Sila Fp v N
Fp=56;
%rozmìry mm
AB=26;
AC=45;
AD=22;
AEx=0;
AEy=44.152;
alfa=10*pi()/180;%rad


% %Konkretni hodnota pro packu v poloze 147°
% %ksi=147*pi()/180; viz vyse
% 
%  Fip=1-ksi;%zmena dle vykresu
%  Fif=0;%dle schématu na papíru
%  Fis=1+ksi;%zmena dle vykresu
%     
%  F_147stupnu=((AB*sin(ksi-alfa)*S*cos(Fis))-(AB*cos(ksi-alfa)*S*sin(Fis))-...
%         (AD*cos(ksi)*Fp*sin(Fip))-(AD*sin(ksi)*Fp*cos(Fip)))/((AC*cos(ksi)*sin(Fif))+...
%         (AC*sin(ksi)*cos(Fif)))
    
 %rozsah Fiz v [rad] !!!
Fiz=((51:51+45).*pi())./180;   
%funkce sily F
a=size(Fiz);
for i=1:a(2)
    
    Fip(i)=180-atan((AD*cos(Fiz(i))-AEx)/(AEy+AD*sin(Fiz(i))));%
    Fif(i)=Fiz(i)-(90*pi()/180);%dle schematu na papíru
    Fis(i)=Fiz(i)-alfa;%

    F(i)=((AB*sin(Fiz(i)-alfa)*S*cos(Fis))-(AB*cos(Fiz(i)-alfa)*S*sin(Fis))-...
        (AD*cos(Fiz(i))*Fp*sin(Fip))-(AD*sin(Fiz(i))*Fp*cos(Fip)))/((AC*cos(Fiz(i))*sin(Fif))+(AC*sin(Fiz(i))*cos(Fif)));
    
    
end

Fiz_plot=Fiz*180/pi();
figure(1)
plot(Fiz_plot,F)

