function [upr_vec,upreps_vec,Fdef_vec,Fdefeps_vec] = ...
    getDefaultFcts(epshat,epsstar_nd,epsbar_d,A,B_nd,B_d,params)

gaMMa=params.gaMMa;
m=params.logmueps;
s=params.logsigeps;
nodes=params.nodes;
weights=params.weights;

f_epshat_nd=integEpsGL(epshat,epsstar_nd,A,B_nd,m,s,gaMMa,nodes,weights);
f_epsbar_d=integEpsGL(0,epsbar_d,A,B_d,m,s,gaMMa,nodes,weights);

upr_vec=[f_epshat_nd(:,2),f_epsbar_d(:,2)];
upreps_vec=[f_epshat_nd(:,3),f_epsbar_d(:,3)];

Fdef_vec=zeros(size(A,1),4);
Fdef_vec(:,1)=logncdf(epsstar_nd,m,s)-logncdf(epshat,m,s);
Fdef_vec(:,2)=logncdf(epsbar_d,m,s);
Fdef_vec(:,3)=1-logncdf(epsstar_nd,m,s);
Fdef_vec(:,4)=logncdf(epshat,m,s)-logncdf(epsbar_d,m,s);

Fdefeps_vec=zeros(size(A,1),2);
s2 = s.^2;
Fdefeps_vec(:,1)=exp(m + s2/2) .* (normcdf((log(epsstar_nd) - m - s2)./s) ...
                                  - normcdf((log(epshat) - m - s2)./s));
Fdefeps_vec(:,2)=exp(m + s2/2) .* normcdf((log(epsbar_d) - m - s2)./s);

end