clc //initialisation of variables g=264//quartz p=0.39//percent //CALCULATIONS S=(1-p)*(g-1)//in //RESULTS printf('the hydralic gradient and seepage velocity=% f in',S)