def METHOD_NAME():
    i: i32; j: i32; k: i32; l: i32; size: i32;
    array1d: f32[10] = empty(10, dtype=float32)
    sinh1d: f32[10] = empty(10, dtype=float32)
    for i in range(10):
        array1d[i] = f32(f64(i)/10.0)
    sinh1d = sinh(sinh(array1d))
    verify1d(array1d, sinh1d, 10)
    arraynd: f64[40, 10, 16, 2] = empty((40, 10, 16, 2), dtype=float64)
    sinhnd: f64[40, 10, 16, 2] = empty((40, 10, 16, 2), dtype=float64)
    size = 40 * 10 * 16 * 2
    for i in range(40):
        for j in range(10):
            for k in range(16):
                for l in range(2):
                    arraynd[i, j, k, l] = float(i + 2*j + 3*k + 4*k)/float(size)
    sinhnd = (sinh(arraynd) + 2.0)/2.0
    verifynd(arraynd, sinhnd, 40, 10, 16, 2)