When I add a non-binding constraint to a maximization problem the result changes. I don't understand why.

Below you can find the code that I am using.

For the maximization A, H = 6837.66
For the maximization B (equal to A + constraint of H less or equal to 14000), H = 5716.08

If the solution in A is already below 14000, the solution should not change. Why is this happening?

Thank you!

```
$Assumptions -> {{β, θ, α, r, p, ϵ, σ, δ} > 0, Element [{H, β, θ, α, r, p, ϵ, σ, δ}, Reals], {β, θ, α, r, p, ϵ, σ, δ} < 1, H > 0, ϵ > σ, α < r};

Y = 1000;
β = 0.95;
θ = 10;
α = 0.05;
r = 0.05;
ϵ = 1.2;
σ = 0.2;
δ = 0.3;
p = 0.51;

g1[H_] := Max[0, ((((1 - δ)*Y)/(r + α)) - H)]*((((1 - δ)*Y)/(r + α)) - H)^(-1)

h1[H_] :=  Min[0, ((((1 - δ)*Y)/(r + α)) - H)]*((((1 - δ)*Y)/(r + α)) - H)^(-1)

A = Maximize[{(1 + β * (1 - 
          p))*(((Y - r*H)^(1 - 1/ϵ) + (θ*H)^(1 - 
             1/ϵ))^((1 - 1/σ)/(1 - 
            1/ϵ))) /(1 - 1/σ) +  β*p *
     g1[H] *(((Y - (r + α)*H)^(1 - 1/ϵ) + (θ*
             H)^(1 - 1/ϵ))^((1 - 1/σ)/(1 - 
            1/ϵ)) ) /(1 - 1/σ) + β*p *
     h1[H] *(((δ*Y)^(1 - 1/ϵ) + (θ*5000)^(1 - 
             1/ϵ))^((1 - 1/σ)/(1 - 
            1/ϵ)) ) /(1 - 1/σ), 
   H ⩾ 0, (Y - r*H) ⩾ 0}, {H}]



B = Maximize[{(1 + β * (1 - 
          p))*(((Y - r*H)^(1 - 1/ϵ) + (θ*H)^(1 - 
             1/ϵ))^((1 - 1/σ)/(1 - 
            1/ϵ))) /(1 - 1/σ) +  β*p *
     g1[H] *(((Y - (r + α)*H)^(1 - 1/ϵ) + (θ*
             H)^(1 - 1/ϵ))^((1 - 1/σ)/(1 - 
            1/ϵ)) ) /(1 - 1/σ) + β*p *
     h1[H] *(((δ*Y)^(1 - 1/ϵ) + (θ*5000)^(1 - 
             1/ϵ))^((1 - 1/σ)/(1 - 
            1/ϵ)) ) /(1 - 1/σ), 
   H ⩾ 0, (Y - r*H) ⩾ 0, 
   H <= 14000}, {H}]

```