function Z = milnica_film(a,n,fun,slik)
% funkcija Z= milnica_film(a,n,fun,slik,metoda) prikaze film  
% oblikovanja milnice na kvadratu [-a,a]^2, 
% ce so podane robne vrednosti
% n+2 je stevilo delilnih tock na eni koordinatni osi
% fun je funkcija spremenljivke x in a, ki vrne vrednosti na robu
% in sicer matriko s stirimi vrsticami [levo; zgoraj;desno;spodaj]
% parametre maxk in tol podamo funkciji laplace
% slik je stevilo slik v filmu
% uporabimo Jacobijevo iteracijo

% Dolocimo robne vrednosti
x = linspace(-a,a,n+2);
rp = feval(fun,x,a);

Up=rand(n+2);
Up(:,1) = rp(1,:)';
Up(:,n+2) = rp(3,:)';
Up(1,2:n+1) = rp(4,2:n+1);
Up(n+2,2:n+1) = rp(2,2:n+1);
k=slik;
Z=zeros(n+2,n+2,k);
for i=1:k
    [U,k]=laplace_jac(Up,1e-3,i);
    Z(:,:,i)=U;
end
clf
for i=1:slik
    surf(x,x,Z(:,:,i));
    shading flat;
    view([-37,30]);
    pause(0.05);
end
       
