### tp108r.apm.m4

```Model tp108r
! Source version 1

! The present file has to be drawn through the m4 macro processor
! at first, with or without `-Drevisedhs'. With the macro
! defined, the feasible domain is reduced in comparison with the H+S
! one such that the solution becomes unique.

ifdef(`revisedhs',`define(`stricths',0)',`define(`stricths',1)')

Parameters
mfcorrhs = 0
End Parameters

Variables
x[1:8] = 1
x[9:9] = 1, >= 0
obj
End Variables

Intermediates
c[ 1] = 1 - x^2 - x^2
c[ 2] = 1 - x^2
c[ 3] = 1 - x^2 - x^2
c[ 4] = 1 - x^2 - (x - x)^2
c[ 5] = 1 - (x - x)^2 - (x - x)^2
c[ 6] = 1 - (x - x)^2 - (x - x)^2
c[ 7] = 1 - (x - x)^2 - (x - x)^2
c[ 8] = 1 - (x - x)^2 - (x - x)^2
c[ 9] = 1 - x^2 - (x - x)^2
c = x*x - x*x
c = x*x
c = (-1)*x*x
c = x*x - x*x
mf = (-1/2)*(x*x - x*x + &
x*x - x*x + &
x*x - x*x)
mfcorrr = (3*x - 2*x)^2 + (3*x - 2*x)^2
End Intermediates

Equations
c[1:13] >= 0

! the problem is pretty ill-posed; though the best known
! objective is quite stable the solution is capricious
obj = mf + ifelse(stricths,1,`mfcorrhs',`mfcorrr')

! best known objective = (-1/2)*sqrt(3) = -0.8660254037844386
! aux1 = (2/13)*13^(1/2)
! aux2 = (3/13)*13^(1/2)
! aux3 = (1/13)*13^(1/2) - (3/26)*13^(1/2)*3^(1/2)
! aux4 = (3/26)*13^(1/2) + (1/13)*13^(1/2)*3^(1/2)
! begin of best known solution belonging to the revised case
! x = aux1 =  0.5547001962252291
! x = aux2 =  0.8320502943378437
! x = aux3 = -0.4432265940102775
! x = aux4 =  0.8964096085841832
! x = aux1 =  0.5547001962252291
! x = aux2 =  0.8320502943378437
! x = aux3 = -0.4432265940102775
! x = aux4 =  0.8964096085841832
! x = 0
! end of best known solution belonging to the revised case
End Equations
End Model
```

Stephan K.H. Seidl