tp107.apm


Model tp107
  ! Source version 1

  Parameters
    c = (48.4/50.176)*sin(.25)
    d = (48.4/50.176)*cos(.25)
  End Parameters

  Variables
    x[1] =  .8,    >= 0
    x[2] =  .8,    >= 0
    x[3] =  .2
    x[4] =  .2
    x[5] = 1.0454, >=  .90909, <= 1.0909
    x[6] = 1.0454, >=  .90909, <= 1.0909
    x[7] = 1.0454, >=  .90909, <= 1.0909
    x[8] = 0
    x[9] = 0
    obj
  End Variables

  Intermediates
    y1 = sin(x[8])
    y2 = cos(x[8])
    y3 = sin(x[9])
    y4 = cos(x[9])
    y5 = sin(x[8] - x[9])
    y6 = cos(x[8] - x[9])
    ec[1] = .4 - x[1] + 2*c*x[5]^2  &
          - x[5]*x[6]*(d*y1 + c*y2) &
          - x[5]*x[7]*(d*y3 + c*y4)
    ec[2] = .4 - x[2] + 2*c*x[6]^2  &
          + x[5]*x[6]*(d*y1 - c*y2) &
          + x[6]*x[7]*(d*y5 - c*y6)
    ec[3] = .8 + 2*c*x[7]^2         &
          + x[5]*x[7]*(d*y3 - c*y4) &
          - x[6]*x[7]*(d*y5 + c*y6)
    ec[4] = .2 - x[3] + 2*d*x[5]^2  &
          + x[5]*x[6]*(c*y1 - d*y2) &
          + x[5]*x[7]*(c*y3 - d*y4)
    ec[5] = .2 - x[4] + 2*d*x[6]^2  &
          - x[5]*x[6]*(c*y1 + d*y2) &
          - x[6]*x[7]*(c*y5 + d*y6)
    ec[6] = (-.337) + 2*d*x[7]^2    &
          - x[5]*x[7]*(c*y3 + d*y4) &
          + x[6]*x[7]*(c*y5 - d*y6)
    mf = 3000*x[1] + 1000*x[1]^3    &
       + 2000*x[2] + 666.667*x[2]^3
  End Intermediates

  Equations
    ec[1:6] = 0

    obj = mf

    ! best known objective = 5055.011803543816
    ! begin of best known solution
    ! x[1] =  0.6670127623153279
    ! x[2] =  1.022384717050832
    ! x[3] =  0.2282871446675199
    ! x[4] =  0.1848217561297505
    ! x[5] =  1.0909
    ! x[6] =  1.0909
    ! x[7] =  1.069036003142822
    ! x[8] =  0.1066106224065641
    ! x[9] = -0.3387876190699065
    ! end of best known solution
  End Equations
End Model

Stephan K.H. Seidl