Nastran Load case combination

Hello,

I am setting up a FEM model to perform optimization. We have a total of 604 load cases which are combination of individual loads. Insteading of modifying manually in the input file or graphically creating all the load cases which is laboric I would like to create a matlab code to generate combination of load cases. The nastran format is shown as below, where 10 is the load case no, 2,4,7,8 are the individual loadcases with 1.000000 as the weight of each load: LOAD 101.0000001.000000 21.000000 41.000000 7 + 1.000000 8 Can anyone send me if somebody has already written code for this purpose (or) help me in writing the code.

Thank you