function f=pressure_gradient(v,y,t,V1,V2,mu,rho,h)
% function PRESSURE_GRADIENT for FLUID_FLOW program

f=-(12*mu/h^2)*((V1+V2)/2)/rho;