clc
clear all
format long
ind = menu ('Auswahl', 'AWP 1_RungeKutta', 'AWP 1_ode45', 'AWP 2_m2=0.1','AWP 2_m2=0.2', 'AWP 2_m2=0.3', 'AWP 2_m2=0.4','AWP 2_m2=0.5','AWP 2_m2=0.6', 'AWP 2_m2=0.7', 'AWP 2_m2=0.8', 'AWP 2_m2=0.9', 'AWP 2_m2=1')

switch(ind)
    
    case (1)
        
        a = 0;
        b = 3;
        t = [a:b];
        ya = 1;
        h = 0.01;
        f = @(t,y) 2.*t.*y;
        [t, y] = RungeKutta (f, a, b, ya, h);
        

    case (2)
    
        xa = 0;
        xe = 3;
        y0 = 1;
        f = @(x,y) 2.*x.*y;
        [x, y] = ode45 (f, [xa, xe], y0);
    
    case (3)
        
        m1 = 5;
        m2 = 0.1;
        c1 = 10;
        c2 = 1;
        b1 = 0.1;
        b2 = 0.2;
        F1 = 1;
        Omega = 1.414;
        
        M = [m1 , 0; 0, m2];
        C = [c1+c2, -c2; -c2 c2];
        B = [b1+b2, -b2; -b2, b2];
        Fhat = [F1; 0];
        
        a = 0;
        b = 200;
        t = [a:b];
        
        h = 0.01;
        
        ya = [0 0 0 0];
        
        [t,y] = RungeKutta (@(t,y) fd2(t,y,M,B,C,Fhat,Omega),a,b,ya,h);
        
        
        %[t, y] = ode45 (@(t,y) fd(t,y,M,B,C, Fhat, Omega), [0 200], [0 0 0 0]);
        
        x1 = y(1,:);
        x2 = y(2,:);
        
        plot(t, x1, t, x2);
        xlabel('t [s]','FontSize',12);
        ylabel('x [m]','FontSize',12);
        legend('Schwinger', 'Tilger')
        
        
        index = find(t>= 100);
        xSpaet = x1(index);
        amplitude2d = max(abs(xSpaet))
        amplitudeTilger = max(abs(x2(index)))
        maxTilger = max(abs(x2))
        
        
    case (4)
        
        m1 = 5;
        m2 = 0.2;
        c1 = 10;
        c2 = 1;
        b1 = 0.1;
        b2 = 0.2;
        F1 = 1;
        Omega = 1.414;
        
        M = [m1 , 0; 0, m2];
        C = [c1+c2, -c2; -c2 c2];
        B = [b1+b2, -b2; -b2, b2];
        Fhat = [F1; 0];
        
        [t, y] = ode45 (@(t,y) fd(t,y,M,B,C, Fhat, Omega), [0 200], [0 0 0 0]);
        x1 = y(:,1);
        x2 = y(:,2);
        
        plot(t, x1, t, x2);
        xlabel('t [s]','FontSize',12);
        ylabel('x [m]','FontSize',12);
        legend('Schwinger', 'Tilger')
        
      
        index = find(t>= 100);
        xSpaet = x1(index);
        amplitude2d = max(abs(xSpaet))
        amplitudeTilger = max(abs(x2(index)))
        maxTilger = max(abs(x2))
        
        
    case (5)
        
        m1 = 5;
        m2 = 0.3;
        c1 = 10;
        c2 = 1;
        b1 = 0.1;
        b2 = 0.2;
        F1 = 1;
        Omega = 1.414;
        
        M = [m1 , 0; 0, m2];
        C = [c1+c2, -c2; -c2 c2];
        B = [b1+b2, -b2; -b2, b2];
        Fhat = [F1; 0];
        
        [t, y] = ode45 (@(t,y) fd(t,y,M,B,C, Fhat, Omega), [0 200], [0 0 0 0]);
        x1 = y(:,1);
        x2 = y(:,2);
        
        plot(t, x1, t, x2);
        xlabel('t [s]','FontSize',12);
        ylabel('x [m]','FontSize',12);
        legend('Schwinger', 'Tilger')
        
        
        index = find(t>= 100);
        xSpaet = x1(index);
        amplitude2d = max(abs(xSpaet))
        amplitudeTilger = max(abs(x2(index)))
        maxTilger = max(abs(x2))
        
       
     case (6)
         
        m1 = 5;
        m2 = 0.4;
        c1 = 10;
        c2 = 1;
        b1 = 0.1;
        b2 = 0.2;
        F1 = 1;
        Omega = 1.414;
        
        M = [m1 , 0; 0, m2];
        C = [c1+c2, -c2; -c2 c2];
        B = [b1+b2, -b2; -b2, b2];
        Fhat = [F1; 0];
        
        [t, y] = ode45 (@(t,y) fd(t,y,M,B,C, Fhat, Omega), [0 200], [0 0 0 0]);
        x1 = y(:,1);
        x2 = y(:,2);
        
        plot(t, x1, t, x2);
        xlabel('t [s]','FontSize',12);
        ylabel('x [m]','FontSize',12);
        legend('Schwinger', 'Tilger')
        
        
        index = find(t>= 100);
        xSpaet = x1(index);
        amplitude2d = max(abs(xSpaet))
        amplitudeTilger = max(abs(x2(index)))
        maxTilger = max(abs(x2))
        

     case (7)
         
        m1 = 5;
        m2 = 0.5;
        c1 = 10;
        c2 = 1;
        b1 = 0.1;
        b2 = 0.2;
        F1 = 1;
        Omega = 1.414;
        
        M = [m1 , 0; 0, m2];
        C = [c1+c2, -c2; -c2 c2];
        B = [b1+b2, -b2; -b2, b2];
        Fhat = [F1; 0];
        
        [t, y] = ode45 (@(t,y) fd(t,y,M,B,C, Fhat, Omega), [0 200], [0 0 0 0]);
        x1 = y(:,1);
        x2 = y(:,2);
        
        plot(t, x1, t, x2);
        plot(t, x1, t, x2);
        xlabel('t [s]','FontSize',12);
        ylabel('x [m]','FontSize',12);
        legend('Schwinger', 'Tilger')
        
        
        index = find(t>= 100);
        xSpaet = x1(index);
        amplitude2d = max(abs(xSpaet))
        amplitudeTilger = max(abs(x2(index)))
        maxTilger = max(abs(x2))
        
     case (8)
         
        m1 = 5;
        m2 = 0.6;
        c1 = 10;
        c2 = 1;
        b1 = 0.1;
        b2 = 0.2;
        F1 = 1;
        Omega = 1.414;
        
        M = [m1 , 0; 0, m2];
        C = [c1+c2, -c2; -c2 c2];
        B = [b1+b2, -b2; -b2, b2];
        Fhat = [F1; 0];
        
        [t, y] = ode45 (@(t,y) fd(t,y,M,B,C, Fhat, Omega), [0 200], [0 0 0 0]);
        x1 = y(:,1);
        x2 = y(:,2);
        
        plot(t, x1, t, x2);
        plot(t, x1, t, x2);
        xlabel('t [s]','FontSize',12);
        ylabel('x [m]','FontSize',12);
        legend('Schwinger', 'Tilger')
       
        
        index = find(t>= 100);
        xSpaet = x1(index);
        amplitude2d = max(abs(xSpaet))
        amplitudeTilger = max(abs(x2(index)))
        maxTilger = max(abs(x2))
        
       
     case (9)
            
        m1 = 5;
        m2 = 0.7;
        c1 = 10;
        c2 = 1;
        b1 = 0.1;
        b2 = 0.2;
        F1 = 1;
        Omega = 1.414;
        
        M = [m1 , 0; 0, m2];
        C = [c1+c2, -c2; -c2 c2];
        B = [b1+b2, -b2; -b2, b2];
        Fhat = [F1; 0];
        
        [t, y] = ode45 (@(t,y) fd(t,y,M,B,C, Fhat, Omega), [0 200], [0 0 0 0]);
        x1 = y(:,1);
        x2 = y(:,2);
        
        plot(t, x1, t, x2);
        plot(t, x1, t, x2);
        xlabel('t [s]','FontSize',12);
        ylabel('x [m]','FontSize',12);
        legend('Schwinger', 'Tilger')
       
        
        index = find(t>= 100);
        xSpaet = x1(index);
        amplitude2d = max(abs(xSpaet))
        amplitudeTilger = max(abs(x2(index)))
        maxTilger = max(abs(x2))
       
        
    case (10)
            
        m1 = 5;
        m2 = 0.8;
        c1 = 10;
        c2 = 1;
        b1 = 0.1;
        b2 = 0.2;
        F1 = 1;
        Omega = 1.414;
        
        M = [m1 , 0; 0, m2];
        C = [c1+c2, -c2; -c2 c2];
        B = [b1+b2, -b2; -b2, b2];
        Fhat = [F1; 0];
        
        [t, y] = ode45 (@(t,y) fd(t,y,M,B,C, Fhat, Omega), [0 200], [0 0 0 0]);
        x1 = y(:,1);
        x2 = y(:,2);
        
        plot(t, x1, t, x2);
        plot(t, x1, t, x2);
        xlabel('t [s]','FontSize',12);
        ylabel('x [m]','FontSize',12);
        legend('Schwinger', 'Tilger')
        
        index = find(t>= 100);
        xSpaet = x1(index);
        amplitude2d = max(abs(xSpaet))
        amplitudeTilger = max(abs(x2(index)))
        maxTilger = max(abs(x2))
        
      
    case (11)
            
        m1 = 5;
        m2 = 0.9;
        c1 = 10;
        c2 = 1;
        b1 = 0.1;
        b2 = 0.2;
        F1 = 1;
        Omega = 1.414;
        
        M = [m1 , 0; 0, m2];
        C = [c1+c2, -c2; -c2 c2];
        B = [b1+b2, -b2; -b2, b2];
        Fhat = [F1; 0];
        
        [t, y] = ode45 (@(t,y) fd(t,y,M,B,C, Fhat, Omega), [0 200], [0 0 0 0]);
        x1 = y(:,1);
        x2 = y(:,2);
        
        plot(t, x1, t, x2);
        plot(t, x1, t, x2);
        xlabel('t [s]','FontSize',12);
        ylabel('x [m]','FontSize',12);
        legend('Schwinger', 'Tilger')
        
        
        index = find(t>= 100);
        xSpaet = x1(index);
        amplitude2d = max(abs(xSpaet))
        amplitudeTilger = max(abs(x2(index)))
        maxTilger = max(abs(x2))
        
       
    case (12)
       
        m1 = 5;
        m2 = 1;
        c1 = 10;
        c2 = 1;
        b1 = 0.1;
        b2 = 0.2;
        F1 = 1;
        Omega = 1.414;
        
        M = [m1 , 0; 0, m2];
        C = [c1+c2, -c2; -c2 c2];
        B = [b1+b2, -b2; -b2, b2];
        Fhat = [F1; 0];
        
        [t, y] = ode45 (@(t,y) fd(t,y,M,B,C, Fhat, Omega), [0 200], [0 0 0 0]);
        x1 = y(:,1);
        x2 = y(:,2);
        
        plot(t, x1, t, x2);
        plot(t, x1, t, x2);
        xlabel('t [s]','FontSize',12);
        ylabel('x [m]','FontSize',12);
        legend('Schwinger', 'Tilger')
       
        
        index = find(t>= 100);
        xSpaet = x1(index);
        amplitude2d = max(abs(xSpaet))
        amplitudeTilger = max(abs(x2(index)))
        maxTilger = max(abs(x2))
       
end

amp = tilger(param,ind);
param = [m2,c2,b2];
max(abs(eta(:,1)))