"""
        Para X - A
    """
    for i in range(Ny):
        for j in range(Nx):
            a1x[i,j] = max( abs((uv[i,j]/fi)*df(u_plus1x[i,j],ex)) , abs((uv[i,j]/fi)*df(u_minus1x[i,j],ex)) )
            a2x[i,j] = max( abs((uv[i,j]/fi)*df(u_plus2x[i,j],ex)) , abs((uv[i,j]/fi)*df(u_minus2x[i,j],ex)) )

    """
        Para Y - A
    """
    for j in range(Nx):
        for i in range(Ny):
            a1y[i,j] = max( abs((vv[i,j]/fi)*df(u_plus1y[i,j],ex)) , abs((vv[i,j]/fi)*df(u_minus1y[i,j],ex)) )
            a2y[i,j] = max( abs((vv[i,j]/fi)*df(u_plus2y[i,j],ex)) , abs((vv[i,j]/fi)*df(u_minus2y[i,j],ex)) )


    """
        Para X - H
    """
    for i in range(Ny):
        for j in range(Nx):
            H1x[i,j] = ( f(u_plus1x[i,j],ex) + f(u_minus1x[i,j],ex))/(2.0*fi) - a1x[i,j]*(u_plus1x[i,j]-u_minus1x[i,j])/2.0

    for i in range(Ny):
        for j in range(Nx):
            H2x[i,j] = (f(u_plus2x[i,j],ex)+f(u_minus2x[i,j],ex))/(2.0*fi) - a2x[i,j]*(u_plus2x[i,j] - u_minus2x[i,j])/2.0

    duX = (H2x - H1x)/h


    """
        Para Y - H
    """
    for j in range(Nx):
        for i in range(Ny):
            H1y[i,j] = ( f(u_plus1y[i,j],ex) + f(u_minus1y[i,j],ex))/(2.0*fi) - a1y[i,j]*(u_plus1y[i,j]-u_minus1y[i,j])/2.0

    for j in range(Nx):
        for i in range(Ny):
            H2y[i,j] = (f(u_plus2y[i,j],ex)+f(u_minus2y[i,j],ex))/(2.0*fi) - a2y[i,j]*(u_plus2y[i,j] - u_minus2y[i,j])/2.0

    duY = (H2y - H1y)/h