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

c=0.5;
u=@(x,t) 1-2*cos(3*x).*cos(3*c*t); %solucio a visualitzar

a=0; b=pi; hx=0.01; x=a:hx:b;
xmin=a-0.3; xmax=b+0.3;
umin=-1.5; umax=3.5;
t0=0; tf=5; ht=0.05;

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

if( guardaQ )
  close(V)
end
