函数未声明(首先使用此函数）

我正在研究矩阵类，但我无法弄清楚为什么会出现此错误，因为我已经定义了有问题的函数（Dopisi（））。这是我到目前为止所拥有的：

这是我到目前为止所拥有的：

#include<IOStream>
using namespace std;
class Matrica{
private:int Duzina;
    int Sirina;
    int** pok;
public:
    Matrica(){//default constructor
        Duzina=0;
        Sirina=0;
        int i,j;
        pok=new int* [Duzina];
        for(i=0;i<Duzina;i++){
            pok[i]=new int[Sirina];
        }
    }
    Matrica(int duz,int sir){//constructor
        Duzina=duz;
        Sirina=sir;
        int i,j;
        pok=new int* [Duzina];
        for(i=0;i<Duzina;i++){
            pok[i]=new int[Sirina];
            for(j=0;j<Sirina;j++){
                pok[i][j]=0;
            }
        }
    }
    void Dodaj(int i,int j){//add one element
        cout<<"Unesite element ["<<i<<"]["<<j<<"]"<<endl;
        cin>>pok[i][j];
    }
    void Popuni(){//fill entire matrix
        int i,j;
        for(i=0;i<Duzina;i++)
            for(j=0;j<Sirina;j++)
                Dodaj(i,j);
    }
    Matrica(Matrica& mat){//copy constructor
        Duzina=mat.Duzina;
        Sirina=mat.Sirina;
        int i,j;
        pok=new int* [Duzina];
        for(i=0;i<Duzina;i++){
            pok[i]=new int[Sirina];
            for(j=0;j<Sirina;j++){
                pok[i][j]=mat.pok[i][j];
            }
        }
    }
    ~Matrica(){}
    int max(int a,int b){//max value
        if(a>b)return a;
        else
            return b;
    }
    void Ispis(){//cout
        int i,j;
        for(i=0;i<Duzina;i++){
            for(j=0;j<Sirina;j++){
                cout<<pok[i][j]<<" ";
            }
            cout<<endl;
        }
    }
    Matrica& operator =(const Matrica& x){
        for(int i=0;i<Duzina;i++)
            delete pok[i];
        Duzina=x.Duzina;
        Sirina=x.Sirina;
        for(int i=0;i<Duzina;i++){
            for(int j=0;j<Sirina;j++){
                pok[i][j]=x.pok[i][j];
            }
        }
        return *this; 
    }
    Matrica Dopisi(const Matrica a,const Matrica b);
};
Matrica Matrica::Dopisi(const Matrica a,const Matrica b){//function must take two      matrices 
    int i,j,duz,sir;
    duz=max(a.Duzina,b.Duzina);
    sir=a.Sirina+b.Sirina;
    Matrica temp(duz,sir);
    for(i=0;i<a.Duzina;i++){
        for(j=0;j<a.Sirina;j++){
            temp.pok[i][j]=a.pok[i][j];
        }
    }
    for(i=0;i<b.Duzina;i++){
        for(j=0;j<b.Sirina;j++){
            temp.pok[i][j+a.Sirina]=b.pok[i][j];
        }
    }
    temp.Ispis();
    // return temp;
}
int main(){
    Matrica mat1(3,3);
    mat1.Popuni();
    mat1.Ispis();
    Matrica mat2(4,4);
    mat2.Popuni();
    mat2.Ispis();
    mat2=mat1;
    mat2.Ispis();
    Matrica mat3;
    mat3=Dopisi(mat1,mat2);
    system("pause");
    return 0;
}

P.S.变量和方法的名称是塞尔维亚语，希望你能弄清楚。提前感谢