Parsing NASA thermodynamic data

Objective: CODE: R=8.314 f=fopen(\'THERMO.dat\',\'r\'); Header=fgetl(f); global_temp_range=fgetl(f); global_temp_range=str2num(global_temp_range); line1=fgetl(f); line2=fgetl(f); line3=fgetl(f); for i=1:53 line0=fgetl(f); chemical_element=strtok(line0); a=strfind(line0,\'G \'); local_temp=line0(a+1:end-1); local_temp=str2num(local_temp);…