h = -4/(Lx*Lx)*(x-Lx)*x*np.sin(math.pi*y/Lx)