#include "mex.h" 
void mexFunction(int nlhs, mxArray *plhs[],int nrhs, const mxArray *prhs[]) 
{ 
 double *inData; 
 double *outData; 
 int M,N; 
 int i,j;
 //異常處理 
 if(nrhs!=1) 
   mexErrMsgTxt("USAGE: b=reverse(a)\n"); 
 if(!mxIsDouble(prhs[0])) 
   mexErrMsgTxt("the Input Matrix must be double!\n"); 
 inData=mxGetPr(prhs[0]); 
 M=mxGetM(prhs[0]); 
 N=mxGetN(prhs[0]); 
 plhs[0]=mxCreateDoubleMatrix(M,N,mxREAL); 
 outData=mxGetPr(plhs[0]); 
 for(i=0;i<M;i++) 
   for(j=0;j<N;j++) 
    outData[j*M+i]=inData[(N-1-j)*M+i]; 
}