function Jud = Judderiv2_2est(spec1,step1, step21, step22, step23, step24, step25, step26, step_pf, step_pt, step_pb, bests,xDs,xBs)
global Col Specs 
    load_bests_2est
    Col.ColumnFlowsheet.Run
    MV1 = Specs.Item(spec1).CurrentValue;
    Specs.Item(spec1).GoalValue = MV1 * (1 + step1);
    Col.ColumnFlowsheet.Run
    J_1u = Goal_2est(xDs,xBs);
    
    load_bests_2est
    Col.ColumnFlowsheet.Run
    J_1d = DistGain_2est([],[],[],step21, step22, step23, step24, step25, step26, step_pf, step_pt, step_pb, 1);
    Jud = DistGain_2est(J_1u, J_1d, step1, step21, step22, step23, step24, step25, step26, step_pf, step_pt, step_pb, 0);
