/*****************C=A+B******************/
void add(double *A,double *B,double *C,int M,int N) 
{ 
	int m,n; 
	
	for(m=0;m<M;m++) 
		for(n=0;n<N;n++) 
		{ 
			C[m*N+n]=A[m*N+n]+B[m*N+n]; /*相加a+b*/ 
		} 
}