function [Gamma, phi, Length, Force, eps, Maxiters] = getParams
  % Given data:
  Gamma = 18;     % kN/m^3, unit weight of soil
  phi   = pi/6;   % 30deg,  measure of soil friction
  Length     = 3;      % m,      height of pile above ground cut
  Force     = 30;     % kN/m,   force that soil exerts on pile, per unit m
  eps   = 0.01;
  Maxiters = 1000;
