integrate (u*b*(1-((y-h)/h)^2)) dy from y=0 to y=h