Выделение памяти для массива

void matrix::read(ifstream &f)
{
    int i=0;
    int j=0;
    while(f.peek()!=EOF)
    {
        a=new double *;
        j=0;
        while (f.peek()!='\n')
        {
            a=new double;
            f>>b[j];
            j++;
        }
        a[i]=b;
        i++;
        f.ignore();
    }
    n=i;
    m=j;
}

#include <fstream>
using namespace std;
class matrix
{
private:
    int n;
    int m;
    double **a;
public:
    //конструкторы
    matrix(int a, int b); //создает нулевую матрицу размером nxm
    matrix();//создает нулевую матрицу
    matrix (int k,int l, bool E);//создает матрицу размером nxm , нулевую если E=0 , единичную если E=1
    void scanmatrix(); //ввод матрицы с клавиатуры
    void scanmatrix(int, int);
    void randmatrix();
    void printmatrix();
    void record(ofstream &);
    void read(ifstream &);
 
    double getitem(int ,int);
    void putitem(double, int, int);
 
    double mean();
 
    matrix operator + (matrix);
    matrix operator - (matrix);
    matrix operator * (matrix);
    matrix operator * (double k);
    void transp();
};