--- a
+++ b/Application/Examples/BikeSpring/Model/parameterstudy.m
@@ -0,0 +1,47 @@
+% Systematic variation of two parameters
+
+% Design variables as they appear in the AnyScript model
+    % AnyVar L0 = 0.24;             //Spring slack length
+    % AnyVar F1 = 100;              //Spring force at 50% strain
+    % AnyVar SpringRad = 0.15;      //Radius of spring attachment on the crank
+    % AnyVar SpringAngle = 45;      //Angle of spring attachment points on the crank
+    % AnyVar SpringX = -0.10;       //X coordinate of spring attachment on frame
+    % AnyVar SpringY = 0.30
+    
+% Initialize 
+%        L0        F1    SpringRad  SpringAngle  SpringY
+xmin = [0.1       0.0      0.01        0.0         0.1];   % Lower bound
+x =    [0.2897  114.5789    0.1673   49.2547    0.4311];   % Initial guess
+xmax = [0.4    1000.0      0.40      180.0         0.5];   % Upper bound
+
+% Variable choice
+n1 = 4;
+n2 = 5;
+nStep = 8;
+
+% Initialize temporary variable xx
+d1 = (xmax(n1)-xmin(n1))/(nStep-1);
+d2 = (xmax(n2)-xmin(n2))/(nStep-1);
+
+% Fill abscissa arrays
+for i=0:nStep-1
+    p1(i+1) = xmin(n1) + i*d1;
+    p2(i+1) = xmin(n2) + i*d2;
+end
+
+% Fill ordinate array
+count = 1;
+xx = x;
+for i=1:nStep
+    xx(n1) = p1(i);
+    for j=1:nStep
+        xx(n2) = p2(j);
+        result(i,j) = objective(xx);
+        count = count + 1;
+    end
+end
+ 
+%Plotting
+surf(p1,p2,result);
+ 
+clear;