共轭梯度法C语言(西安交大)(共3页).docx

精选优质文档-倾情为你奉上#include#include#define N 10 /*定义矩阵阶数*/void main()int i,j,m,ANN,BN;double XN,akvN,dkaN,rkN,dkN,pk,pkk,ak,bk;for(i=0;iN;i+) /*输入系数矩阵A*/for(j=0;jN;j+)if(i=j)Aij=-2;else if(abs(i-j)=1) Aij=1;else Aij=0; for(i=0;iN;i+) /*输入矩阵B*/ if(i=0|i=N-1)