function [filename3]=v25_call_solving(filename3,solvepar)
%Load the specification file
load([filename3 '1']);
%Check if all or just a limited number of the specifications should be used
if ~isfield(solvepar,'startat'),solvepar.startat = 1;end
if ~isfield(solvepar,'endat'),solvepar.endat = num;end

if solvepar.startat>1;load(['results/' saveasfilename]);load([filename3 '1']);end;

%Solving only:
for i=solvepar.startat:solvepar.endat
    %Display which number is running
    disp(['Solving scenario ' num2str(i) ' out of ' num2str(num) ', i.e. ' num2str(round(i/num*1*100)) '%'])
    if i>1; load([filename3 num2str(i)]); end;%Always load the appropriate dataset
    %Call solving routine
    [interp(i,:),model(i,:),s(i).x,x(i).x,z(i).x]=v25_solve(filename,synchron,varnum,usesolver,loop,method,sim1,sim2,sim_start,shock_nodes,approx_order,state_min,state_max,t_convergence,Mu,sigma1,par);
end;
clear solvepar;
save results/results_all_last;save(['results/' saveasfilename]);%Save results
