clear; close all; clc; syms pe; fprintf('***入力パラメータ***\n'); Pc = input('燃焼室圧力 Pc [Mpa]: '); e = input('ノズル開口比 e=Ae/At: '); k = input('比熱比 k: '); pb = input('背圧 pb [MPa] (デフォルト=大気圧: 0.1013MPa): '); if isempty(…