Objects q = periodic End Objects Connections s = q.x End Connections Constants eps = 0.7 End Constants Parameters d p End Parameters Variables s >= 0 , = 100 stored recovery vx >= 0 # slack variable for energy storage representation vy >= 0 # slack variable for energy recovery representation End Variables Equations minimize p p + recovery/eps - stored >= d p - d = vx- vy stored = p-d + vy recovery = d- p + vx $s = stored - recovery/ eps stored * recovery <= 0 End Equations File *.plt New Trend p s d End File