急！！！怎样在函数体内将一个临时申请的对象作为返回值传出？

本人在毕业设计中用Ｃ＋＋实现一个算法遇到了一点小问题，大致问题描述就是：怎样在函数体内将一个临时申请的对象作为返回值传出．我卡在这个问题上停滞不前，求各位高手帮我提出一点解决办法，下面是我将主要问题突出的一段简化的代码，程序运行时报了内存出错，大家帮我看看怎样改正一下，谢谢了！！！
matrix.h头文件如下：
class matrix{
public:
matrix(double**,int,int);
matrix();
~matrix();
matrix operator+(matrix &s);//重载两个数组相加
private:
double** p;//指向二维数组的二维指针
int m,n;//m为二维数组的行数，n为二维数组的列数
};
matrix.cpp文件如下：
#include"matrix.h"
#include<iostream.h>
matrix::matrix(double **A,int k,int l)
{
m=k;
n=l;
p=new double*[m];
for(int i=0;i<m;i++)
p[i]=new double[n];
for(i=0;i<m;i++)
for(int j=0;j<n;j++)
p[i][j]=A[i][j];
}
matrix::~matrix()
{
for (int i = 0; i < m; i++)
delete [n]p[i];
delete [m] p;
m=n=0;
}
matrix matrix::operator+(matrix &s)
{
if(s.m==m&&s.n==n)
{
double **r;
r=new double*[m];
for(int i=0;i<m;i++)
r[i]=new double[n];
for(i=0;i<m;i++)
for(int j=0;j<n;j++)
r[i][j]=p[i][j]+s.p[i][j];
matrix sum(r,m,n);
return sum;
}
else
cout<<"error!两个数组不能相加"<<endl;
return *this;

}
main主函数如下：
#include"matrix.h"
#include<iostream.h>
int main()
{
double **a,**b;
for(int i=0;i<4;i++)
for(int j=0;j<4;j++)
a[i][j]=2;
for(i=0;i<4;i++)
for(int j=0;j<4;j++)
b[i][j]=2;
matrix amatrix(a,4,4),bmatrix(b,4,4);
amatrix+bmatrix;
return 0;
}
由于操作数组比较大，所以用的是二维指针，大家帮忙看下，由于毕业设计时间比较紧，请各位早点答复，再次谢谢大家了！！！