guardaQ=0; %valors: 1 o 0, segons si volem guardar fitxer de video o no
if( guardaQ )
  V=VideoWriter('edpdalembert1');          %fitxer AVI
  %V=VideoWriter('edpdalembert1','MPEG-4'); %fitxer MP4
  V.FrameRate=10;
  open(V)
end

u=@(x,t) cos(x)-sin(t)-cos(x).*cos(t)+t; %solucio a visualitzar

a=-10; b=10; hx=0.01; x=a:hx:b;
umin=-2; umax=16;
t0=0; tf=12; ht=0.05;

for t=t0:ht:tf
  uu=@(x) u(x,t);
  z=uu(x);
  plot(x,z,'b','LineWidth',2.5)
  axis([a,b,umin,umax])
  grid
  hold on
  plot([a,b],[0,0],'k')
  plot([0,0],[umin,umax],'k')
  legend('u(x,t)')
  text(-9,-1,sprintf('t=%.2f',t))
  hold off
  drawnow
  pause(0.05)
  if( guardaQ )
    writeVideo(V,getframe)
  end
end

if( guardaQ )
  close(V)
end
