function dh=Fvecwithproduction_nop(x,y)
global rho sigma sigmasqr lambda gamma A phi delta kappa etah etal eentry m maversion

q=y(1,:);qprime=y(2,:);%q=y(1,:);qprime=y(2,:);
theta=max(q./x,1/(1-lambda));
i=(q-1)/kappa;c=A-delta-i-kappa.*(i.^2)./2;
F=q./x-m.*theta.*(qprime);
G=c.*kappa.*F+gamma.*q.*qprime.*m.*(1-theta).*q;
eta=etal.*(x<eentry)+etah.*(x>=eentry);

sigmae=x.*q.*sigma.*(m.*theta-1)./(q-x.*m.*theta.*(qprime));

%%up to here

%a11=pprime.*c.*kappa./G.*(-m.*(1-theta).*gamma.*0.5.*q.*sigmae.^2./c./kappa.*(p+q)+m.*theta.*0.5.*sigmae.^2)...
%    +p.*gamma.*(q.*qprime.*m.*theta.*0.5.*sigmae.^2+F.*sigmae.^2*0.5.*q)./G;
%a12=pprime.*(c.*kappa.*m.*theta.*0.5.*sigmae.^2)./G+0.5.*sigmae.^2+p.*gamma.*q.*qprime.*m.*theta.*0.5.*sigmae.^2./G;
a=qprime.*(c.*kappa./G.*(-m.*(1-theta).*gamma.*0.5.*q.*sigmae.^2./c./kappa.*q+m.*theta.*0.5.*sigmae.^2))...
    +0.5.*sigmae.^2+gamma.*q./G.*(q.*qprime.*m.*theta.*0.5.*sigmae.^2+F.*q.*sigmae.^2.*0.5);
%a22=qprime.*(c.*kappa.*m.*theta.*0.5.*sigmae.^2)./G+gamma.*(q.^2).*qprime.*m.*theta.*0.5.*sigmae.^2./G;
mm=(m.*(1-theta).*(rho+gamma.*i-gamma.*0.5.*(qprime.^2.*sigmae.^2./(c.*kappa))...
    -0.5.*gamma.*(1+gamma).*(c.*sigma-q.*qprime.*sigmae./kappa).^2./c.^2)-i-eta).*q...
    +m.*theta.*(-q.*delta+A);
nn=(rho+gamma.*i).*c.*kappa.*F...
    -gamma.*q.*qprime.*((-i-eta).*q+m.*theta.*(-q.*delta+A))...
-0.5.*gamma.*F.*qprime.^2.*(sigmae.^2)-0.5.*gamma.*(1+gamma).*(F.*kappa.*(c.*sigma-q.*qprime.*sigmae./kappa).^2)./c;
%nn2=(rho+i-delta).*(A-i-kappa.*(i.^2)./2).*kappa.*((p+q)./x-theta.*(qprime+pprime))...
%    -q.*qprime.*(((1-theta).*delta-i-eta).*(p+q)+theta.*(A./theta+(p-(1-theta)./theta).*i-(1-theta)./theta*0.5*kappa.*i.^2))...
%-0.5*((p+q)./x-theta.*(qprime+pprime)).*qprime.^2.*(sigmae.^2)...
%-(((p+q)./x-theta.*(qprime+pprime)).*kappa.*((A-i-kappa.*(i.^2)./2).*sigma-q.*qprime.*sigmae./kappa).^2)./(A-i-kappa.*(i.^2)./2);

b=(q.*sigma+qprime.*sigmae).*sigma.*maversion.*theta.*((q-x.*qprime)./(x.*F))...
    -qprime.*c.*kappa./G.*(mm)...
    -A+q.*delta+q./G.*nn;
q2prime=b./a;
dh=[y(2,:)
    q2prime];