import numpy as np

def meshFiber_rectangle(h,b,nFibersInRow,nFibersInCol):
    """
    Create a mesh of fibers in a rectangle.
    """

    nTot = nFibersInRow * nFibersInCol  
    coordinates = np.zeros((nTot, 2))

    for i in range(nFibersInRow):
        y = h / 2 - h / (2 * nFibersInRow) - h / nFibersInRow * i
        for j in range(nFibersInCol):
            z = -b / 2 + b / (2 * nFibersInCol) + b / nFibersInCol * j
            coordinates[nFibersInCol * i + j, 0] = y
            coordinates[nFibersInCol * i + j, 1] = z

    Afib = (h * b) / nTot  
    Anf = Afib * np.ones(nTot) 

    return coordinates, Anf