function [ output_args ] = maxima( data )
%MAXIMA Bestimmung der Maxima mittels Gauß-Lorentz-Fit
%1. Import der Molekül-Dateien



%% Import der Molekül-Dateien

A=1;n=0;initial_parameters_=0;length_=0;

while A~=0
n=n+1;

[filename, pathname] = uigetfile({'tab/*.xlsx'},'Öffnen');          %Auswahl der Daten
fullpathname = strcat(pathname,filename); 


v = genvarname('initial_parameters_', who);
   eval([v ' = xlsread(fullpathname)']);

w = genvarname('length_', who);
   eval([w ' = length(xlsread(fullpathname))']);


    % Construct a questdlg with two options
    choice = questdlg('Load another Molecule?', ...
    	'Answer', ...
    	'Yes','No','No');
    % Handle response
    switch choice
        case 'Yes'
        A=A+1;
        
        case 'No'
        A=0;      
    end

end

summe = 0;
for z=1:n
summe = summe + eval(['length_',num2str(z)]);
end


% n = Anzahl der Dateien
% initial_parameters_1 bis initial_parameters_n Dateien
% length_1 bis length_n Länge der Dateien
% 'summe': Anzahl theoretischer Peaks

%% Fit Funktion

Fun_=0;Fun=0;

for z=1:summe    
fun = [['a',num2str(z)],'*exp(-((x-',['b',num2str(z)],')/',['c',num2str(z)],').^2)'];
q = genvarname('Fun_', who);
   eval([q '= fun']);
   

   Fun=[[Fun,'+'],eval(['Fun_',num2str(z)])];    
end


%%

Var=[];
for z=1:summe
   
Var = [Var,['a',num2str(z),', '],['b',num2str(z),', '],['c',num2str(z)],', '];
% eval([['a',num2str(z)] '= []']);
% eval([['b',num2str(z)] '= []']);
% eval([['c',num2str(z)] '= []']);
end
Var =[Var, 'x'];


g= ['fittype( @(',Var,') ' ,Fun,' );'];


%%






end