求解PDE变量系数

Question

我正在尝试使用PDEtoolbox解决某个PDE（2个变量）。

等式为-(laplacian)u + V*u = 1，椭圆等式c和f等于1，a=V是一个2变量函数，定义于m文件。 我对V的定义存在一些问题。这是文件的代码（M是已经定义的全局矩阵）

function v = V(x,y)

global M

[i,j]=size(M)

intx=int32(floor(i*x)+1)

inty=int32(floor(j*y)+1)

v = M(intx,inty)

可以绘制函数，但是当我尝试将V(x,y)放在＆＃39; PDE规范＆＃39;中的a字段时，当我要求解决方程时，我会得到以下错误消息：

'Expression evaluates to wrong size. Must be scalar or row vector. In a system case, pass first or second row; for example u(2,:). In expression: '-ux' when evaluating pde coefficients.'

欢迎任何帮助，解释为什么这不起作用或给我另一种解决方程的方法。