První částí práce je odvození rovnice difuse, přičemž byla za tímto účelem odvozena Gaussova věta. Další část je věnována velmi krátce Fourierovým řadám, aby bylo možno pochopit a popsat Fourierovu metodu pro řešení této parciální diferenciální rovnice s Neumannovými a Dirichletovými okrajovými podmínkami. Dále je v této práci popsána metoda sítí (Explicitní metoda a Implicitní metoda v 1D i 2D a Crankovo - Nicolsonovo schéma v 1D), vše naprogramováno v Matlabu. Některé výsledky jsou srovnány s analytickým řešením.The first part of the thesis adresses the derivation of the diffusion equation and for this purpose the Gauss Theorem is derived. The next part adresses shortly the Fourier Series so as to understand and describe The Fourier method for solution of this partial differential equation with Neumann and Dirichlet boundary conditions. Then the Finite difference method is described (explicit method, implicit method in 1D and 2D and Crank-Nicolson method in 1D), all programmed in Matlab. Some results are compared with analytical solution.