%CME434:  Solution to HW4G by TSD
clear all;
status = fclose('all')		%Make sure all files are closed

global   nn ne nm ndim nep nqn
global   nq nf nmpc
global   XYZ CTAB AE QTAB FTAB MATERIAL
global   FIDI_NAME ns
global   k l m
global   K L F P Q

FIDI_NAME = 'hw40data.txt'	%Name of input data file
%FIDI_NAME = 'hw41data.txt'
%FIDI_NAME = 'hw42data.txt'
%FIDI_NAME = 'hw43data.txt'

hw4g_sub;	%call functions