function r=force(x,y)
    r=8*sin(6*x).*cos(5*y);