求助：二阶PDE的NR数值解法# Computation - 科学计算
L*g
1 楼
边界值问题的2阶PDE，Numerical Recipe给的程序是mglin()，具体操作本人一时难以掌
握，希望有经验的人指点一二。具体如下：
方程u''(x)+u''(y)=rho, 边界条件f(x,y)=0
1. mglin()里说一开始输入的u(n,n)里应该包含rho，具体应该怎么包含？
2. 边界条件f(x,y)=0在程序的什么地方输入？
3. resid()中的res[i][1]=res[i][n]=...=0是否要根据我的边界条件修改？
握，希望有经验的人指点一二。具体如下：
方程u''(x)+u''(y)=rho, 边界条件f(x,y)=0
1. mglin()里说一开始输入的u(n,n)里应该包含rho，具体应该怎么包含？
2. 边界条件f(x,y)=0在程序的什么地方输入？
3. resid()中的res[i][1]=res[i][n]=...=0是否要根据我的边界条件修改？