Project 1 - Parsing NASA thermodynamic data

Command window = 


low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 2.56942078E+00-8.59741137E-05 4.19484589E-08-1.00177799E-11 1.22833691E-15    2'


o =

    '-6.12806624E-09 2.11265971E-12 2.91222592E+04 2.05193346E+00                   4'


low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 3.28253784E+00 1.48308754E-03-7.57966669E-07 2.09470555E-10-2.16717794E-14    2'


o =

    '-9.68129509E-09 3.24372837E-12-1.06394356E+03 3.65767573E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 2.50000001E+00-2.30842973E-11 1.61561948E-14-4.73515235E-18 4.98197357E-22    2'


o =

    ' 2.30081632E-18-9.27732332E-22 2.54736599E+04-4.46682853E-01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 3.33727920E+00-4.94024731E-05 4.99456778E-07-1.79566394E-10 2.00255376E-14    2'


o =

    ' 2.01572094E-08-7.37611761E-12-9.17935173E+02 6.83010238E-01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 3.09288767E+00 5.48429716E-04 1.26505228E-07-8.79461556E-11 1.17412376E-14    2'


o =

    '-3.88113333E-09 1.36411470E-12 3.61508056E+03-1.03925458E-01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 3.03399249E+00 2.17691804E-03-1.64072518E-07-9.70419870E-11 1.68200992E-14    2'


o =

    '-5.48797062E-09 1.77197817E-12-3.02937267E+04-8.49032208E-01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 4.01721090E+00 2.23982013E-03-6.33658150E-07 1.14246370E-10-1.07908535E-14    2'


o =

    '-2.42763894E-08 9.29225124E-12 2.94808040E+02 3.71666245E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 4.16500285E+00 4.90831694E-03-1.90139225E-06 3.71185986E-10-2.87908305E-14    2'


o =

    '-2.15770813E-08 8.62454363E-12-1.77025821E+04 3.43505074E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 2.49266888E+00 4.79889284E-05-7.24335020E-08 3.74291029E-11-4.87277893E-15    2'


o =

    '-7.32234889E-10 2.66521446E-13 8.54438832E+04 4.53130848E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 2.87846473E+00 9.70913681E-04 1.44445655E-07-1.30687849E-10 1.76079383E-14    2'


o =

    ' 3.16217327E-09-1.40609067E-12 7.07972934E+04 2.08401108E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 2.87410113E+00 3.65639292E-03-1.40894597E-06 2.60179549E-10-1.87727567E-14    2'


o =

    '-3.85091153E-09 1.68741719E-12 4.60040401E+04 1.56253185E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 2.29203842E+00 4.65588637E-03-2.01191947E-06 4.17906000E-10-3.39716365E-14    2'


o =

    '-6.68815981E-09 1.94314737E-12 5.04968163E+04-7.69118967E-01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 2.28571772E+00 7.23990037E-03-2.98714348E-06 5.95684644E-10-4.67154394E-14    2'


o =

    '-6.87117425E-09 2.54385734E-12 1.64449988E+04 1.60456433E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 7.48514950E-02 1.33909467E-02-5.73285809E-06 1.22292535E-09-1.01815230E-13    2'


o =

    '-4.84743026E-08 1.66693956E-11-1.02466476E+04-4.64130376E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 2.71518561E+00 2.06252743E-03-9.98825771E-07 2.30053008E-10-2.03647716E-14    2'


o =

    ' 9.07005884E-10-9.04424499E-13-1.43440860E+04 3.50840928E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 3.85746029E+00 4.41437026E-03-2.21481404E-06 5.23490188E-10-4.72084164E-14    2'


o =

    ' 2.45919022E-09-1.43699548E-13-4.83719697E+04 9.90105222E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 2.77217438E+00 4.95695526E-03-2.48445613E-06 5.89161778E-10-5.33508711E-14    2'


o =

    '-1.33144093E-08 4.33768865E-12 3.83956496E+03 3.39437243E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 1.76069008E+00 9.20000082E-03-4.42258813E-06 1.00641212E-09-8.83855640E-14    2'


o =

    '-3.79285261E-08 1.31772652E-11-1.43089567E+04 6.02812900E-01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 3.69266569E+00 8.64576797E-03-3.75101120E-06 7.87234636E-10-6.48554201E-14    2'


o =

    '-1.04532012E-08 4.36967278E-12-3.19391367E+03 5.47302243E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.00'


mid_temp =

    '1000.000'


high_temp =

    '3000.00'


m =

    ' 0.03770799E+02 0.07871497E-01-0.02656384E-04 0.03944431E-08-0.02112616E-12    2'


o =

    '-0.07377636E-07 0.02075610E-10 0.09786011E+04 0.13152177E+02                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 1.78970791E+00 1.40938292E-02-6.36500835E-06 1.38171085E-09-1.17060220E-13    2'


o =

    '-7.10806889E-08 2.61352698E-11-2.56427656E+04-1.50409823E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 3.16780652E+00 4.75221902E-03-1.83787077E-06 3.04190252E-10-1.77232770E-14    2'


o =

    ' 2.94791045E-08-1.09331511E-11 6.68393932E+04 6.22296438E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 4.14756964E+00 5.96166664E-03-2.37294852E-06 4.67412171E-10-3.61235213E-14    2'


o =

    ' 2.80152437E-08-8.50072974E-12 2.64289807E+04 1.39397051E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 3.01672400E+00 1.03302292E-02-4.68082349E-06 1.01763288E-09-8.62607041E-14    2'


o =

    '-3.57657847E-08 1.47150873E-11 3.48598468E+04 8.51054025E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 2.03611116E+00 1.46454151E-02-6.71077915E-06 1.47222923E-09-1.25706061E-13    2'


o =

    '-6.91588753E-08 2.69884373E-11 5.08977593E+03 4.09733096E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 1.95465642E+00 1.73972722E-02-7.98206668E-06 1.75217689E-09-1.49641576E-13    2'


o =

    '-5.99126606E-08 2.30509004E-11 1.28416265E+04 4.70720924E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 1.07188150E+00 2.16852677E-02-1.00256067E-05 2.21412001E-09-1.90002890E-13    2'


o =

    '-7.08466285E-08 2.68685771E-11-1.15222055E+04 2.66682316E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '3500.000'


m =

    ' 4.51129732E+00 9.00359745E-03-4.16939635E-06 9.23345882E-10-7.94838201E-14    2'


o =

    ' 9.34397568E-09-2.01457615E-12-7.04291804E+03 1.22156480E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.00'


mid_temp =

    '1000.000'


high_temp =

    '4000.00'


m =

    ' 0.56282058E+01 0.40853401E-02-0.15934547E-05 0.28626052E-09-0.19407832E-13    2'


o =

    ' 0.17275759E-07-0.50664811E-11 0.20059449E+05 0.12490417E+02                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.000'


mid_temp =

    '1000.000'


high_temp =

    '5000.000'


m =

    ' 0.59238291E+01 0.67923600E-02-0.25658564E-05 0.44987841E-09-0.29940101E-13    2'


o =

    ' 0.43137131E-07-0.14014594E-10 0.80316143E+04 0.13874319E+02                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.00'


mid_temp =

    '1000.000'


high_temp =

    '4000.000'


m =

    ' 0.52097030E+01 0.29692911E-02-0.28555891E-06-0.16355500E-09 0.30432589E-13    2'


o =

    '-0.16226120E-08-0.23511081E-12 0.28637820E+05 0.89927511E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.38022392E+01 0.31464228E-02-0.10632185E-05 0.16619757E-09-0.97997570E-14    2'


o =

    ' 0.10092349E-07-0.30089028E-11 0.14712633E+05 0.89164419E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.29792509E+01 0.34944059E-02-0.78549778E-06 0.57479594E-10-0.19335916E-15    2'


o =

    '-0.17137094E-07 0.55454573E-11 0.11548297E+05 0.17498417E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.24159429E+01 0.17489065E-03-0.11902369E-06 0.30226245E-10-0.20360982E-14    2'


o =

    ' 0.00000000E+00 0.00000000E+00 0.56104637E+05 0.41939087E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.37667544E+01 0.28915082E-02-0.10416620E-05 0.16842594E-09-0.10091896E-13    2'


o =

    '-0.21726464E-07 0.79469539E-11 0.28791973E+05 0.29779410E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.48230729E+01 0.26270251E-02-0.95850874E-06 0.16000712E-09-0.97752303E-14    2'


o =

    ' 0.96819806E-08-0.29307182E-11 0.87417744E+04 0.10757992E+02                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.27836928E+01 0.13298430E-02-0.42478047E-06 0.78348501E-10-0.55044470E-14    2'


o =

    ' 0.24816442E-08-0.10356967E-11 0.41880629E+05 0.18483278E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.28347421E+01 0.32073082E-02-0.93390804E-06 0.13702953E-09-0.79206144E-14    2'


o =

    '-0.56115197E-08 0.16440717E-11 0.21885910E+05-0.14184248E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.26344521E+01 0.56662560E-02-0.17278676E-05 0.23867161E-09-0.12578786E-13    2'


o =

    '-0.22808887E-07 0.82638046E-11-0.67417285E+04-0.62537277E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.32606056E+01 0.11911043E-02-0.42917048E-06 0.69457669E-10-0.40336099E-14    2'


o =

    '-0.93361354E-08 0.28035770E-11 0.98446230E+04 0.22808464E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.48847542E+01 0.21723956E-02-0.82806906E-06 0.15747510E-09-0.10510895E-13    2'


o =

    '-0.20475426E-07 0.78350564E-11 0.28966179E+04 0.63119917E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.000'


mid_temp =

    '1382.000'


high_temp =

    '5000.000'


m =

    ' 6.59860456E+00 3.02778626E-03-1.07704346E-06 1.71666528E-10-1.01439391E-14    2'


o =

    ' 4.41432836E-09-7.57521466E-13 1.92990252E+04 1.07332972E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.000'


mid_temp =

    '1368.000'


high_temp =

    '5000.000'


m =

    ' 5.89784885E+00 3.16789393E-03-1.11801064E-06 1.77243144E-10-1.04339177E-14    2'


o =

    ' 5.17195767E-10 1.19360788E-14-2.82698400E+03 5.63292162E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.000'


mid_temp =

    '1478.000'


high_temp =

    '5000.000'


m =

    ' 6.22395134E+00 3.17864004E-03-1.09378755E-06 1.70735163E-10-9.95021955E-15    2'


o =

    '-6.61271298E-10 3.62235752E-13-1.55873636E+04 6.19457727E+00                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.51521845E+01 0.23051761E-02-0.88033153E-06 0.14789098E-09-0.90977996E-14    2'


o =

    ' 0.48016964E-08-0.13313595E-11 0.14682477E+05 0.95504646E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.37459805E+01 0.43450775E-04 0.29705984E-06-0.68651806E-10 0.44134173E-14    2'


o =

    '-0.31516323E-09-0.46430356E-12 0.51708340E+05 0.39804995E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.000'


mid_temp =

    '1000.000'


high_temp =

    '5000.000'


m =

    ' 0.58946362E+01 0.39895959E-02-0.15982380E-05 0.29249395E-09-0.20094686E-13    2'


o =

    ' 0.12125540E-07-0.32357378E-11 0.54261984E+05 0.11675870E+02                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.000'


mid_temp =

    '1000.000'


high_temp =

    '5000.000'


m =

    ' 0.02926640E+02 0.14879768E-02-0.05684760E-05 0.10097038E-09-0.06753351E-13    2'


o =

    ' 0.05641515E-07-0.02444854E-10-0.10208999E+04 0.03950372E+02                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.000'


mid_temp =

    '1000.000'


high_temp =

    '5000.000'


m =

    ' 0.02500000E+02 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00    2'


o =

    ' 0.00000000E+00 0.00000000E+00-0.07453750E+04 0.04366000E+02                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.000'


mid_temp =

    '1000.000'


high_temp =

    '5000.000'


m =

    ' 0.75341368E+01 0.18872239E-01-0.62718491E-05 0.91475649E-09-0.47838069E-13    2'


o =

    '-0.21977499E-07 0.95149253E-11-0.13958520E+05 0.19201691E+02                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.000'


mid_temp =

    '1000.000'


high_temp =

    '5000.000'


m =

    ' 0.77026987E+01 0.16044203E-01-0.52833220E-05 0.76298590E-09-0.39392284E-13    2'


o =

    '-0.19609569E-07 0.93732470E-11 0.10631863E+05 0.21122559E+02                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '200.000'


mid_temp =

    '1000.000'


high_temp =

    '6000.000'


m =

    ' 0.54041108E+01 0.11723059E-01-0.42263137E-05 0.68372451E-09-0.40984863E-13    2'


o =

    '-0.57458611E-07 0.21931112E-10-0.21572878E+05 0.41030159E+01                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 

low_temp =

    '300.000'


mid_temp =

    '1000.000'


high_temp =

    '5000.000'


m =

    ' 0.05975670E+02 0.08130591E-01-0.02743624E-04 0.04070304E-08-0.02176017E-12    2'


o =

    '-0.07158583E-07 0.02867385E-10 0.15214766E+04 0.09558290E+02                   4'

Warning: Directory already exists. 
> In file_writing (line 93) 
>> 

I also put 3 figures related to specific heat, entalpy and entropy.

NASA file parsing is saved in my desktop, I have put screenshots all of them.
(1). nasa3
(2). nasa2
(3). nasa4
(4). nasa5

NASA file parsing = 

clear all
close all
clc

%getting data from prvided file
f1 = fopen('THERMO.dat','r');

for i = 1:5
    fgetl(f1);
end
for i = 1:53


    a = fgetl(f1);
    b = strsplit(a, ' ');
    species = (b{1});


%low_temp = (b{6});
%low_temp = str2num(low_temp);

%mid_temp = (b{7});
%mid_temp = str2num(mid_temp);

%high_temp = (b{8});
%high_temp = str2num(high_temp);


low_temp = b{end-3}
low_temp= str2num(low_temp);

mid_temp = b{end-1}
mid_temp = str2num(mid_temp);

high_temp = b{end-2}
high_temp = str2num(high_temp);


m = fgetl(f1)
A = strfind(m, 'E');



% To calculate coefficients
t1 = m(1:A(1)+3);
t1 = str2num(t1);
t2 = m(A(1)+4:A(2)+3);
t2 = str2num(t2);
t3 = m(A(2)+4:A(3)+3);
t3 = str2num(t3);
t4 = m(A(3)+4:A(4)+3);
t4 = str2num(t4);
t5 = m(A(4)+4:A(5)+3);
t5 = str2num(t5);
n = fgetl(f1);
B = strfind(n, 'E');
t6 = n(1:B(1)+3);
t6 = str2num(t6);
t7 = n(B(1)+4:B(2)+3);
t7 = str2num(t7);
t8 = n(B(2)+4:B(3)+3);
t8 = str2num(t8);
t9 = n(B(3)+4:B(4)+3);
t9 = str2num(t9);
t10 = n(B(4)+4:B(5)+3);
t10 = str2num(t10);
o = fgetl(f1)
C = strfind(o, 'E');
t11 = o(1:C(1)+3);
t11 = str2num(t11);
t12 = o(C(1)+4:C(2)+3);
t12 = str2num(t12);
t13 = o(C(2)+4:C(3)+3);
t13 = str2num(t13);
t14 = o(C(3)+4:C(4)+3);
t14 = str2num(t14);

T = linspace(low_temp, high_temp,500);
R = 8.314;

% specific heat, Enthalpy and Entropy calculation
for i = 1:length(T)
    if T(i)>mid_temp
        cp(i) = (t1+t2*T(i)+t3*T(i).^2+t4*T(i).^3+t5*T(i).^4)*R;
        H(i) = R*T(i).*(t1+((t2*T(i)/2))+((t3*T(i).^2)/3)+((t4*T(i).^3)/4)+((t5*T(i).^4)/5)+((t6./T(i))));
        S(i) = (t1.*log(T(i))+t2.*T(i)+t3.*(T(i)^2)/2+t4.*(T(i)^3)/3+t5.*(T(i)^4)/4+t7)*R;
    else
        cp(i) = (t8+t9.*T(i)+t10.*T(i).^2+t11.*T(i).^3+t12.*T(i).^4)*R;
        H(i) = R*T(i).*(t8+((t9*T(i)/2))+((t10*T(i).^2)/3)+((t11*T(i).^3)/4)+((t12*T(i).^4)/5)+((t13./T(i))));
        S(i) = (t8.*log(T(i))+t9.*T(i)+t10.*(T(i)^2)/2+t11.*(T(i)^3)/3+t12.*(T(i)^4)/4+t14)*R;
    end
end
mkdir(species)                                                                          %mkdir = making directory where the
cd(species)                                                                                      %the program save
% To plot variation of cp with T
figure(1)
plot(T, cp, 'linewidth',2,'Color','r')
xlabel('Temperature in K')
ylabel('specific heat in KJ/kg.K')                                                      %cp = specific heat
grid on
title(sprintf('variation of specific heat with temperature for %s',species))            %s  = string
saveas(figure(1),'cp.png');

% To plot variation of enthalpy with temperature
figure(2)
plot(T, H, 'linewidth',2,'Color','g')                                                   %H = enthalpy
xlabel('Temperature in K')
ylabel('Enthalpy in KJ/kg')
grid on
title(sprintf('variation of Enthalpy with temperature for %s',species))
saveas(figure(2),'H.png');


% To plot variation of Entropy with temperature
figure(3)
plot(T, S, 'linewidth',2,'Color','b')
xlabel('Temperature in K')
ylabel('Entropy in KJ/K')                                                                % S = entropy
grid on
title(sprintf('variation of Entropy with temperature for %s',species))
saveas(figure(3),'S.png');
cd ..
end