% Testroutine zum Programm
% TPSplVal, Blatt 10

disp(' ')
disp('Figur 1 zeigt die Flaeche als Drahtgittermodell');
disp('Figur 2 zeigt die Flaeche schattiert');
disp(' ')
disp('Klicken Sie in Figur 2 auf das Lampensymbol');
disp('und probieren Sie dann die verschiedenen Buttons aus');
disp(' ')

% Daten erzeugen
S = -5:5;
T = -5.5:5.5;

s = (1-eps)*linspace(-3,3,50);
t = (1-eps)*linspace(-2.5,2.5,40);

[X,Y] = meshgrid(s,t);
X = X';
Y = Y';

P = peaks(8);

% Programm aufrufen
Z = TPSplVal(P,S,T,s,t);

% Drahtgittermodell plotten
figure(1)
clf reset
h=mesh(X,Y,Z);
set(h,'EdgeColor','b')
cameratoolbar
figure(1)

% Schattierte Flaeche plotten
figure(2)
clf reset
h=surf(X,Y,Z);
shading interp
%colormap(copper)
cameratoolbar
colorbar
figure(2)