D2=8*0.0254:2*0.0254:20*0.0254;
Fr2=6804.5*(D2).^2+561.0094*(D2).^3;
plot(D2,Fr2)

xlabel('angle of inclination')
ylabel('resultant force')
title('case#2 fr')