% Calculation of minimum deflection
dely = ((HY ((xd 2) + (HL A2)) A.5) cos (( atan (HL/xd)))) / (cos ((pi/180) 45) -
(cos (( atan ( HL/xd )))));
xd (1)= xd;


fori =1:N
forj =2:((HL-xd)/.01)
xd(j)= xd (j- 1)+.01;
% Calculation of minimum stiffness
k () =(( 3 (mp 9.81) (i)) + (6 ((mts + mbs) 9.81) (i 1)) + 6 mts 9.81 .5 sin ((2 *
atan (HL /xd (j) ))))/ (( 6 (((HY ((xd (j) A2) + (HL 2)) A.5) + dely) cos (( atan
(HL/xd(j)))))));
end
kk (i)= max (k());
A (i) = (.4 (HY HL)) kk(i) .7071;
end
fori =1: N
if (A(i) > fxd)
nk (i)=k (i);
ndely (i) = (.4 (HY HL));
elseif(A (i) <fxd)
nk (i)= .166 fxd/ ((.4 (HY-HL)) .7071);
ndely (i) = (.4 (HY HL));
YD (i)= ndely (i) + (HY HL);
LS (i) = HY + ndely (i);
end
end


% Results
disp ('The number of layers on each side ');
disp (N)