//example 4.5
clc; funcprot(0);
c=50;
//from table
Ncq=6.3;
FS=4;
qu=c*Ncq;
qall=qu/4;
disp(qall,"allowed shear stress in kN/m^2");