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: loadGeoTerrain.m
    path of file: ./repos/monster/layout
    
    the code of the file until where you have to start completion: function [terrain, trajectories] = loadGeoTerrain(Monster, numUsers)
	%% About
	% Utility to load terrain information from shapefile files
	% 
	% The function loads a shapefile with roads, then drops a user in a random
	% road. A walk is simulated along the road where the user has
	% been dropped for up to the simulation time

%% Load geo data

seed = Monster.Runtime.seed;
velocity = 1.5; % m/s
time = Monster.Runtime.totalRounds*10e-4; % s
roadsFileName = Monster.Config.Terrain.roadsFile;
buildingsFileName = Monster.Config.Terrain.buildingsFile;
roadsShapeFile = shaperead(roadsFileName);
buildingsShapeFile = shaperead(buildingsFileName);
% Extract road coordinates from the shape file
% convert from geodesic WGS84 to ECEF
% elevation is given as an average elevation of the area
averageElevation = 33; % in m

% Create geodesic ellipsoid
wgs84 = wgs84Ellipsoid();


    
function [terrain, trajectories] = loadGeoTerrain(Monster, numUsers)
	%% About
	% Utility to load terrain information from shapefile files
	% 
	% The function loads a shapefile with roads, then drops a user in a random
	% road. A walk is simulated along the road where the user has
	% been dropped for up to the simulation time

%% Load geo data

seed = Monster.Runtime.seed;
velocity = 1.5; % m/s
time = Monster.Runtime.totalRounds*10e-4; % s
roadsFileName = Monster.Config.Terrain.roadsFile;
buildingsFileName = Monster.Config.Terrain.buildingsFile;
roadsShapeFile = shaperead(roadsFileName);
buildingsShapeFile = shaperead(buildingsFileName);
% Extract road coordinates from the shape file
% convert from geodesic WGS84 to ECEF
% elevation is given as an average elevation of the area
averageElevation = 33; % in m

% Create geodesic ellipsoid
wgs84 = wgs84Ellipsoid();

for iRoad = 1: length(roadsShapeFile)
	geoLat = roadsShapeFile(iRoad).Y;
	geoLon = roadsShapeFile(iRoad).X;
	geoH = averageElevation * ones(1, length(geoLat));
	[x, y, ~] = geodetic2ecef(wgs84, geoLat, geoLon, geoH);
	roads(iRoad) = struct(...,
		'x',x,...
		'y', y);
end