
% Benytter seg av trim rutinen til å finne steadystate-tilstand

%clear all
q=ones(12,1);
%optu0;
%opty0;
%basetilstand;
i=1

[X,U,Y,DX]=trim('opttelin',x0,u0,y0,[],[5 9 12],[5 7 8 9 12 15 31 38 39])

for j=1:20,
 
  x0=X
 if i>1
  sum(abs(((u0-U)./U)))
 end
  i=i+1
  [X,U,Y,DX]=trim('opttelin',x0,u0,y0,[],[5 9 12],[5 7 8 9 12 15 31 38 39])
         
end

X-x0
Y-y0
U-u0

%Linearisere modellen om stasonærpunktet

[A,B,C,D]=linmod('opttelin',X,U) 



