function z = uexact3(x,y)

    z = (x*x+y*y)*sin(x);