clear ;close;clc;

global m d rho                         %globale Variablen
prompt={'Pendelmasse [kg]'
        'd [1/s]'
        'rho [N/mm]'};
def={'3.841e-27','15','0'};
title='Eingabe-Fenster';
lineNo=1;

answer=inputdlg(prompt,title,lineNo,def);



m=str2double(answer{1});   %masse in [g]?????
d=str2double(answer{2});           %Frequenz in [1/s]??????
rho=str2double(answer{3});          %Dämpfungskonstante in [N/mm]????


%%          EINGABE DER ANFANGSBEDINGUNGEN

prompt   = {'Pendel 1 Anfangsauslenkung [nm]','Pendel 1 Anfangsgeschwindigkeit [m/s]'...
            'Pendel 2 Anfangsauslenkung [nm]','Pendel 2 Anfangsgeschwindigkeit [m/s]'...
            'Pendel 3 Anfangsauslenkung [nm]','Pendel 3 Anfangsgeschwindigkeit [m/s]'...
            'Pendel 4 Anfangsauslenkung [nm]','Pendel 4 Anfangsgeschwindigkeit [m/s]'...
            'Pendel 5 Anfangsauslenkung [nm]','Pendel 5 Anfangsgeschwindigkeit [m/s]'...
            'Pendel 6 Anfangsauslenkung [nm]','Pendel 6 Anfangsgeschwindigkeit [m/s]'};
                            

name     = 'Anfangsbedingungen 1:';                                                                                   % Überschrift des Dialogfensters
numlines = 1;                                                                                                       % Menge der Eingabezeilen
def      = {'1','0','0','0','0','0','0','0','0','0','0','0'};   % Default-Einstellungen der AB                                         

answer2  = inputdlg(prompt,name,numlines,def); 


%%          ZUWEISEN DER VARIABLEN DER AB
U1x_0    = (str2double(answer2{1}));%*10^-9;                               
U1v_0    = str2double(answer2{2});                                                       
U2x_0    = (str2double(answer2{3}));%*10^-9;                               
U2v_0    = str2double(answer2{4});                                                      
U3x_0    = (str2double(answer2{5}));%*10^-9;                               
U3v_0   = str2double(answer2{6});  
U4x_0    = (str2double(answer2{7}));%*10^-9;                               
U4v_0    = str2double(answer2{8});
U5x_0    = (str2double(answer2{9}));%*10^-9;                               
U5v_0   = str2double(answer2{10});
U6x_0    = (str2double(answer2{11}));%*10^-9;                               
U6v_0    = str2double(answer2{12});

%%




prompt   = {'Pendel 7 Anfangsauslenkung [nm]','Pendel 7 Anfangsgeschwindigkeit [m/s]'...
            'Pendel 8 Anfangsauslenkung [nm]','Pendel 8 Anfangsgeschwindigkeit [m/s]'...
            'Pendel 9 Anfangsauslenkung [nm]','Pendel 9 Anfangsgeschwindigkeit [m/s]'...
            'Pendel 10 Anfangsauslenkung [nm]','Pendel 10 Anfangsgeschwindigkeit [m/s]'...
            'Pendel 11 Anfangsauslenkung [nm]','Pendel 11 Anfangsgeschwindigkeit [m/s]'...
            'Pendel 12 Anfangsauslenkung [nm]','Pendel 12 Anfangsgeschwindigkeit [m/s]'};
            
                            

name     = 'Anfangsbedingungen 2:';                                                                                   % Überschrift des Dialogfensters
numlines = 1;                                                                                                       % Menge der Eingabezeilen
def      = {'0','0','0','0','0','0','0','0','0','0','0','0'};   % Default-Einstellungen der AB                                              

answer3  = inputdlg(prompt,name,numlines,def); 




U7x_0    = (str2double(answer3{1}));%*10^-9;                               
U7v_0    = str2double(answer3{2});
U8x_0    = (str2double(answer3{3}));%*10^-9;                               
U8v_0    = str2double(answer3{4});
U9x_0    = (str2double(answer3{5}));%*10^-9;                               
U9v_0    = str2double(answer3{6});
U10x_0    = (str2double(answer3{7}));%*10^-9;                               
U10v_0    = str2double(answer3{8});
U11x_0    = (str2double(answer3{9}));%*10^-9;                               
U11v_0    = str2double(answer3{10});
U12x_0    = (str2double(answer3{11}));%*10^-9;                               
U12v_0    = str2double(answer3{12});



U_0=[U1x_0 U1v_0 U2x_0 U2v_0 U3x_0 U3v_0 U4x_0 U4v_0 U5x_0 U5v_0 U6x_0 U6v_0 U7x_0 U7v_0 U8x_0 U8v_0 U9x_0 U9v_0 U10x_0 U10v_0 U11x_0 U11v_0 U12x_0 U12v_0];  

Zeit=[0,50];     

Genauigkeit=odeset('RelTol',1e-10,'AbsTol',1e-10);

[t,U]=ode45('oszitotal_dgl',Zeit,U_0,Genauigkeit);