求一个三维变分的程序