function [res,xnew]    = callsteady(newpar,fullpar);

I            = fullpar(1);
xguess       = fullpar(2:3*I+I*I);
iparam       = fullpar(I*I+3*I+1:13+14*I+8*I*I);
datastats    = fullpar(14+14*I+8*I*I:13+15*I+10*I*I);
iparind      = fullpar(14+15*I+10*I*I:13+16*I+12*I*I);
nind         = fullpar(14+16*I+12*I*I);
dataind      = fullpar(15+16*I+12*I*I:length(fullpar));

pnew         = iparam;
ii           = iparind(dataind);
pnew(ii)     = newpar;
xnew         = secant('steady',xguess,pnew,1e-5,100);
[r,modstats] = steady(xnew,pnew);
%res          = datastats(dataind)./modstats(dataind)-1;

res          = datastats(dataind)-modstats(dataind);


%disp(sum(res))
%disp([datastats(dataind),modstats(dataind),res])
%disp([datastats,modstats]);
