hold off                          
x=0:.01:1;                        
plot(x,2*exp(2*x),";2*exp(2x);")  % analytische Lösung
hold on                      
h=1/4;                            % Wahl der Schrittweite
xh=0:h:1;                         
u=2;                              % Anfangswert
for i=1:length(xh)-1              
u=[u 1/(1-2*h)*u(length(u))];     
end                               
plot(xh,u,";Differenzen h=1/4;")                        
                                   %Ausgabe in diff.ps
gset key left                    
gset term postscript             
gset output "diff.ps"            
replot