You are an expert in writing code in many different languages. Your goal is to perform code completion for the following code, keeping in mind the rest of the code and the file meta data.
    Metadata details:
    filename: runSim.m
    path of file: ./repos/AUVForwardDynamics
    
    the code of the file until where you have to start completion: 
  function runSim()

% Input function for specifying which maneuver to perform

% u =1,1.5,2m/s
% CASE1: TURNING DEL_R_ORDERED=10DEGREE
% CASE2: TURNING DEL_R_ORDERED=15DEGREE
% CASE3: HORIZONTAL ZIG-ZAG MOTION: DEL_R_ORDERED=10DEGREE INITIALLY
% CASE4: HORIZONTAL ZIG-ZAG MOTION: DEL_R_ORDERED=-10DEGREE INTIALLY
% CASE5: HORIZONTAL ZIG-ZAG MOTION: DEL_R_ORDERED=20DEGREE INITIALLY
% CASE6: HORIZONTAL ZIG-ZAG MOTION: DEL_R_ORDERED=-20DEGREE INITIALLY
% CASE7: VERTICAL   ZIG-ZAG MOTION: DEL_ST_ORDERED= 10DEGREE INITIALLY
% CASE8: VERTICAL   ZIG-ZAG MOTION: DEL_ST_ORDERED=-10DEGREE INITIALLY
% CASE9: STEERING PD CONTROLLER
% CASE10:  DIVING PD CONTROLLER
% CASE 11: TURNING DEL_R_ORDERED = -10 DEGREE, 
% CASE 12: TURNINGDEL_R_ORDERED = -15 DEGREE, 
% CASE 13: VERTICAL ZIG-ZAG MOTION: DEL_ST_ORDERED = 20DEGREE INITIALLY,
% CASE 14: VERTICAL ZIG-ZAG MOTION: DEL_ST_ORDERED = -20 DEGREE INITIALLY
% CASE 16: Roll heel angle simulation. Adjust roll_heel in forwarddyanics2.m
% CASE 16: Pitch heel angle simulation. Adjust pitch_heel in forwarddyanics2.m

% INPUT : An integer specifying case number to be run and timespan 
%         
% OUTPUT: Trajectory of bot
% 
% function runsim() 


  clc;
  close all;
  addpath('actuator dynamics');
  addpath('utils');  
  addpath('PDcontrol');
  addpath('Solvers'
    
function runSim()

% Input function for specifying which maneuver to perform

% u =1,1.5,2m/s
% CASE1: TURNING DEL_R_ORDERED=10DEGREE
% CASE2: TURNING DEL_R_ORDERED=15DEGREE
% CASE3: HORIZONTAL ZIG-ZAG MOTION: DEL_R_ORDERED=10DEGREE INITIALLY
% CASE4: HORIZONTAL ZIG-ZAG MOTION: DEL_R_ORDERED=-10DEGREE INTIALLY
% CASE5: HORIZONTAL ZIG-ZAG MOTION: DEL_R_ORDERED=20DEGREE INITIALLY
% CASE6: HORIZONTAL ZIG-ZAG MOTION: DEL_R_ORDERED=-20DEGREE INITIALLY
% CASE7: VERTICAL   ZIG-ZAG MOTION: DEL_ST_ORDERED= 10DEGREE INITIALLY
% CASE8: VERTICAL   ZIG-ZAG MOTION: DEL_ST_ORDERED=-10DEGREE INITIALLY
% CASE9: STEERING PD CONTROLLER
% CASE10:  DIVING PD CONTROLLER
% CASE 11: TURNING DEL_R_ORDERED = -10 DEGREE, 
% CASE 12: TURNINGDEL_R_ORDERED = -15 DEGREE, 
% CASE 13: VERTICAL ZIG-ZAG MOTION: DEL_ST_ORDERED = 20DEGREE INITIALLY,
% CASE 14: VERTICAL ZIG-ZAG MOTION: DEL_ST_ORDERED = -20 DEGREE INITIALLY
% CASE 16: Roll heel angle simulation. Adjust roll_heel in forwarddyanics2.m
% CASE 16: Pitch heel angle simulation. Adjust pitch_heel in forwarddyanics2.m

% INPUT : An integer specifying case number to be run and timespan 
%         
% OUTPUT: Trajectory of bot
% 
% function runsim() 


  clc;
  close all;
  addpath('actuator dynamics');
  addpath('utils');  
  addpath('PDcontrol');
  addpath('Solvers');
  addpath('excel data');
  addpath('excel data/u_1m_s');
  addpath('excel data/u_1dot5m_s');
  addpath('excel data/u_2m_s');
  addpath(genpath('datapoints'));
  addpath(genpath('sensorModels'));
  addpath(genpath('helperFunctions'));
  addpath(genpath('stateEstimation'));
    
  
  r2d =180/pi;
  lbl = {'t(s)', 'u(m/s)','v(m/s)','w(m/s)',...
         'p(deg/s)','q(deg/s)','r(deg/s)'...
          'X(m)','Y(m)','Z(m)','phi(deg)','theta(deg)','psi(deg)'...
          'Del_s(deg)','Del_bp(deg)','Del_bs(deg)','Del_r(deg)',...
          'Del_delb',...
          'n(RPS)','ax(m/s2)','ay(m/s2)','az(m/s2)',...
          'alphax(deg/s2)','alphay(deg/s2)','alphaz(deg/s2)'};
          
  

%% ************* Code for taking input from user**********************
 
%Initialise state vector X here
global tinc;
tinc = 0.1;
X = init_state();

AUVsensors;
EKFparams;


prompt = 'input case no(input 0 to quit):';
temp = input(prompt);
  
  
if(temp~=0)
     
     caseNo=num2str(temp,'%0d') ;
     prompt = 'input simulation time(sec):';
     end