% Max-periaate (Vrt. Maple: L6.mws)
close all
u=inline('t+1/2*x.^2','x','t')   % k=1, C=0
x=linspace(-1,1,30);
t=linspace(0,1,20);
[X,T]=meshgrid(x,t);
surfc(x,t,u(X,T))
colorbar

U=u(X,T);  
figure
plot(x,U)  % Aikasiivut
figure
plot(t,U)  % "Paikkasiivut"


umax=max(max(U))

X(U==umax)
%ans =

 %   -1
 %    1


T(U==umax)

%ans =

 %    1
 %    1

 % Tässä max-periaatteen ilmentymä: max pisteissä (-1,1) ja (1,1).    
 
 % Vastaavasti min:
 umin=min(U(:))   % Toinen tapa, jonoutetaan U vektoriksi.
 xmin=X(U==umin)
 tmin=T(U==umin)
 
