void multiply(double*p1,double*p2,double*pmatrix,int n)
{
	    double t;
		//pmatrix=new double[n*n];
		for(int e=0;e<n;e++)
		{
			for(int f=0;f<n;f++)
			{
                t=0;
				for(int g=0;g<n;g++)
					t+=p1[e*n+g]*p2[g*n+f];
					pmatrix[e*n+f]=t;
			}
		}
		/*for(int h=0;h<n;h++)
		{
		    for(int w=0;w<n;w++)
			{
			    cout.precision(6); 
                cout.setf(ios::scientific,ios::floatfield); 
                cout<<pmatrix[h*n+w]<<" ";
			}
			cout<<endl;
		}
		delete[n*n] pmatrix;*/ 
}