omega_u=linspace(1,7,1000);
%b=linspace(10,100,1000);


F=1900;
m=10;
k=100;
omega_0=sqrt(k/m);
b_c=2*sqrt(k*m)
b=(b_c)/10;
b_sim=linspace(1,b_c-1,1000);
iplot=3;
N=n1;
%A=A_functie(F,m,omega_u,omega_0,b);


%plotten van de 10 grafieken met verschillende b
%-------------------------------------------------

Amaximas=zeros(n1,3); % 1ekolom = omega_u, 2ekolom = A_max, 2ekolom = b's

    A=A_functie(F,m,omega_u,omega_0,b);
    phi=fase_functie(b,omega_u,m,omega_0);
    V=A.*omega_u;
    
    figure(1);
    plot(omega_u,A,'b',omega_u,phi,'r',omega_u,V/3,'g');
    zoom on;
    ylabel('Amplitude (blauw), Faseverschil (rood) en Snelheid (groen)');
    xlabel('\omega_u (pulsatie van aandrijvende kracht)');
    title('Amplitude- en fasecurve');
    zoom on;
    hold on;


saveas(gcf,'Extra','tif');		%opslaan van figuur voor import in andere programma's
%saveas(gcf,'BijkomendeOpg1','fig');