%% 1.)Definition 
%Parameters
G=[0, 10 ,20 ,30 ,40, 50, 60, 70, 80,  100]; %distance from the front  
S=[0, 7, 16, 25, 34, 44, 54, 64, 74, 86, 100]; 
h_g=[35, 77, 73, 54, 50, 47, 45, 44, 43, 0]; % pressure read from manometer
h_s=[35, 61, 64, 64, 73, 72, 71, 70, 69, 67, 0] ;

b=0.113; %[m]
L=0.1; %[m]
%atmosphere, manometer
h_01=4;
h_02=19;
h_03=27;
dh_g=zeros(1,10);
dh_s=zeros(1,11) ;

%% 2.) Calculations
for i=1:3
dh_g(i)=(h_02-h_g(i));
end
for k= 4:9
    dh_g(k)=(h_01-h_g(k));
end
for l=1:4
    dh_s(l)=(h_02-h_s(l));
end
for m=5:10
    dh_s(m)=(h_03-h_s(m));
end

y_g=10*dh_g;
y_s=10*dh_s;
%% 3.) Plot

plot(G,abs(y_g));
hold on ;
plot (S,y_s);
grid on;
xlabel('x[mm]');
ylabel('p[Pa]');
title('Rozkład ciśnienia wokół profilu płata')

%% 4.) Calculation of bearing force
%trapezoidal rule
P_g=trapz(G,abs(y_g));
P_s=trapz(S,abs(y_s));
%difference of the fields
D_mm=abs(P_g-P_s);
D_cm=D_mm/100;
%1cm^2 - 10 Pa
D_p=D_cm*10;
%Bearing force
R_y=b*L*D_p; %[N]