d2=8:2:20
D2=d2*0.0254

a=3.8608+15.4432*D2+4.25*(D2).^2
b=61.7728+16*D2;
yR2=a./b
plot(D2,yR2)
xlabel('angle of inclination')
ylabel('resultant force')
title('case#2')