def make_full_report(list_output_results): list_ultimate_points = [] list_yield_points = [] list_stiffness = [] # Creating subsequent force displacement curves for i in range(0, len(list_output_results)): plots = Plotter(list_output_results[i], "plot" + str(i)) plots.plot_xy_data(r"Rotation $\theta$ ($deg$)", r"Torque ($Nm$)") # storing the ultimates (maximum) and offset yields list_ultimate_points.append(list_output_results[i].ultimate_point.y) list_yield_points.append(list_output_results[i].yield_point.y) list_stiffness.append(list_output_results[i].stiffness) mean_ultimate = np.mean(np.asarray(list_ultimate_points)) mean_yield = np.mean(np.asarray(list_yield_points)) std_ultimate = np.std(np.asarray(list_ultimate_points)) std_yield = np.std(np.asarray(list_yield_points)) print("mean of ultimates = ", mean_ultimate) print("mean of yields = ", mean_yield) return "done"