% This M file uses function thickness to estimate
% various pavement thicknesses for CBR values
% ranging from and initial value (ini_cbr_value) to a
% final CBR value (fin_cbr_value)

P = 60000;		% load (pounds)
A = 254;		% tire contact area

in_cbr_value 	= 10;
fin_cbr_value	= 30;
npoints			= 100;

cbr_vector = linspace(in_cbr_value,fin_cbr_value,npoints);

t_vector = thickness(P,A,cbr_vector);

plot(cbr_vector,t_vector)
xlabel('CBR')
ylabel('Pavement Thickness (in)')
grid
