использование .h файлов в описании класса С++

[shir994]

Как правильнее описывать класс:полностью в .h файле или в .cpp,а в .h выносить только описания функций?.И вот второй вопрос,почему прога вылетает,говоря,что происходит попытка записи в занятую память.
Ошибка явно где-то в описании+=,но не могу понять где,вроде проверял,значение всех переменных верные.

h.h
#include<iostream>
using namespace std;

class mer
{
public:
	double x,y,z,t;
	mer(double x1=0,double x2=0,double x3=0,double x4=0)
		:x(x1),y(x2),z(x3),t(x4)
	{}
	mer(const mer& q):x(q.x),y(q.y),z(q.z),t(q.t)
	{}
	//~mer()
	//{}
    mer& operator=(const mer& q)
	{   
		x=q.x;y=q.y;z=q.z;t=q.t;
	
		return *this;
	};
	mer operator +=(const mer& q)
	{  x=x+q.x;
	   y=y+q.y;
	   z=z+q.z; 
       t=t+q.t;
	   return(*this);
	}
      mer operator-=(const mer& q)
	{  x=x-q.x;
	   y=y-q.y;
	   z=z-q.z; 
       t=t-q.t;
	   return(*this);
	}
	mer operator+(const mer& q)
	{
		mer p(*this);
		p+=q;
		return p;
	}
		mer operator-(const mer& q)
	{
		mer p(*this);
		p-=q;
		return p;
	}
	mer operator*(const double& q)
	{
		mer p(*this);
		p.x=x*q;
		p.y=y*q;
		p.z=z*q;
		p.t=t*q;
		return p;
	}
	mer operator/(const double& q)
	{
        mer p(*this);
		p.x=x/q;
		p.y=y/q;
		p.z=z/q;
		p.t=t/q;
		return p;
	}
};

ostream& operator<<(ostream& os,const mer& q)
{
	cout<<"(" << q.x << ","<<q.y<<","<<q.z<<","<<q.t<<")"<<flush;
    return os;
}

istream& operator>> (istream& is,mer& q)
{
  double x,y,z,t; is >>x ;
  char c;       is >> c;
  if (c==',') is>> y;
  is >> c;
  if (c==',') is>> z;
  is >> c;
  if (c==',') is>> t;
  q=mer(x,y,z,t);
  return is;
}
main.cpp
#include<iostream>
#include"h.h"
using namespace std;
class vektor
{
protected:
	mer* p;
	int N;
public:
	vektor (int n=0):p(NULL),N(0)
	{
		if(n<=0) return;
		p=new mer[n];
		if(p) N=n;
	}
	~vektor()
	{
		delete p;
		cout<<"rip";
	}
	mer& operator [](int i)
	{
		if (i<0) {cout<<"error";}
		return p[i];
	}
    
	vektor (const vektor& s)
	{
		p=new mer(s.N);
		N=s.N;
		for (int i=0;i<N;i++)
			p[i]=s.p[i];
	}
	int size()const {return N;}
	void insert(int n,mer x)
	{
	mer*	p1=new mer[N+1];
	int i;
	for(i=N;i>n;i--)
		p1[i]=p[i-1];
	p1[i]=x;
	for(i=n-1;i>=0;i--)
		p1[i]=p[i];
	delete p;
	p=p1;
	N=N+1;
	}
	void erase(int n)
	{
		mer*	p1=new mer[N-1];
	int i;
	for(i=N-1;i>n;i--)
		p1[i-1]=p[i];
	for(i=n-1;i>=0;i--)
		p1[i]=p[i];
	delete p;
	p=p1;
	N=N-1;
	}
	vektor operator+=(const vektor q)
	{
		for(int i=0;i<1;i++)
		p[i]+=q.p[i];
		return *this;
	}
};



ostream& operator<<(ostream& os,vektor& q)
{
	for(int i=0;i<q.size();i++)
    cout<<q[i]<<' ';
    return os;
}


void main()
{
	int n;
	cout<<"vvedite razmernost"<<endl;
	cin>>n;
	vektor x(n),y(n);
	mer k;
	for(int i=0;i<n;i++)
	{
		k.x=i;k.y=i;k.z=i;k.t=i;
		x[i]=k;
		y[i]=k;
	}
x+=y;
    cout<<x<<endl<<y;
}