function [output_args] = readData (input_args)

% Read in training data
fidX = fopen('data/train-images-idx3-ubyte');
fread(fidX, 4, '*int32'); % skip over headers
X = (fread(fidX, [28*28,inf], 'uint8'))'; % Load data into matrix
fclose(fidX);
% Read in training labels
fidLabels = fopen('data/train-labels-idx1-ubyte');
fread(fidLabels, 2, '*int32'); % skip over headers
Labels = fread(fidLabels, 'uint8'); 
fclose(fidLabels);

% Verify it has the correct dimensions
sizeX = size(X)
lengthLabels = length(Labels)


% Read in test data
fidX = fopen('data/t10k-images-idx3-ubyte');
fread(fidX, 4, '*int32'); % skip over headers
testX = (fread(fidX, [28*28,inf], 'uint8'))';
fclose(fidX);
% Read in test labels
fidLabels = fopen('data/t10k-labels-idx1-ubyte');
fread(fidLabels, 2, '*int32'); % skip over headers
testLabels = fread(fidLabels, 'uint8');
fclose(fidLabels);

% Verify it has the correct dimensions
sizeTestX = size(testX)
lengthLabelsTest = length(testLabels)