diff --git a/Case1_ParticleSettling/DeterministicCalibration/PSUADE/1_ResponseSurface/Quadratic/PSbuildRS_2.in b/Case1_ParticleSettling/DeterministicCalibration/PSUADE/1_ResponseSurface/Quadratic/PSbuildRS_2.in
new file mode 100644
index 0000000000000000000000000000000000000000..4d34a1e700230cb485e3219b7ab99d943677be98
--- /dev/null
+++ b/Case1_ParticleSettling/DeterministicCalibration/PSUADE/1_ResponseSurface/Quadratic/PSbuildRS_2.in
@@ -0,0 +1,33 @@
+load ../../psData_OLH_n120_i6_o1_y2 # load data file
+ana_expert      # psuade> analysis expert mode on
+printlevel 3    # psuade> printlevel set from 0 to 3
+rscheck         # Check for best data-fitted surrogate model
+2               # Select Linear regression as response surface method
+1               # Quantity of interest no
+y               # Perform cross-cross_validation
+120             # Enter the number of groups to validate : (2 - 120)
+n               # Random selection of leave-out groups ? (y or n)
+sys mv RSFA_CV_err.m R2_QUAD_120F_RSFA_CV_err.m # Rename the CV error file before next one
+sys mv RSFA_training_err.m R2_QUAD_RSFA_training_err.m # Rename the training error file before next one
+rs2
+256             # Grid resolution ? (32 - 256)
+2               # Available response surface tools:  select Quadratic 
+6               # Enter the input for x axis (1 - 7) :
+5               # Enter the input for y axis (1 - 7), not 1 : 2
+y               # Set other nominal values automatically ? (y or n)
+1               # Enter output number (1 - 1) :
+n               #  Ymin and Ymax found = 2.037723e-01 3.421112e-01. Set lower threshold ? (y or n) :
+n               # Set upper threshold ? (y or n) : n
+sys mv matlabrs2.m R2_QUAD_rs2_x6x5.m    # Rename the response surface constructed
+rs2
+256             # Grid resolution ? (32 - 256)
+2               # Available response surface tools:  select LINEAR
+6               # Enter the input for x axis (1 - 7) :
+3               # Enter the input for y axis (1 - 7), not 1 : 2
+y               # Set other nominal values automatically ? (y or n)
+1               # Enter output number (1 - 3) :
+n               #  Ymin and Ymax found = 2.037723e-01 3.421112e-01. Set lower threshold ? (y or n) :
+n               # Set upper threshold ? (y or n) : n
+sys mv matlabrs2.m R2_QUAD_rs2_x6x3.m    # Rename the response surface constructed
+sys ls -alrt R2_QUAD*  # Check most one created after mv
+q