--- a
+++ b/Application/Examples/BikeSpring/Model/ObjectivePenal.m
@@ -0,0 +1,15 @@
+% Computes the pseudo-objective function penalized for side constraint
+% violations
+
+function [Fps] = ObjectivePenal(xmin,x,xmax)
+
+penalty = 2;  % exterior penalty factor
+Fps = objective(x);
+for j=1:length(x)
+    if xmin(j) > x(j)
+        Fps = Fps*(1 + penalty*(xmin(j)-x(j))^2);
+    end
+    if x(j) > xmax(j)
+        Fps = Fps*(1 + penalty*(xmin(j)-x(j))^2);
+    end
+end