%#######################################################################
% Octave-Skript zum Erzeugen eines simplen 1D und eines 2D-Plots
%#######################################################################

% alte Variablen wegputzen
clear

% alle offenen Plot-Fenster schliessen
close all
% 'close' alleine schliesst nur das zuletzt geoeffnete Fenster

%#################################################
% Plot 1: Wir wollen ein Polynom auf einem Gitter
% auswerten.
%#################################################

% erzeuge ein aequidistantes Gitter zwischen 1.2 und 2.8 mit Schrittweite 0.01
x=1.2:0.01:2.8;

% Auswerten eines Polynoms an dern Gitterpunkten
% Man beachte: Die folgende Auswertung kommt ohne Schleifen aus!
y = 4.*(x-1.4).*(x-2).*(x-2.2).*(x-2.5);

% Zeichnen der Funktionswerte. Standardmaessig werden diese linear verbunden
%plot(x,y)
%plot(x,y,'r')
%plot(x,y,'b')
plot(x,y,'*k-')
% plot(x,y,'k')

% aktuelles Plotfenster aktiv halten, damit nachfolgende Plot-Befehle wieder
% in dieses Fenster gezeichnet werden
hold on
plot(y,x)

% ab sofort sollen fuer neue Plots auch wieder neue Fenster geoeffnet werden
hold off


% Formatieren der Achsen
% Syntax: axis ([xStart xEnde yStart yEnde]); (2 Argumente pro Dimension).
axis([1.2 2.8 -0.2 0.2]);


% #################################################
% # Plot 2: Wir wollen die Funktion sin(x)*cos(y)
% # auf dem Quadrat [-2*pi,2*pi]x[-1*pi,2*pi]
% # auswerten.
% ##################################################

% alle Fenster zu
%close all;

% Wir erzeugen diesmal gleich zu Beginn ein leeres Plot-Fenster.
figure

% Gitter erzeugen
% Wir unterteilen den Bereich von -2*pi bis 2*pi in 50 gleiche Teile.
x=linspace(-2*pi,2*pi,50);

SIN=sin(x);
COS=cos(x);

[X Y] = meshgrid(x,x);

% fuer Oberflaechenplots
mesh(X,Y,SIN'*COS)
%contour(X,Y,SIN'*COS)

figure

% eine Spiralkurve Kurve in 3D
t = 0:0.1:10*pi;
r = linspace (0, 1, numel (t));
z = linspace (0, 1, numel (t));
plot3 (r.*sin(t), r.*cos(t), z);
%hold on
