function house=getModel

fileName = 'houseModel';
order = [7 3 5]; %ny nu nx

%get parameters and initial states
[x0,par]=getData;

% Initial states
initialStates = x0;

flag=0;

% Constructing model object
Ts = 0;
house = idnlgrey(fileName,order,flag,initialStates,Ts,'Name','House heating model');

inputName = {'Heat to the floor','Heat to the room','Air flow'};
inputUnit = {'KW','KW','Kg/s'};
outputName = {'Temperature floor','Temeprature room','air mass','Co2 Level','Accumulated Energy Cost','Outdoor Temperature','Price'};
outputUnit = {'K','K','Kg','%','NOK','K','NOK/KW'};
stateName =  {'Temperature floor','Temeprature room','air mass','Co2 Level','Accumlated Energy Cost'};
stateUnit = {'K','K','Kg','%','NOK'};

house.Algorithm.SimulationOptions.Solver='ode23s';
house.Algorithm.SimulationOptions.RelTol=1e-6;
house.Algorithm.SimulationOptions.AbsTol=1e-6;

set(house, 'TimeVariable','s')
set(house, 'InputName',inputName, 'InputUnit',inputUnit);
setinit(house, 'Name',stateName)
setinit(house,'Unit',stateUnit)
setinit(house, 'Fixed',{true true true true true});
set(house, 'OutputName',outputName,'OutputUnit',outputUnit);
%setpar(nlcstr, 'Fixed',{true true true true true true true true true true true true true true true});
 setpar(house,'Name','param')
 setpar(house,'Unit','')
 setpar(house, 'Fixed',{true});

 
present(house)