% Funktionen mit "handle" definieren"
f=@(x) sin(x);   % fiel mir grad so ein
g=@(x) cos(x);   % Ableitung, selber im Kopf bestimmt

x=-pi:pi/100:pi; % Definitionsbereicht, Knotenpunkte
y=f(x);          % Funktionswerte - zusammen mit x - Wertetabelle

% figure(1)      % erstes Abbilung, einfach malen und Achsen dazu
% plot(x,y)
% axis tight     % Koordinatensystem so knapp wie möglich
% line(x,0*x)    % Linienbefehl, es müssen genauso viel y- wie x-Werte sein
% line(0*y,y)
% 
% figure(2)      % noch mehr Spielchen zum "Nicen" von Abbildungen  
% x=pi:pi/100:3*pi;
% y=sin(x);
% plot(x,y,'m','linewidth',3,'marker','*')
% axis tight
% axis equal     % gleiche Einheiten auf beiden Achsen 
% xlabel('x')    % Achsenbeschriftungen
% ylabel('y','rotation',0)
% 
% figure(3)      % diskretisierte Ableitungen
h=x(2)-x(1);
dydx=diff(y)/h;
% plot(x(1:200),dydx)
% axis tight

figure(4)       % mehr Effekte 
plot((x(1:200)+x(2:201))/2,dydx,'g', x(1:200),dydx,'r','linewidth',3)
set(gca,'xlim',[x(1) x(end)],'ylim',[-1 1])
set(gca,'plotboxaspectratio',[2*pi 2 1],'ylim',[-1 1])
set(gcf,'position',[50 250 1450 500])
set(gca,'fontsize',20)
xlabel('x')
ylabel('y','rotation',0)
set(line(x,x/2),'color','k','linewidth',3,'linestyle','--')

% lösen einer nichtlinearen Gleichung mit einer Unbekannten
x0=fzero(@(x) g(x)-x/2, 2);
y0=g(x0);
set(legend('cos(x)','cos(x+h/2)','x/2'),'fontsize',16) % Legende erstellen

% Schnittpunkt markieren, Linie aus genau einem markierten Punkt
set(line(x0,y0),'marker','+','markersize',20,'color','c','linewidth',5)

% Ausgabe als png-Grafikdatei
print -dpng sinusitis

