能帮忙找找这段程序出错的原因吗？分数可加

#include <iostream.h>
//定义矩阵类
class matrix
{
short rows, cols;//表示矩阵的行、列
double * elems;//double数组存放矩阵中的各元素public:
matrix(short rows,short cols);
~matrix();

//重载运算符“（）”，用来返回元素值
double operator() (short row,short col);
//给元素赋值
void setelem(short row,short col,double val);
//重载运算符“+”，实现距阵相加
friend matrix operator + (matrix p,matrix q);

//输出矩阵中的个元素
void print();
};//定义构造函数
matrix::matrix(short rows, short cols)
{
matrix::rows = rows;
matrix::cols = cols;
elems = new double[rows * cols];
}//定义析构函数
inline matrix::~matrix()//inline表示函数是内涵的
{
delete elems;
}double matrix::operator()(short row,short col)
{
return ((row>=1 && row<=rows && col>=1 && col<=cols) ? elems[(row-1)*cols+(col-1)]:0.0);
}//为矩阵赋值
void matrix::setelem(short row,short col,double val)
{
if(row>=1 && row<=rows && col>=1 && col<=cols)
elems[(row-1)*cols + (col-1)]=val;//数组顺序存放，计算出元素下标
}//“+”
matrix operator + (matrix p, matrix q)
{
matrix m(p.rows,q.cols);
if(p.rows != q.rows || p.cols !=q.cols)
return m;
for(int r=1; r<=p.rows; r++)
for(int c=1; c<=p.cols; c++)
m.setelem(r,c,(p(r,c) + q(r,c)));
return m;
}//输出元素各值
void matrix::print()
{
for(int r=1; r<=this->rows; ++r)
{
for(int c=1; c<=this->cols; c++)
cout<<(*this)(r,c)<<" ";
cout<<"\n";//调用重载运算符“（）”
}
}//主函数
void main()
{
matrix a(2,3),b(2,3),d(2,3);//创建矩阵对象
a.setelem(1,1,1.0);
a.setelem(1,2,2.0);
a.setelem(1,3,3.0);
a.setelem(2,1,4.0);
a.setelem(2,2,5.0);
a.setelem(2,3,6.0);
b.setelem(1,1,1.0);
b.setelem(1,2,2.0);
b.setelem(1,3,3.0);
b.setelem(2,1,4.0);
b.setelem(2,2,5.0);
b.setelem(2,3,6.0);
d=a+b;
d.print();}