武汉理工大学数值分析实验代码（java）

package shuzhifenxi; import java.util.Scanner; public class n超松弛法 { public static void main(String[] args) { // TODO 自动生成的方法存根 Scanner in=new Scanner (System.in); System.out.println("请输入系数矩阵行和列数:"); int n=in.nextInt(); System.out.println("请输入系数矩阵:"); double a[][]=new double[n+1][n+1]; for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) a[i][j]=in.nextDouble(); } System.out.println("请输入方程组右侧向量:"); double b[]=new double[n+1]; for(int i=1;i<=n;i++) { b[i]=in.nextDouble(); } System.out.println("请输入精度e:"); double e=in.nextDouble(); System.out.println("请输入松弛因子w:"); double w=in.nextDouble(); System.out.println("请输入步数N:"); int N=in.nextInt(); int k=1; double x[]=new double[n+1]; for(int i=1;i<=n;i++) {//初始化x x[i]=0; } double y[]=new double[n+1]; for(int i=1;i<=n;i++) {//初始化y y[i]=0; } double Y[]=new double[n+1]; for(int i=1;i<=n;i++) {//初始化y Y[i]=0; } for(k=1;k<=N;k++) { double max=0; for(int i=1;i<=n;i++) { double s=0; for(int j=1;j<=n;j++) { if(j!=i) s=s+a[i] [j]*Y[j]; } y[i]=(b[i]-s)/a[i][i]; Y[i]=w*y[i]+(1-w)*x[i]; double t= Math.abs(Y[i]-x[i]); if(t>max) max=t; } if(max<e) { for(int r=1;r<=n;r++) { System.out.println("x"+r+"="+Y[r]); } break; } for(int q=1;q<=n;q++) { x[q]=Y[q]; } } if(k>N) System.out.println("迭代失败！"); } }