Ошибка: нарушение прав доступа при чтении

#include <iostream>
#include <math.h>
#include <vector>
 
using std::cout;
using std::cin;
using std::endl;
using std::vector;
 
class macierz
{
    void create_table();
    int rows;
    int columns;
public:
    int getc(){return columns;}
    int getr(){return rows;}
    macierz(macierz,macierz*,int);
    macierz(int,int);
    macierz();
    ~macierz();
    double **cells;
    void clear_macierz();
    double det(vector<bool>,int,int);
    double det();
    //double det(macierz*,macierz*,int);
};
 
void macierz::clear_macierz()  //очистка матрицы
{
    for(int i=0;i<rows;i++)
        for(int j=0;j<columns;j++)
            cells[i][j]=0;
}
 
macierz::macierz(int a,int b)   //создание матрицы если дано а и б
{
    rows=a;
    columns=b;
    create_table();
    clear_macierz();
}
 
macierz::macierz()            //создание матрицы вручную
{
    cout<<"ilosc wierszy macierzy = ";
    cin>>rows;
    cout<<"ilosc kolumn macierzy = ";
    cin>>columns;
    create_table();
    for(int i=0;i<rows;i++){
        for(int j=0;j<columns;j++){
            cout<<"a("<<i+1<<","<<j+1<<") = ";
            cin>>cells[i][j];
    }}
}
 
macierz::macierz(macierz A,macierz* B,int u)     //создание матрицы из заменой столбцов (метод Крамера)
{
    columns=rows=A.getr();
    create_table();
    for(int i=0;i<rows;i++)
    {
        for(int j=0;j<rows;j++)
        {
            if(j!=u)
            {this->cells[i][j]=A.cells[i][j];}
            else{this->cells[i][j]=B->cells[i][0];}
        }
    }
}
 
macierz::~macierz()
{
    for(int i=0;i<rows;i++)
        delete []cells[i];
    delete []cells;
}
 
void macierz::create_table()     //выделение памяти
{
    cells=new double*[rows];
    for(int i=0;i<rows;i++)
        cells[i]=new double[columns];
    }
 
double macierz::det()    //возвращает детерминант
{
    vector<bool> a(columns);
    int ii=0,aj=0;
    return det(a,ii,aj);
}
 
double macierz::det(vector<bool> a,int ii,int aj)
{
    double wyz=0;
    static bool aa=0;
    if(aa){a[aj]=1;}else{aa=1;}
    int oo=0;
    for(int j=0;j<columns;j++)
    {cout<<"a[]= "<<a[0]<<a[1]<<endl;
        if(!a[j])
        {
            if(ii==(columns-1)){return cells[ii][j];}
            else{wyz+=cells[ii][j]*pow((double)-1,oo)*det(a,ii+1,j);oo++;}  
        }
    }
    return wyz;
}
 
macierz* mnoz(macierz A,macierz B)   //множение двух матриц
{
    macierz *C=new macierz(A.getr(),B.getc());
    for(int i=0;i<A.getr();i++)
        for(int j=0;j<B.getc();j++)
            for(int u=0;u<A.getc();u++)
                C->cells[i][j]=C->cells[i][j]+A.cells[i][u]*B.cells[u][j];
    return C;
}
 
 
void kramer(macierz A,macierz *B)     //нахождение неизвесных с помощью метода крамера
{
    //int j=0;
    double wyz=A.det();
    for(int j=0;j<A.getc();j++)
    {
        //X=new macierz(*A,B,j);
 
        double wy=macierz(A,B,j).det();
        cout<<"x["<<j+1<<"] = "<<wy/wyz<<endl;
        cout<<"wy["<<j+1<<"] = "<<wy<<endl;
        cout<<"wyz["<<j+1<<"] = "<<wyz<<endl;
        //delete X;
    }
}
 
 
double det(macierz *A,macierz*B,int u)     //возвращает детерминант матрицы А, в которой заменен столбец u матрицей В
{
    vector<double> vec(A->getr());
    for(int i=0;i<A->getr();i++)
    {
        vec[i]=A->cells[i][u];
        //cout<<A->getr()<<"  "<<i<<endl<<A->cells[i][u];
        A->cells[i][u]=B->cells[i][0];}//cout<<A->cells[i][u]<<endl;}
    cout<<A->cells[0][0]<<" "<<A->cells[0][1]<<" "<<A->cells[0][2]<<endl
        <<A->cells[1][0]<<" "<<A->cells[1][1]<<" "<<A->cells[1][2]<<endl
        <<A->cells[2][0]<<" "<<A->cells[2][1]<<" "<<A->cells[2][2]<<endl;
    double wyzn=A->det();
    cout<<"wyzn=A->det() = "<<wyzn<<endl;
    for(int i=0;i<A->getr();i++)A->cells[i][u]=vec[i];
    return wyzn;
}
 
void main()
{
    //setlocale(0,"En");
    macierz *A,*B;//,*X;
    cout<<"Macierz A:"<<endl;
    A=new macierz;
    cout<<"Macierz B:"<<endl;
    B=new macierz;
    //cout<<det(A,B,0);
 
    /////double wyz=A->det();
    /////cout<<"wyz=A->det() = "<<wyz<<endl;
    /////for(int j=0;j<(A->getc());j++)
    /////   cout<<"x["<<j+1<<"] = "<<det(A,B,j)/wyz<<endl;
    
    //if(A->det()!=0)
    kramer(*A,B);
    //macierz *C=mnoz(*A,*B);
    //macierz *A;
    //cout<<"Macierz A:"<<endl;
    //A=new macierz;
    
    //double wyz=A->det();
    //cout<<wyz<<endl;
    //cout<<"x[1] = "<<(macierz(A,B,0).det())<<endl;
    
    //macierz *X=new macierz(A,B,1);
    
    /*
    X=new macierz(A,B,2);
    cout<<X->cells[0][2]<<endl;
    cout<<X->cells[1][2]<<endl;
    cout<<X->cells[2][2]<<endl;
    cout<<"detX = "<<macierz(A,B,2).det();
    */
    //system("PAUSE");
getchar();getchar();getchar();
}