Three types of weight in both energy saving and optimum room temperature environment were changed to obtain the optimal control solution of the ground thermal storage system in a greenhouse. The relation diagram between the optimal solution of a performance function, and the state constraints and control function constraints was created in consideration of the energy term in a control function value area and the room temperature environment. As a result, the whole image of the performance function could be grasped in consideration of the energy term with inequality constraints and the room temperature environmental term in this study. The rate of a weighting factor in the performance function significantly influences the optimal solution. The influence on the optimal solution also changes when the optimal room temperature schedule differs. The influence that three types of rising algorithm exert on the convergence and converging speed was investigated. Superiority or inferiority occurs according to the space properties of a performance function. A zigzag method is most disadvantageous. The constraints can be converged to the optimal solution using an SUMT outer point method irrespective of the initial value. 6 refs., 6 figs., 4 tabs.