julia
ut=−uy−x−2+ϵ(uxx+uyy) for ϵ=1 in [−1,1]2, subject to a homogeneous Dirichlet boundary condition. Use m=n=30