% input_file_name = 'sample_dahlquist_data_fixed_P_free_b';
input_file_name = 'full_schade_data';
igraph          = 1;
iestimate       = 1;
[output]        = general_parameter_estimation_function(input_file_name,igraph,iestimate);

% figHandles  = findobj('Type','figure');
% nfig        = max(figHandles);
% 
% for ii = 3:nfig
%     
%     eval(['figure(' num2str(ii) ')']);
%     eval(['print -djpeg figure_' num2str(ii-2)]);
%     eval(['print -dill figure_' num2str(ii-2)]);
%     
% end