slvsm2.c

来自「适合大型数值计算代码 现在网络上已经找不到了 购买需要20$」· C语言 代码 · 共 14 行

#include <math.h>

void slvsm2(double **u, double **rhs)
{
	void fill0(double **u, int n);
	double disc,fact,h=0.5;

	fill0(u,3);
	fact=2.0/(h*h);
	disc=sqrt(fact*fact+rhs[2][2]);
	u[2][2] = -rhs[2][2]/(fact+disc);
}