function hyp1
disp('Initial Value Problem of Hyperbolic PDE')
disp('Model Equation u_t + c u_x = 0')
disp('Initial condition u(x,0) = exp(-50*x^2), x in [-1,2] ')
c = input('Wave speed c ?  (default = 1)  ');
if isempty(c)
   c = 1.;
   end

time = input('Stopping time ?  (default = 1)  ');
if isempty(time)
   time = 1.;
   end

dx = 1./50;
x = -1:dx:2.;
dt = time/20;
t = 0:dt:time;
%
clf
for i=1:length(t)
    u = exp(-50.*(x-c*t(i)).^2);
    plot(x,u)
%    disp(['Solution at time t = ',num2str(t(i))]);
    title(['Solution at time t = ',num2str(t(i))]);
    M(i) = getframe;
%    fname = ['tmp.',num2str(i)]
%    print('-djpeg',fname)
end
movie(M,3)

