很简单的题目，就是得不出结果，大家帮帮忙，纠下错，或者直接写一个出来，题目在图片上
NDSolve[{D[u[x, y, t], x] + D[v[x, y, t], y] == 0,
20 D[u[x, y, t], t] == -D[P[x, y, t], x] +
0.2[D[u[x, y, t], x, x] + D[v[x, y, t], y, y]],
20 D[v[x, y, t], t] == -D[P[x, y, t], y] +
0.2[D[u[x, y, t], x, x] + D[v[x, y, t], y, y]],
u[0, y, t] == 1 - (y - 1)^2, u[x, 0, t] == 0, u[x, 2, t] == 0,
v[0, y, t] == 2 - y^2, v[x, 0, t] == 0, v[x, 2, t] == 0,
P[10, y, t] == 0, u[x, y, 0] == 0, v[x, y, 0] == 0,
P[x, y, 0] == 0}, {u, v}, {x, 0, 10}, {y, 0, 2}, {t, 0, 60}]

我试了试也不行。。。
估计是书上打印错误了？？
那个第一式等与0，为啥不先带入第二式？
不是我试过了，貌似还是不行。。。
要么你拿Maple试试。。。MMA有点写地方还是超级弱鸡的。。。

边界条件v(0,y,t)=2-y^2和初始条件v(x,y,0)=0不相容，随便取个x=0,y=1,t=0就能看出来。

对外链的限制好像有所放松，试一下：https://mathematica.stackexchange.com/q/212441/1871