def METHOD_NAME():
    solution = sin(3*pi/2 * x)**2 * sin(3*pi/2 * y)**2
    f        = laplace(laplace(solution))
    dir_zero_boundary = get_boundaries(1, 3)
    l2_error, h1_error, h2_error = run_biharmonic_2d_dir(solution, f,
            dir_zero_boundary, ncells=[2**3, 2**3], degree=[3, 3])
    print()
    print(l2_error)
    print(h1_error)
    print(h2_error)
    print()
    assert False