function M=film(lambda0,mu0,parametri,w0,w1,n);
  %parametri=[L,H,ro,g];
  L=parametri(1);
  H=parametri(2);
  M=moviein(n);
  w=linspace(w0,w1,n);
  for j=1:n
    clf;hold on;
    axis([-L/10,L,-L/10,L])
    beta=parametri(4)/(w(j)^2*H);
    l=L/H;
    %X=fsolve(@vrni_g,[lambda0;mu0],optimset('fsolve'),l,beta);  
    %lambda=X(1);
    %mu=X(2);
    [lambda,mu]=newton(@vrni_gj,@vrni_g,lambda0,mu0,l,beta,1e-3)
    [t,u]=ode45(@sistem1,[0;l],[lambda;mu;0;0],[],beta);
    plot(u(:,3)*H,u(:,4)*H,'b');
    %axis equal;
    M(:,j)=getframe;
  end
