% Delivered to the U.S. Government in Accordance with Contract FA8650-10-C-7080.
%
%
% (c) 2011 United Technologies Corp.
% Developed with the sponsorship of the Defense Advanced Research Projects Agency (DARPA).
% Permission is hereby granted, free of charge, to any person obtaining a copy of this data,
% including any software or models in source or binary form, as well as any drawings, 
% specifications, and documentation (collectively "the Data"), to deal in the Data without 
% restriction, including without limitation the rights to use, copy, modify, merge, publish,
% distribute, sublicense, and/or sell copies of the Data, and to permit persons to whom the Data
% is furnished to do so, subject to the following conditions:
%  
% The above copyright notice and this permission notice shall be included in all copies or
% substantial portions of the Data.
%  
% THE DATA IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING
% BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
% NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS, SPONSORS, DEVELOPERS, CONTRIBUTORS, OR
% COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
% OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE DATA OR THE
% USE OR OTHER DEALINGS IN THE DATA.

% Authors:  Fernando Rodriguez                      John Finn
%           Sr Engineer/Scientist                   Contractor
%           Power Electronics Systems               Embedded Systems & Networks
%           United Technologies Research Center     United Technologies Research Center
%           411 Silver Lane                         411 Silver Lane
%           East Hartford,CT 06108                  East Hartford,CT 06108
%           rodrigf@utrc.utc.com                    ccfinnje@utrc.utc.com
%           tel: (860) 610-7102                     tel: +1 860-610-7696
%%

close all
clear all


nl = 36;    
p = 1;
Gen_Bus = 3;

[A, B, Cr, Ci, B_No] = statespace_formulation_rev08_orig(nl);

loc = -1;
for i = 1: length(B_No)
    if B_No(i) == Gen_Bus
        loc = i;
    end
end
D = zeros(1,length(B_No));
D(loc) = 1;

vrms = 115;
f = 400;
Ki = 20;
Kp = .1;
dec = 1;

start = 0;
stop = 1;
reltol = 1e-8;
abstol = 1e-8;
model = 'solverABCD_man_file_C_feed.mdl';

delete('output.mat')
tic
sim(model,'StartTime', num2str(start),...
        'StopTime', num2str(stop), ...
         'Abstol', num2str(abstol), ...
         'Reltol', num2str(reltol));

rt = toc;

fprintf('\n\nMean Runtime: %f seconds\n\n', rt);

close all
if p
    plot_data(nl, B_No)
end