Update for 1.6

2025-04-28 19:48:28 +02:00
1 parent 3b2c67d5ed
commit 41ae21469e
@@ -382,6 +382,7 @@ damaging-speed = 250
 blowout-speed = 300
 maneuvering_flaps = 0
 lift_scalar = 0.5
+// KH From 0.75
 drag_scalar = 0.475
 pitch_scalar = 1
 max_on_ground_position = 1
@@ -397,6 +398,7 @@ damaging-speed = 250
 blowout-speed = 300
 maneuvering_flaps = 0
 lift_scalar = 0.5
+// KH From 0.75
 drag_scalar = 0.475
 pitch_scalar = 1
 max_on_ground_position = 1
@@ -1,8 +1,7 @@
-//TFDi Designs MD-11 (GE CF6-80C2D1F) vRC1.6.4_GE_F 10DEC24
+//TFDi Designs MD-11 (GE CF6-80C2D1F) vRC1.6.5_GE 04APR25

 [VERSION]
-//KH Bumped major
-major = 4
+major = 3
 minor = 1

 [GENERALENGINEDATA]
@@ -11,7 +10,7 @@ fuel_flow_scalar =0.031
 min_throttle_limit =-0.25
 master_ignition_switch =0
 starter_type =0
-max_contrail_temperature =-30
+max_contrail_temperature =-25
 Engine.0 = -3.00, -27.75, -8.50
 ThrustAnglesPitchHeading.0 = 0, 0
 Engine.1 = -103.50, 0.00, 15.00
@@ -29,8 +28,7 @@ DisablePropellerControls =1
 DisableAutopilotControls=1

 [TURBINEENGINEDATA]
-//KH From 1, as per SDK warning
-fuel_flow_gain=0
+fuel_flow_gain=1
 low_idle_n1=24.2
 low_idle_n2=62.9
 inlet_area=9.92
@@ -40,19 +38,18 @@ reverser_available=1
 reverser_mach_controlled=0
 afterburner_available=0
 afterburner_throttle_threshold=0
-//KH Unchanged, 0.304 to .344 according to GE, leaning towards .344 for the 80C2D1F
 ThrustSpecificFuelConsumption=0.362
 use_commanded_Ne_Table=1
 use_corrected_N2_from_FF_table=1
 use_n2_to_n1_table=1
 use_old_fuelflow_simvar=1
-//KH START Modified
 n2_to_n1_table=0:0:0.85,0:0:0,10:1.1:1,20:2.7:2.4,30:5.7:4.7,40:11.1:9.78,50:16.35:14.9,64.5:25.34:25.39,70:38:37.16,80:59.2:58.56,90:80.4:79.96,99.6:110.75:110.5,105:112:112
-mach_0_corrected_commanded_ne_table=0:1:5.42,0:64:64.5,0.143:69.214:69.857,0.286:74.428:75.214,0.429:79.642:80.571,0.572:84.856:85.928,0.715:90.07:91.285,0.858:95.284:96.642,1:102.5:104
-mach_hi_corrected_commanded_ne_table=0.88:1:5.42,0:66:67.5,0.143:71.143:72.571,0.286:76.286:77.642,0.429:81.429:82.713,0.572:86.572:87.784,0.715:91.715:92.855,0.858:96.858:97.926,1:108:109
-corrected_n2_from_ff_table=0:17,0.00991:64.5,0.01992:70,0.02993:75,0.03994:80,0.04996:85,0.05997:90,0.06998:95,0.07999:100,0.09:106.5
-n1_and_mach_on_thrust_table=0:0:0.25:0.39:0.49:0.5:0.67:0.74:0.8:0.82:0.85,0:0:0:0:0:0:0:0:0:0:0,20:0.0304:0.0598:0.0598:0.0598:0.0205:0.0282:0.0404:0.0411:0.0224:0.0297,25:0.0408:0.0704:0.0704:0.1304:0.1291:0.1235:0.1278:0.1262:0.1291:0.0393,30:0.0635:0.0821:0.0821:0.0335:0.0547:0.0847:0.0847:0.1039:0.1039:0.1039,35:0.0924:0.0963:0.0963:0.0425:0.0601:0.0901:0.0901:0.1053:0.1053:0.1053,40:0.1212:0.1134:0.1134:0.1007:0.092:0.102:0.102:0.1156:0.1161:0.1176,45:0.147:0.1278:0.1278:0.13:0.3795:0.4095:0.4095:0.4235:0.4235:0.4235,50:0.1816:0.1448:0.1448:0.1518:0.3853:0.4153:0.4153:0.4389:0.4389:0.4389,55:0.2388:0.1711:0.1632:0.1715:0.3906:0.4206:0.4206:0.442:0.442:0.442,60:0.2726:0.2079:0.1902:0.2076:0.3943:0.4243:0.4243:0.446:0.446:0.446,65:0.3294:0.2532:0.2336:0.2367:0.4002:0.4302:0.4302:0.4542:0.4542:0.4542,70:0.3807:0.3116:0.2786:0.2768:0.409:0.439:0.439:0.4642:0.4642:0.4642,75:0.4477:0.3938:0.3318:0.3282:0.4274:0.4574:0.4574:0.4705:0.4703:0.4706,80:0.5395:0.4915:0.4124:0.4138:0.4321:0.4621:0.4621:0.484:0.4813:0.4805,85:0.6524:0.6349:0.5626:0.5072:0.5044:0.5355:0.5309:0.5785:0.5275:0.5356,90:0.7545:0.7735:0.6998:0.61:0.61:0.6173:0.6574:0.6638:0.6044:0.6297,95:0.8364:0.8395:0.708:0.69:0.701:0.7296:0.8207:0.8673:0.7222:0.7682,103.5:0.9779:0.9811:0.7157:0.7212:0.744:0.7453:0.8651:0.9026:0.9209:0.9312,110:0.9853:0.9826:0.7775:0.83:0.85:0.8257:0.8749:0.9122:0.9322:0.9317,120:1:0.9945:0.7877:0.84:0.88:0.8339:0.8818:0.9131:0.9473:0.9404
-//KH END Modified
+mach_0_corrected_commanded_ne_table=0:1:5.42,0:64:64.5,0.143:69.214:69.857,0.286:74.428:75.214,0.429:79.642:80.571,0.572:84.856:85.928,0.715:90.07:91.285,0.858:95.284:96.642,1:100.5:102
+mach_hi_corrected_commanded_ne_table=0.88:1:5.42,0:66:67.5,0.143:71.143:72.571,0.286:76.286:77.642,0.429:81.429:82.713,0.572:86.572:87.784,0.715:91.715:92.855,0.858:96.858:97.926,1:102:103
+corrected_n2_from_ff_table=0:17,0.0099:63.7,0.02:70.5,0.03:73.1,0.04:79.4,0.05:83.6,0.06:89.5,0.07:94.8,0.08:100.1,0.09:106.5
+//KH START Modified
+n1_and_mach_on_thrust_table=0.000000:0.000000:0.206803:0.280661:0.354519:0.428377:0.502235:0.576093:0.649951:0.723810:0.797668,0.000000:0.000000:0.000000:0.000000:0.000000:0.000000:0.000000:0.000000:0.000000:0.000000:0.000000,20.000000:0.020334:0.018926:0.018365:0.017768:0.017149:0.016515:0.015863:0.015208:0.014534:0.013860,25.000000:0.029434:0.027422:0.026632:0.025797:0.024940:0.024059:0.023152:0.022235:0.021306:0.020375,30.000000:0.033986:0.030672:0.028919:0.026653:0.025928:0.025184:0.025802:0.021374:0.018702:0.023291,35.000000:0.070581:0.064962:0.062424:0.059648:0.056670:0.053536:0.050308:0.046924:0.043444:0.039910,40.000000:0.115706:0.107129:0.103436:0.099453:0.095210:0.090734:0.086072:0.081361:0.076653:0.071841,45.000000:0.164968:0.152931:0.147828:0.142341:0.136513:0.130440:0.124498:0.118283:0.111870:0.105388,50.000000:0.218546:0.202746:0.196107:0.188982:0.181423:0.174048:0.166338:0.158287:0.149967:0.142050,55.000000:0.276676:0.256828:0.248552:0.239686:0.230419:0.221263:0.211681:0.201686:0.191443:0.181712,60.000000:0.339274:0.315126:0.305110:0.294428:0.283152:0.272084:0.260539:0.248522:0.236179:0.224430,65.000000:0.406072:0.377362:0.365531:0.352959:0.339697:0.326347:0.312745:0.298645:0.284078:0.270104,70.000000:0.476658:0.443151:0.429436:0.414899:0.399618:0.383785:0.368032:0.351772:0.335030:0.318542,75.000000:0.550568:0.512033:0.496350:0.479778:0.462414:0.444307:0.426075:0.407561:0.388596:0.369500,80.000000:0.627228:0.583453:0.565725:0.547049:0.527537:0.507273:0.486477:0.465611:0.444349:0.422704,85.000000:0.705886:0.656705:0.636874:0.616041:0.594344:0.571897:0.548783:0.525433:0.501810:0.477902,90.000000:0.785094:0.730712:0.708789:0.685821:0.661969:0.637381:0.612132:0.586376:0.560384:0.534209,95.000000:0.863065:0.803680:0.779839:0.754920:0.729110:0.702557:0.675366:0.647478:0.619241:0.590907,100.000000:0.936071:0.872277:0.847042:0.820755:0.793551:0.765559:0.736891:0.707114:0.677066:0.646923,105.000000:0.994083:0.928996:0.903986:0.877967:0.850934:0.822172:0.792494:0.762245:0.731543:0.700495,110.000000:0.974676:0.921881:0.905331:0.888737:0.871094:0.851457:0.828909:0.802711:0.772336:0.737555
+//KH START Modified
 corrected_airflow_table=0:0:0.9,0:0:0,20:4.29:14.08,85:4.62:24,90:5.75:25,95:7.23:27,100:8.36:28,105:9.29:29,120:10.94:30
 JET_density_on_FF_table=0.000462:28.74,0.000587:32.93,0.000738:35.02,0.000891:29.46,0.001066:31.02,0.001267:32.56,0.001496:34.54,0.0017565:35.11,0.002048:30.8,0.002377:22.55
 epr_max=1.8
@@ -66,16 +63,12 @@ oil_temp_factor_from_rpm=0.00:0.80, 250.00:0.80, 500.00:0.80, 2500.00:1.00, 4000
 oil_press_max=13000
 oil_press_tc=0.8
 oil_press_tuning_constant=1
-//KH From 2141, 600 above EGT
-itt_peak_temperature=2531
-itt_tc=3
+itt_peak_temperature=2141
+itt_tc=2
 itt_tuning_constant=1
-//KH From 1600, 800C
-egt_peak_temperature=1931
-//KH From 0.1, taken from the A300, made sense
-egt_tc=3
-//KH From 0.1
-egt_tuning_constant=1
+egt_peak_temperature=1600
+egt_tc=0.1
+egt_tuning_constant=0.1
 fuel_press_max=2376
 fuel_press_tc=2
 fuel_press_tuning_constant=1
@@ -94,15 +87,11 @@ min_n2_for_apu_bleed_air_cutoff=50
 n1_normal_tc=0.5
 n1_start_max_rate=1
 n1_start_combustion_max_rate=2
-//KH Unchanged, IRL closer to 1200
 idle_fuel_flow=600
 idle_high_fuel_flow=40
 min_n2_for_fuel_flow=17
-//KH Unchanged, IRL 117.5 as per EASA TCDS for CF6-80C2D1F
 high_n1=110
-//KH Unchanged, IRL 112.5 as per EASA TCDS for CF6-80C2D1F
 high_n2=105
-//KH Unchanged, IRL closer to 15500
 high_fuel_flow=9800
 mach_influence_on_n1=10
 fuel_flow_max=10000
@@ -112,22 +101,15 @@ fuel_flow_controller_d=0.1
 fuel_flow_controller_iboundary=10
 fuel_flow_controller_dboundary=100
 max_torque_protection=0
-//KH Disabled
-max_n1_protection=0
-//KH Disabled
-max_n2_protection=0
-//KH Disabled
-max_egt_protection=0
-//KH From 2
-fuel_flow_max_itt_factor=2.5
-//KH From 0.9
-fuel_flow_min_itt_factor=0.74
+max_n1_protection=110
+max_n2_protection=100
+max_egt_protection=2000
+fuel_flow_max_itt_factor=2
+fuel_flow_min_itt_factor=0.9
 n1_cooling_factor=0.01
 n2_from_bleed_air_psi_table=0.0:0.0, 4.0:8.0, 10.0:18.0, 18.0:25.0
 bleed_air_on_n2_tc=1
 N1_to_oil_pressure_table=0.00:0.00, 0.21:0.71, 0.90:0.99, 1.00:1.00
-//KH New
-itt_maxcorrection=400

 [JET_ENGINE]
 thrust_scalar=1.0
@@ -0,0 +1,415 @@
+//TFDi Designs MD-11 (GE CF6-80C2D1F) vRC1.4_GE 29AUG24
+
+[VERSION]
+//KH Bump major
+major = 4
+minor = 1
+
+[WEIGHT_AND_BALANCE]
+max_gross_weight = 630500
+empty_weight = 283975
+reference_datum_position = 0, 0, 0
+empty_weight_CG_position = -13.3, 0, 0
+CG_forward_limit = 0.13
+CG_aft_limit = 0.34
+empty_weight_pitch_MOI = 20897456
+empty_weight_roll_MOI = 6368555
+empty_weight_yaw_MOI = 40761191
+empty_weight_coupled_MOI = 0
+activate_mach_limit_based_on_cg = 0
+activate_cg_limit_based_on_mach = 0
+max_number_of_stations = 21
+station_load.0=170.0, 82.0, -2.0, 4.0, Pilot, 1
+station_load.1=170.0, 82.0, 2.0, 4.0, First Officer, 2
+station_load.2=170.0, 80.0, 3.3, 4.0, Engineer, 0
+station_load.3=850.0, 66.0, 0.0, 4.0, Cabin crew 1, 0
+station_load.4=680.0, 45.0, -5.0, 4.0, Business 1 L, 0
+station_load.5=680.0, 45.0, 0.0, 4.0, Business 1 C, 0
+station_load.6=680.0, 45.0, 5.0, 4.0, Business 1 R, 0
+station_load.7=680.0, 25.0, -5.0, 4.0, Business 2 L, 0
+station_load.8=960.0, 25.0, 0.0, 4.0, Business 2, 0
+station_load.9=680.0, 25.0, 5.0, 4.0, Business 2 R, 0
+station_load.10=5723.0, -20.0, -5.0, 4.0, Economy 3 L, 0
+station_load.11=5724.0, -20.0, 0.0, 4.0, Economy 3 C, 0
+station_load.12=5723.0, -20.0, 5.0, 4.0, Economy 3 R, 0
+station_load.13=9066.0, -50.0, -5.0, 4.0, Economy 4 L, 0
+station_load.14=9067.0, -50.0, 0.0, 4.0, Economy 4 C, 0
+station_load.15=9066.0, -50.0, 5.0, 4.0, Economy 4 R, 0
+station_load.16=850.0, -55.0, 0.0, 4.0, Cabin crew 2, 0
+station_load.17=28648.0, 20.0, 0.0, -4.5, Forward cargo, 0
+station_load.18=38648.0, -50.0, 0.0, -4.5, Rear cargo, 0
+station_load.19=0.0, 5.0, -2.0, -4.5, Left aux tank, 0
+station_load.20=0.0, 5.0, 2.0, -4.5, Right aux tank, 0
+
+[CONTACT_POINTS]
+static_pitch = 0
+static_cg_height = 16.4
+gear_system_type = 1
+tailwheel_lock = 0
+max_number_of_points = 16
+gear_locked_on_ground = 0
+gear_locked_above_speed = -1
+max_speed_full_steering = 25
+min_available_steering_angle_pct = 0.2
+max_speed_decreasing_steering = 45
+max_speed_full_steering_castering = 20
+max_speed_decreasing_steering_castering = 35
+emergency_extension_type = 2
+hyd_need_power_to_function= 0
+min_castering_angle = 0.05236
+max_castering_angle = 1.32645
+allow_stopped_steering = 1
+point.0 = 1, 64.22,  0.00,  -17.31, 1500, 0, 1.67, 70, 1.15, 1.183, 0.9, 49, 50,  0, 200, 250, 8
+point.1 = 1, -19.00, -17.15, -18.02, 2000, 1, 2.25, 0,  1.60, 1.116, 0.9, 51, 52, 2, 200, 250, 8
+point.2 = 1, -19.00,  17.15, -18.02, 2000, 2, 2.25, 0,  1.60, 1.116, 0.9, 53, 54, 3, 200, 250, 8
+point.3 = 1, -19.00,  0.00, -18.02, 2000, 2, 2.25, 0,  1.60, 1.116, 0.9, 55, 56, 3, 200, 250, 8
+point.4 = 2, 93.0, 0.0, 0.0, 787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 9, 0, 0, 1
+point.5 = 2, -48, -82.5, 0.0,  787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 5, 0, 0, 1
+point.6 = 2, -48, 82.5, 0.0,  787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 6, 0, 0, 1
+point.7 = 2, 60, 0.0, -9.0, 787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 9, 0, 0, 1
+point.8 = 2, 60, 0.0, 11.0, 787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 9, 0, 0, 1
+point.9 = 2, 16.0,  -28.00,  -13.00,  787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 4, 0, 0, 1
+point.10 = 2, 16.0,  28.00,  -13.00,  787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 4, 0, 0, 1
+point.11 = 2, -110.0,  0.00,  14.20,  787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 4, 0, 0, 1
+point.12 = 2, -96.5, 0.0, 42.0,  787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 4, 0, 0, 1
+point.13 = 2, -94.5, -30.0, 8.0,  787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 4, 0, 0, 1
+point.14 = 2, -94.5, 30.0, 8.05,  787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 4, 0, 0, 1
+point.15 = 2, -52.0, 0.0, -9.0, 787, 0, 0.0,  0, 0.00, 0.0, 0.0,  0.0,  0.0, 9, 0, 0, 1
+
+[FUEL]
+LeftMain = -20.00, -38.00, -1.0, 8938, 10.00
+RightMain = -20.00, 38.00, -1.0, 8938, 10.00
+Center1 = 10.00, 0.00, -3.00, 13000.00, 10.00
+Center2 = 0.00, 0.00, -1.00, 1642.00, 10.00
+Center3 = 0.00, 0.00, -5.00, 1958.00, 11.00
+LeftTip = -39.00, -74.00, 0.00, 884.00, 10.00
+RightTip = -39.00, 74.00, 0.00, 884.00, 10.00
+External1 = -70.00, 0.00, -3.00, 2371.00, 10.00
+LeftAux = 5.00, -2.00, -3.00, 0.00, 10.00
+RightAux = 5.00, 2.00, -3.00, 0.00, 10.00
+fuel_type = 2
+number_of_tank_selectors = 3
+electric_pump = 0
+engine_driven_pump = 1
+manual_transfer_pump = 0
+manual_pump = 0
+anemometer_pump = 0
+fuel_dump_rate = 0.02
+default_fuel_tank_selector.1 = 1
+default_fuel_tank_selector.2 = 1
+default_fuel_tank_selector.3 = 1
+
+[AIRPLANE_GEOMETRY]
+wing_area = 3648
+wing_span = 170.5
+wing_root_chord = 34.68
+wing_camber = 1
+wing_thickness_ratio = 0.04
+wing_dihedral = 2.2
+wing_incidence = 1.5
+wing_twist = -0.5
+oswald_efficiency_factor = 0.72
+wing_winglets_flag = 1
+wing_sweep = 38
+wing_pos_apex_lon = 1
+wing_pos_apex_vert = -2.75
+htail_area = 920
+htail_span = 59.5
+htail_pos_lon = -88
+htail_pos_vert = 5.0
+htail_incidence = 1.0
+htail_sweep = 36.5
+htail_thickness_ratio = 0.02
+vtail_area = 760
+vtail_span = 28.1
+vtail_sweep = 43
+vtail_pos_lon = -79.5
+vtail_pos_vert = 28
+vtail_thickness_ratio = 0.03
+fuselage_length = 202.16
+fuselage_diameter = 19.75
+fuselage_center_pos = 10, 0, 1.0
+elevator_area = 327.0
+aileron_area = 225.0
+rudder_area = 190.0
+elevator_up_limit = 25
+elevator_down_limit = 15
+aileron_up_limit = 25
+aileron_down_limit = 15
+rudder_limit = 23
+rudder_trim_limit = 20
+elevator_trim_limit = 30
+spoiler_limit = 60
+spoilerons_available = 1
+aileron_to_spoileron_gain = 3
+min_ailerons_for_spoilerons = 10
+min_flaps_for_spoilerons = 0
+spoiler_extension_time = 2.0
+spoiler_handle_available = 1
+auto_spoiler_available = 1
+auto_spoiler_auto_retracts = 0
+auto_spoiler_min_speed = 72
+positive_g_limit_flaps_up = 5.0
+positive_g_limit_flaps_down = 4.0
+negative_g_limit_flaps_up = -4.0
+negative_g_limit_flaps_down = -3.0
+load_safety_factor = 1.5
+elevator_trim_neutral = -1
+aileron_to_rudder_scale = 0
+flap_to_aileron_scale = 0
+fly_by_wire = 0
+elevator_elasticity_table = 0:0.7, 87:0.3, 210:0.5, 315:0.5
+aileron_elasticity_table = 0:1, 300:0.7, 600:0.6, 1500:0.4
+rudder_elasticity_table = 0:1, 300:1, 800:1, 1000:0.7
+elevator_trim_elasticity_table = 0:0.0, 87:0.40, 240:0.73, 290:0.73, 315:0.65
+
+[AERODYNAMICS]
+lift_coef_pitch_rate = -8.09375
+lift_coef_daoa = 2.68359
+lift_coef_delta_elevator = -0.28906
+lift_coef_horizontal_incidence = 0.91357
+lift_coef_flaps = 1
+lift_coef_spoilers = 0
+drag_coef_zero_lift = 0.02539
+drag_coef_flaps = 0.1889
+drag_coef_gear = 0.0549
+drag_coef_spoilers = 0.05207
+side_force_slip_angle = -0.77392
+side_force_roll_rate = 0.15478
+side_force_yaw_rate = 0
+side_force_delta_rudder = -0.14746
+pitch_moment_horizontal_incidence = 4.24755
+pitch_moment_delta_elevator = -2.37011
+pitch_moment_delta_trim = 0
+pitch_moment_pitch_damping = -214.1455
+pitch_moment_aoa_0 = 0.11621
+pitch_moment_daoa = 8.87255
+pitch_moment_flaps = 0.0105
+pitch_moment_gear = 0.0171
+pitch_moment_spoilers = 0
+pitch_moment_delta_elevator_propwash = -2.37011
+pitch_moment_pitch_propwash = -214.1455
+roll_moment_slip_angle = 0.26855
+roll_moment_roll_damping = -3.40869
+roll_moment_yaw_rate = -0.27197
+roll_moment_spoilers = 0.07
+roll_moment_delta_aileron = -0.34521
+roll_moment_delta_rudder = 0.01660
+roll_moment_delta_aileron_trim_scalar = -0.2
+yaw_moment_slip_angle = 0.18066
+yaw_moment_roll = 0.07470
+yaw_moment_yaw_damping = -2.49023
+yaw_moment_yaw_propwash = 0
+yaw_moment_delta_aileron = -0.00683
+yaw_moment_delta_rudder = 0.07128
+yaw_moment_delta_rudder_propwash = 0.075
+yaw_moment_delta_rudder_trim_scalar = -0.2
+compute_aero_center = 0
+aero_center_lift = -13.75
+lift_coef_aoa_table = -3.15:0, -2.36:0.36, -1.57:0, -0.349:-1.2, 0:-0.022, 0.0314:0.6082, 0.221:2.049, 0.244:1.15, 0.349:0.6, 0.436:0.5, 0.698:0.1, 1.963:-0.359, 2.356:-0.5, 3.15:0
+lift_coef_ground_effect_mach_table = 0.06:0.50, 0.227:0.80, 0.242:0.80, 0.3:1.0, 0.9:1.0
+lift_coef_mach_table = 0:1
+lift_coef_delta_elevator_mach_table = 0:0
+lift_coef_daoa_mach_table = 0:0
+lift_coef_pitch_rate_mach_table = 0:0
+lift_coef_horizontal_incidence_mach_table = 0:0
+drag_coef_zero_lift_mach_tab = 0:0, 0.4:0.0054, 0.5:0.0056, 0.6:0.0059, 0.7:0.0062, 0.72:0.0062, 0.8:0.0063, 0.8:0.0063, 1:0.0337, 1.2:0.0337, 1.4:0.0337, 1.6:0.0337, 1.8:0.0337, 2.0:0.0337
+side_force_slip_angle_mach_table = 0:0
+side_force_delta_rudder_mach_table = 0:0
+side_force_yaw_rate_mach_table = 0:0
+side_force_roll_rate_mach_table = 0:0
+pitch_moment_aoa_table = -3.15:-2, -0.436:-0.626, -0.262:-0.434, 0:-0.146, 0.087:-0.05, 0.175:0.0466, 0.262:0.143, 0.314:0.2, 0.349:0.1, 0.436:0.15, 3.15:2
+pitch_moment_delta_elevator_aoa_table = -180:-1, -40:0.05, -20:0.455, -14:0.853, -7:1.007, 0:1, 7:0.839, 14:0.693, 20:0.381, 40:-0.08, 180:-1
+pitch_moment_horizontal_incidence_aoa_table = 0:1
+pitch_moment_daoa_aoa_table = 0:1
+pitch_moment_pitch_alpha_table = 0:1
+pitch_moment_delta_elevator_mach_table = 0:0
+pitch_moment_daoa_mach_table = 0:0
+pitch_moment_pitch_rate_mach_table = 0:0
+pitch_moment_horizontal_incidence_mach_table = 0:0
+pitch_moment_aoa_0_mach_table = 0:0
+roll_moment_aoa_table = 0:0
+roll_moment_slip_angle_aoa_table = 0:1
+roll_moment_roll_rate_aoa_table = 0:1
+roll_moment_delta_aileron_aoa_table = 0:1
+roll_moment_slip_angle_mach_table = 0:0
+roll_moment_delta_rudder_mach_table = 0:0
+roll_moment_delta_aileron_mach_table = 0:0
+roll_moment_yaw_rate_mach_table = 0:0
+roll_moment_roll_rate_mach_table = 0:0
+yaw_moment_aoa_table = 0:0
+yaw_moment_slip_angle_aoa_table = 0:1
+yaw_moment_yaw_rate_aoa_table = 0:1
+yaw_moment_delta_rudder_aoa_table = 0:1
+yaw_moment_slip_angle_mach_table = 0:0
+yaw_moment_delta_rudder_mach_table = 0:0
+yaw_moment_delta_aileron_mach_table = 0:0
+yaw_moment_yaw_rate_mach_table = 0:0
+yaw_moment_roll_rate_mach_table = 0:0
+elevator_scaling_table = -0.785:0.5, -0.524:0.85, -0.175:0.9, 0:1, 0.175:0.9, 0.524:0.85, 0.785:0.5 ; 0:1
+aileron_scaling_table = -0.785:1 ;0:1
+rudder_scaling_table = -0.785:0.5, -0.524:0.6, -0.175:0.9, 0:1, 0.175:0.9, 0.524:0.6, 0.785:0.5 ;0:1
+elevator_lift_coef=3
+rudder_lift_coef=2
+aileron_load_factor_effectiveness_table = 0:1
+lift_coef_at_drag_zero = 0.1
+lift_coef_at_drag_zero_flaps = 0.1
+
+[FLIGHT_TUNING]
+modern_fm_only = 1
+cruise_lift_scalar = 1
+parasite_drag_scalar = 1
+induced_drag_scalar = 1
+flap_induced_drag_scalar = 1
+elevator_effectiveness = 1
+elevator_maxangle_scalar = 1
+aileron_effectiveness = 1
+rudder_effectiveness = 1
+rudder_maxangle_scalar = 1
+pitch_stability = 5
+roll_stability = 1
+yaw_stability = 1
+pitch_gyro_stability = 1
+roll_gyro_stability = 1
+yaw_gyro_stability = 1
+empty_CG_deviation_limit=0
+elevator_trim_effectiveness = 1.5
+aileron_trim_effectiveness = 1
+rudder_trim_effectiveness = 1
+hi_alpha_on_roll = 0
+hi_alpha_on_yaw = 0
+p_factor_on_yaw = 0
+torque_on_roll = 0
+gyro_precession_on_roll = 0
+gyro_precession_on_yaw = 0
+engine_wash_on_roll = 0
+wingflex_scalar = 2
+wingflex_offset = 0
+stall_coef_at_min_weight = 0.71
+ground_crosswind_effect_zero_speed=10
+ground_crosswind_effect_max_speed=150
+ground_high_speed_steeringwheel_static_friction_scalar=1.5
+ground_high_speed_otherwheel_static_friction_scalar=1.2
+//KH START New ground physics
+ground_new_contact_model_gear_flex=0.0044
+ground_new_contact_model_gear_flex_damping=7
+ground_new_contact_model_rolling_stickyness=0.15
+ground_new_contact_model_up_to_speed_lateral=10000
+ground_new_contact_model_up_to_speed_lateral_steering=10000
+ground_new_contact_model_up_to_speed_longitudinal=10000
+enable_high_accuracy_integration=1
+//KH END New ground physics
+
+[REFERENCE SPEEDS]
+full_flaps_stall_speed = 115
+flaps_up_stall_speed = 149
+cruise_speed = 345
+cruise_mach = 0.85
+crossover_speed = 330
+max_mach = 0.9
+max_indicated_speed = 395
+max_flaps_extended = 170
+normal_operating_speed = 325
+airspeed_indicator_max = 562.5
+rotation_speed_min = 145
+climb_speed = 250
+cruise_alt = 41000
+takeoff_speed = 165
+spawn_cruise_altitude = 5000
+spawn_descent_altitude = 500
+best_angle_climb_speed = 0
+approach_speed = 0
+best_glide = 0
+max_gear_extended = 200
+
+[INTERACTIVE POINTS]
+number_of_interactive_points = 3
+interactive_point.0 = 0.4, 44.15, -10.05, -0.567, 0, 0, 0, -88, 35, 0, 35, 0, 2, 33
+interactive_point.1 = 0.4, 75.65, 7.92, -0.567, 0, 0, 0, 75, 33, 2, 33, 2, 2, 33
+interactive_point.2 = 0.4, 33.3, 6.28, -7.18, 1, 0, 0, 90, 0, 0, 0, 0, 0, 0
+
+[STALL PROTECTION]
+stall_protection = 1		; Alpha Protection
+off_limit = 11.5				; Alpha below which the Alpha Protection can be disabled (If also below off_yoke_limit)
+off_yoke_limit = 1			; Yoke position (in Pct) below which the Alpha Protection can be disabled (If also below off_limit)
+on_limit = 14				; Alpha above which the Alpha Protection timer starts
+on_goal = 10				; The alpha that the Alpha Protection will attempt to reach when triggered
+timer_trigger = 1.0			; Duration (in Sec) the alpha must be above on_limit before the Alpha Protection is triggered
+
+;===================== FLAPS =====================
+
+[FLAPS.0]
+type = 1
+system_type = 1
+span-outboard = 0.25
+extending-time = 22
+damaging-speed = 250
+blowout-speed = 300
+maneuvering_flaps = 0
+lift_scalar = 0.1
+drag_scalar = 1.1
+pitch_scalar = 1
+max_on_ground_position = 7
+flaps-position.0 = 0, -1, 0
+flaps-position.1 = 0, -1, 0
+flaps-position.2 = 15, 215, 0.25, 2, 0, 0, 0, 1
+flaps-position.3 = 20, 200, 0.26, 1.1, 0, 0, 0, 1
+flaps-position.4 = 24, 185, 0.28, 1.0, 0, 0, 0, 1
+flaps-position.5 = 28, 177, 0.3, 1.0, 1.5, 0, 0
+flaps-position.6 = 35, 160, 0.32, 1
+flaps-position.7 = 50, 150, 0.33, 1
+
+
+[FLAPS.1]
+type = 1
+system_type = 1
+span-outboard = 0.7
+extending-time = 23
+damaging-speed = 250
+blowout-speed = 300
+maneuvering_flaps = 0
+lift_scalar = 0.1
+drag_scalar = 1.1 ;1.5
+pitch_scalar = 1
+max_on_ground_position = 7
+flaps-position.0 = 0, -1, 0
+flaps-position.1 = 0, -1, 0
+flaps-position.2 = 15, 215, 0.25, 2, 0, 0, 0, 1
+flaps-position.3 = 20, 200, 0.26, 1.1, 0, 0, 0, 1
+flaps-position.4 = 24, 185, 0.28, 1.0, 0, 0, 0, 1
+flaps-position.5 = 28, 177, 0.3, 1.0, 1.5, 0, 0
+flaps-position.6 = 35, 160, 0.32, 1
+flaps-position.7 = 50, 150, 0.33, 1
+
+
+[FLAPS.2]
+type = 2
+system_type = 1
+span-outboard = 0.2
+extending-time = 8
+damaging-speed = 250
+blowout-speed = 300
+maneuvering_flaps = 0
+lift_scalar = 0.5
+// KH From 0.75
+drag_scalar = 0.75
+pitch_scalar = 1
+max_on_ground_position = 1
+flaps-position.0 = 0, -1, 0
+flaps-position.1 = 20, 215, 0.5, 0.5
+
+[FLAPS.3]
+type = 2
+system_type = 1
+span-outboard = 0.7
+extending-time = 21
+damaging-speed = 250
+blowout-speed = 300
+maneuvering_flaps = 0
+lift_scalar = 0.5
+// KH From 0.75
+drag_scalar = 0.75
+pitch_scalar = 1
+max_on_ground_position = 1
+flaps-position.0 = 0, -1, 0
+flaps-position.1 = 20, 215, 0.5, 0.5
@@ -0,0 +1,290 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": 9,
+   "id": "21e0ac82",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import numpy as np\n",
+    "import matplotlib.pyplot as plt\n",
+    "import pandas as pd"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 10,
+   "id": "653d2200",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Read in data\n",
+    "dataNp = np.loadtxt(open(\"N2onMachTest.csv\", \"r\"), delimiter=\",\", skiprows=1)\n",
+    "dataNp = np.delete(dataNp, 1, 1)\n",
+    "dataNp = np.delete(dataNp, 1, 1)\n",
+    "dataNp = np.delete(dataNp, 1, 1)\n",
+    "dataNp = np.delete(dataNp, 1, 1)\n",
+    "dataNp = np.delete(dataNp, 1, 1)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "ecfc5eed",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Constants\n",
+    "max_thrust = 62000\n",
+    "max_EPR = 1.6432\n",
+    "idle_thrust = 2170\n",
+    "idle_epr = 1.02282\n",
+    "k = 0.3\n",
+    "\n",
+    "n1s = np.linspace(20,110,19)\n",
+    "n1s = np.insert(n1s, 0, 0)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 12,
+   "id": "ddbebf95",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Thrust function\n",
+    "thrust = lambda epr, mach: (idle_thrust + (max_thrust - idle_thrust) * ((epr - idle_epr) / (max_EPR - idle_epr))) * (1 - k * mach)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 30,
+   "id": "c65ab937",
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAHHCAYAAABtF1i4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABMEUlEQVR4nO3deVhU9cMF8DMzMMO+CYIoCO47qxCarZSmUdqm4m5qrqm8LeKapmKZ/jQ1LTXNXMvMyswyyswiERDLDTcQF0D2QbaBmfv+gU5OgjII3FnO53nmEe7cO3PmAjPHu3yvRBAEAUREREQikYodgIiIiMwbywgRERGJimWEiIiIRMUyQkRERKJiGSEiIiJRsYwQERGRqFhGiIiISFQsI0RERCQqlhEiIiISFcsIEYnCx8cHzz77rNgxiMgAsIwQmZnNmzdDIpFAIpHgyJEjd90vCAK8vLwgkUgMvix8++23CAwMhJWVFby9vTFv3jxUVlbWalmNRoP3338fvr6+sLKyQrdu3bBjx44GTkxE1WEZITJTVlZW2L59+13Tf/vtN1y9ehUKhUKEVLX3ww8/oH///nBycsKqVavQv39/LFy4EFOmTKnV8rNmzcLbb7+Np556CqtWrYK3tzciIyOxc+fOBk5ORP8l4YXyiMzL5s2bMWrUKLzwwgs4fPgwMjIyYGFhob1/3LhxSEpKQk5ODrp06YJ9+/Y1SA4fH58HevzOnTvD0tISCQkJ2vyzZ8/G4sWLcfr0aXTo0KHGZa9duwZfX1+MGzcOq1evBlC1RejRRx9Famoq0tLSIJPJ6pSLiPTHLSNEZmrw4MHIzc3FwYMHtdNUKhV2796NyMjIapf54IMP0KNHDzRp0gTW1tYICgrC7t27q51369atCAkJgY2NDZydnfHII4/gp59+umu+I0eOICQkBFZWVmjVqhW2bNly3+ynT5/G6dOnMW7cOJ0iNXHiRAiCUGOm27755htUVFRg4sSJ2mkSiQQTJkzA1atXERcXd98MRFR/WEaIzJSPjw/CwsJ0jpP44YcfUFhYiEGDBlW7zMqVKxEQEIAFCxZg8eLFsLCwwMsvv4zvv/9eZ7758+dj2LBhsLS0xIIFCzB//nx4eXnhl19+0ZnvwoULeOmll/DUU09h2bJlcHZ2xsiRI3Hq1Kl7Zj9+/DgAIDg4WGe6p6cnWrRoob3/Xsvb2tqiY8eOOtNDQkJ0Hp+IGofF/WchIlMVGRmJ6OholJaWwtraGtu2bcOjjz4KT0/Pauc/d+4crK2ttd9PnjwZgYGBWL58Ofr16wegqmAsWLAAAwYMwO7duyGV/vt/nv/uFU5JScHhw4fRq1cvAMArr7wCLy8vbNq0CR988EGNuTMyMgAAzZo1u+u+Zs2a4fr16/d83RkZGXB3d4dEIrlrWQD3XZ6I6he3jBCZsVdeeQWlpaXYt28fioqKsG/fvhp30QDQKSL5+fkoLCxEr169kJSUpJ2+d+9eaDQazJ07V6eIALjrw79Tp07aIgIAbm5uaN++PS5dunTP3KWlpQBQ7UG2VlZW2vvvtXxNy975+ETUOLhlhMiMubm5ITw8HNu3b0dJSQnUajVeeumlGufft28fFi5ciOTkZJSXl2un31kyLl68CKlUik6dOt33+b29ve+a5uzsjPz8/Hsud7sU3ZnhtrKyMp3SVNPyNS175+MTUePglhEiMxcZGYkffvgB69atwzPPPAMnJ6dq5/v999/x3HPPwcrKCh999BH279+PgwcPIjIy8q7dL7VV0xkr93u827tTbu+uuVNGRkaNu5nuXD4zM/Ou57n9ePdbnojqF8sIkZkbMGAApFIp/vrrr3vuovnqq69gZWWFH3/8EaNHj8YzzzyD8PDwu+Zr3bo1NBoNTp8+3WCZ/f39AQAJCQk6069fv46rV69q77/X8iUlJThz5ozO9KNHj+o8PhE1DpYRIjNnZ2eHtWvX4p133kFERESN88lkMkgkEqjVau20tLQ07N27V2e+/v37QyqVYsGCBdBoNDr31dewRp07d0aHDh3wySef6ORZu3YtJBKJzq6mwsJCnD17FoWFhdppzz//PCwtLfHRRx/pZFu3bh2aN2+OHj161EtOIqodHjNCRBgxYsR95+nXrx+WL1+OPn36IDIyEjdu3MCaNWvQpk0b/P3339r52rRpg1mzZuHdd99Fr1698MILL0ChUODYsWPw9PRETExMvWReunQpnnvuOTz99NMYNGgQTp48idWrV2PMmDE6p+x+/fXXGDVqFDZt2oSRI0cCAFq0aIFp06Zh6dKlqKioQPfu3bF37178/vvv2LZtGwc8I2pkLCNEVCtPPPEENm7ciCVLlmDatGnw9fXFe++9h7S0NJ0yAgALFiyAr68vVq1ahVmzZsHGxgbdunXDsGHD6i3Ps88+iz179mD+/PmYMmUK3NzcMHPmTMydO7dWyy9ZsgTOzs74+OOPsXnzZrRt2xZbt269564qImoYHA6eiIiIRMVjRoiIiEhULCNEREQkKpYRIiIiEhXLCBEREYmKZYSIiIhExTJCREREojKKcUY0Gg2uX78Oe3v7u676SURERIZJEAQUFRXB09Pzrqt438koysj169fh5eUldgwiIiKqgytXrqBFixY13m8UZcTe3h5A1YtxcHAQOQ0RERHVhlKphJeXl/ZzvCZGUUZu75pxcHBgGSEiIjIy9zvEggewEhERkahYRoiIiEhULCNEREQkKqM4ZqQ2NBoNVCqV2DGogcjl8nueFkZERMbLJMqISqVCamoqNBqN2FGogUilUvj6+kIul4sdhYiI6pnRlxFBEJCRkQGZTAYvLy/+79kE3R70LiMjA97e3hz4jojIxBh9GamsrERJSQk8PT1hY2MjdhxqIG5ubrh+/ToqKythaWkpdhwiIqpHRr8ZQa1WAwA335u42z/f2z9vIiIyHXqXkcOHDyMiIgKenp6QSCTYu3fvfZc5dOgQAgMDoVAo0KZNG2zevLkOUe+Nm+5NG3++RESmS+8yUlxcDD8/P6xZs6ZW86empqJfv354/PHHkZycjGnTpmHMmDH48ccf9Q5LREREpkfvMvLMM89g4cKFGDBgQK3mX7duHXx9fbFs2TJ07NgRkydPxksvvYT//e9/eoc1J4cOHYJEIkFBQYHYUe4pNjYWHTt2fKDdJwcOHIC/vz/PhiIiMlMNfsxIXFwcwsPDdab17t0bcXFxNS5TXl4OpVKpczMlEonknrd33nlH1Hw+Pj5YsWJFreZ96623MHv2bMhkMgDA8ePHERAQADs7O0RERCAvL087b2VlJYKCghAfH6/zGH369IGlpSW2bdtWb6+BiIiMR4OXkczMTLi7u+tMc3d3h1KpRGlpabXLxMTEwNHRUXvz8vJq6JiNKiMjQ3tbsWIFHBwcdKa98cYbdXrcxh707ciRI7h48SJefPFF7bQxY8bgiSeeQFJSEgoLC7F48WLtfcuWLUPPnj0REhJy12ONHDkSH374YaPkJiIyZkVlFdBoBLFj1CuDPJsmOjoahYWF2tuVK1fEjlSvPDw8tDdHR0dIJBKdaXZ2dtp5ExMTERwcDBsbG/To0QMpKSna+9555x34+/tjw4YN8PX1hZWVFYDqt2z4+/trt7gIgoB33nkH3t7eUCgU8PT0xOuvvw4AeOyxx3D58mVMnz5du6WmJjt37sRTTz2lfV4AOHPmDMaOHYt27dph8ODBOHPmDADg0qVL2LhxIxYtWlTtY0VERCAhIQEXL16s/YokIjIz24+mw3/BQYTGxGLW1//g9/PZqFAb/y7uBh9nxMPDA1lZWTrTsrKy4ODgAGtr62qXUSgUUCgUdXo+QRBQWiHO6Z/WlrJ6P+tj1qxZWLZsGdzc3DB+/HiMHj0af/zxh/b+Cxcu4KuvvsKePXu0u0ru56uvvsL//vc/7Ny5E507d0ZmZiZOnDgBANizZw/8/Pwwbtw4jB079p6P8/vvvyMyMlJnmp+fHw4ePIg2bdogNjYW3bp1AwCMHz8e77//Puzt7at9LG9vb7i7u+P3339H69ata/U6iIjMyYbfL2Hh91X/wcsuKse2o+nYdjQdDlYWCO/kjt6dPfBIWzdYy2v3WWBIGryMhIWFYf/+/TrTDh48iLCwsAZ5vtIKNTrNFedMndMLesNGXr+rdNGiRXj00UcBADNmzEC/fv1QVlam3RqhUqmwZcsWuLm51fox09PT4eHhgfDwcFhaWsLb21u768TFxQUymQz29vbw8PC45+NcvnwZnp6eOtM2bNiAiRMn4oMPPkDPnj0RHR2Nzz//HDY2NujevTt69+6NixcvYtCgQVi4cKHOsp6enrh8+XKtXwcRkblY/ct5fPDTOQDAa4+0Qs82rjhwKhM/ncpEzk0V9iRdw56ka7C2lOGx9m7o08UDj3doCgcr4xgkUu9Pzps3b+LChQva71NTU5GcnAwXFxd4e3sjOjoa165dw5YtWwBU/Y949erVeOuttzB69Gj88ssv+OKLL/D999/X36swYbe3LABAs2bNAAA3btyAt7c3AKBly5Z6FREAePnll7FixQq0atUKffr0Qd++fREREQELC/1+HUpLS3V20QBA586d8dtvv2m/z83Nxbx583D48GFMmTIFPXr0wJ49e9C9e3eEhoYiIiJCO6+1tTVKSkr0ykBEZMoEQcDSH1Pw0aGqXdhRT7XDlCfaQCKR4JF2bnj3+S5ISs/HgZOZOHAyE9cKSvHDyUz8cDITljIJerZxRZ/OHniqkzua2NVtj0Nj0LuMJCQk4PHHH9d+HxUVBQAYMWIENm/ejIyMDKSnp2vv9/X1xffff4/p06dj5cqVaNGiBTZs2IDevXvXQ/y7WVvKcHpBwzx2bZ67vt059PntXUB3ngJra2t71zJSqRSCoHtwU0VFhfZrLy8vpKSk4Oeff8bBgwcxceJELF26FL/99pteQ627uroiPz//nvNERUVh2rRpaNGiBQ4dOoSLine truncated
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAHHCAYAAABtF1i4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABPUElEQVR4nO3dd1zU9eMH8NcN7thDpixxL1QQBNHMSvqamWXlNldpmSOVluZMc5RlljnSTMttrsrKMsrSXy5AVBT3wsFUOPbB3fv3h3p5iQoIfG68no8HD+Fz78/d6z6Oe/kZ749MCCFAREREJBG51AGIiIjIurGMEBERkaRYRoiIiEhSLCNEREQkKZYRIiIikhTLCBEREUmKZYSIiIgkxTJCREREkmIZISIiIkmxjBCRkaCgIDzzzDNSxyAiK8IyQmQiVq5cCZlMBplMhj179tz1uBACAQEBkMlkJl8WfvjhB7Ru3Rq2trYIDAzE1KlTUVpa+sD1Tpw4gXfeeQchISFwcnJC7dq10bVrV8TFxZU5/sqVK+jVqxdcXV3h7OyM5557DufOnStzbFpaGl577TX4+fnB1tYWQUFBeOWVV+4a9/vvv+Pxxx+Hh4cHXF1dERERgVWrVpX5nMuXL0fTpk1ha2uLhg0bYsGCBXeN2bp1Kzp37gxfX1+o1Wr4+/ujR48eSEpKeuD2ILIWSqkDEJExW1tbrF27Fo888ojR8r/++guXL1+GWq2WKFn5/PLLL+jevTsee+wxLFiwAEePHsUHH3yA9PR0LF68+L7rfvXVV1i+fDlefPFFjBgxAjk5Ofjyyy/Rtm1b7NixA9HR0YaxeXl5ePzxx5GTk4P33nsPNjY2+PTTT9GxY0ckJibC3d3dMDYlJQXt27cHAAwfPhx+fn64evUqDhw4YPT6P/zwA7p3746oqChMmzYNMpkMGzduxMCBA5GZmYlx48YZxn755ZcYPnw4XnzxRcTExGD37t144403UFBQgHfffdcw7ujRo3Bzc8OYMWPg4eGB1NRUfP3114iIiMDevXvRqlWrh9reRBZBEJFJWLFihQAgXnjhBeHh4SFKSkqMHh82bJgICwsTderUEV27dq22HA/7/M2aNROtWrUyyj9x4kQhk8lEcnLyfdeNi4sTubm5RssyMzOFp6enaN++vdHyDz/8UAAQBw4cMCxLTk4WCoVCTJgwwWhsly5dRN26dUVmZuZ9X//JJ58Uvr6+oqioyLCspKRE1K9fX7Rs2dKwrKCgQLi7u9+1nfr37y8cHBzE9evX7/s6qampQqlUitdee+2+44isBQ/TEJmYvn37IisrCzt37jQs02q12LRpE/r161fmOh9//DHatWsHd3d32NnZISwsDJs2bSpz7OrVqxEREQF7e3u4ubnh0UcfxW+//XbXuD179iAiIgK2traoV68evv322wdmP378OI4fP45XX30VSuW/O15HjBgBIcQ9M90WFhYGR0dHo2Xu7u7o0KEDkpOTjZZv2rQJbdq0QZs2bQzLmjRpgk6dOmHjxo2GZSdOnMAvv/yCt99+G+7u7igqKkJJSUmZr6/RaODm5ma090mpVMLDwwN2dnaGZX/++SeysrIwYsQIo/VHjhyJ/Px8/PTTT/d9n15eXrC3t0d2dvZ9xxFZC5YRIhMTFBSEqKgorFu3zrDsl19+QU5ODvr06VPmOp999hlCQ0Mxffp0zJo1C0qlEj179rzrQ/H999/HgAEDYGNjg+nTp+P9999HQEAA/vjjD6NxZ86cQY8ePfDkk0/ik08+gZubGwYPHoxjx47dN/uhQ4cAAOHh4UbLfX194e/vb3i8olJTU+Hh4WH4Wa/X48iRI3e9DgBERETg7NmzyM3NBXDzHBAA8Pb2RqdOnWBnZwc7Ozt06dIFFy5cMFr3sccew7FjxzB58mScOXMGZ8+exYwZMxAXF4d33nnnge8zLCwMcrm8zPeZnZ2NjIwMHD16FEOHDoVGo0GnTp0qtT2ILI7Uu2aI6Kbbh2kOHjwovvjiC+Hk5CQKCgqEEEL07NlTPP7440KIsg+j3B53m1arFcHBweKJJ54wLDt9+rSQy+Xi+eefFzqdzmi8Xq83fF+nTh0BQPz999+GZenp6UKtVos333zzvu9h7ty5AoC4dOnSXY+1adNGtG3b9r7rl+Xvv/8WMplMTJ482bAsIyNDABDTp0+/a/zChQsFAHHixAkhhBBvvPGGACDc3d3FU089JTZs2CDmzp0rHB0dRf369UV+fr5h3by8PNGrVy8hk8kEAAFA2Nvbi23bthm9xsiRI4VCoSgzr6enp+jTp89dyxs3bmx4TkdHRzFp0qS7fh+IrBX3jBCZoF69eqGwsBDbt29Hbm4utm/ffs9DNACMDiHcuHEDOTk56NChAxISEgzLt23bBr1ejylTpkAuN/6rL5PJjH5u1qwZOnToYPjZ09MTjRs3vueVKrcVFhYCQJkn2dra2hoeL6/09HT069cPdevWNdoz8aDXuXNMXl4eAMDHxwc//fQTevXqhbfeegvLli3D2bNnsXbtWsO6arUajRo1Qo8ePbBu3TqsXr0a4eHheOmll7Bv3z6j11epVGVmvtf7XLFiBXbs2IFFixahadOmKCwshE6nq9D2ILJUvJqGyAR5enoiOjoaa9euRUFBAXQ6HXr06HHP8du3b8cHH3yAxMREFBcXG5bfWTLOnj0LuVyOZs2aPfD1AwMD71rm5uaGGzdu3He926Xozgy3FRUVGZWmB8nPz8czzzyD3Nxc7Nmzx+hckge9zp1jbv/aq1cvoxLWs2dPDBgwAP/88w+GDh0KABg1ahT27duHhIQEw9hevXqhefPmGDNmDPbv3294Tq1WW2bue73PqKgow/d9+vRB06ZNAdw834fI2nHPCJGJ6tevH3755RcsWbIEXbp0gaura5njdu/ejWeffRa2trZYtGgRfv75Z+zcuRP9+vWDEKJSr61QKMpc/qDnq127NgDg2rVrdz127do1+Pr6luv1tVotXnjhBRw5cgTff/89goODjR6vVasW1Gr1PV8HgOG1bv/q7e1tNE6hUMDd3d1QsLRaLZYvX46uXbsalRYbGxt06dIFcXFxhgJSu3Zt6HQ6pKen35U7Kyvrge/Tzc0NTzzxBNasWfPAbUFkDVhGiEzU888/D7lcjn379t33EM3mzZtha2uLX3/9FS+//DK6dOliNB/HbfXr14der8fx48erLXNISAgA3DVJ2dWrV3H58mXD4/ej1+sxcOBAxMbGYu3atejYseNdY+RyOVq0aFHmZGj79+9HvXr14OTkBODmSaXAzQnS7qTVapGZmQlPT08AQFZWFkpLS8s8dFJSUgK9Xm947F7vMy4uDnq9vlzvs7CwEDk5OQ8cR2QNWEaITJSjoyMWL16MadOmoVu3bvccp1AoIJPJjD5EL1y4gG3bthmN6969O+RyOaZPnw69Xm/0WGX3oPxX8+bN0aRJEyxdutQoz+LFiyGTyYwONeXk5ODEiRN3fSCPHj0aGzZswKJFi/DCCy/c87V69OiBgwcPGhWCkydP4o8//kDPnj0Nyx577DF4eXlhzZo1hkM4wM0Zb3U6HZ588kkANy+3dXV1xdatW40OweTl5eHHH39EkyZNDIdfnnjiCdSqVeuuSdwWL14Me3t7dO3a1bDsv3tPgJu/P7GxsWVeDURkjXjOCJEJGzRo0APHdO3aFfPmzcNTTz2Ffv36IT09HQsXLkSDBg1w5MgRw7gGDRpg4sSJmDFjBjp06IAXXngBarUaBw8ehK+vL2bPnl0lmefOnYtnn30W//vf/9CnTx8kJSXhiy++wNChQw3nSQA3p0kfMmQIVqxYgcGDBwMA5s+fj0WLFiEqKgr29vZYvXq10XM///zzcHBwAHBz7pJly5aha9eueOutt2BjY4N58+bB29sbb775pmEdtVqNuXPnYtCgQXj00UcxYMAAXLp0CZ999plhOwA3S91bb72FSZMmoW3bthg4cCB0Oh2WL1+Oy5cvG2Wxs7PDjBkzMHLkSPTs2ROdO3fG7t27sXr1asycORO1atUyjG3RogU6deqEkJAQuLm54fTp01i+fDlKSkowZ86cKtnmRGZP4qt5iOiWOy/tvZ+yLu1dvny5aNiwoVCr1aJJkyZixYoVYurUqaKsv+Jff/21CA0NFWq1Wri5uYmOHTuKnTt33vf5hRCiY8eOomPHjuV6L1u3bhUhISFCrVYLf39/MWnSJKHVast8vytWrDAsGzRokOHy17K+zp8/b/QcKSkpokePHsLZ2Vk4OjqKZ555Rpw+fbrMTOvWrROtWrUSarVaeHt7i1GjRgmNRnPXuDVr1oiIiAjh6uoq7OzsRGRkpNi0aVOZz7l06VLRuHFjoVKpRP369cWnn35qdJm0EEJMnTpVhIeHCzc3N6FUKoWvr6/o06ePOHLkSDm2JJF1kAlRRftniYiIiCqB54wQERGRpFhGiIiISFIsI0RERCQplhEiIiKSFMsIERERSYplhIiIiCRlFpOe6fV6XL16FU5OTnfdXZSIiIhMkxACubm58PX1vetu4XcyizJy9epVBAQESB2DiIiIKiElJQX+/v73fNwsysjtG16lpKTA2dlZ4jRERERUHhqNBgEBAYbP8XsxizJy+9CMs7MzywgREZGZedApFjyBlYiIiCTFMkJERESSYhkhIiIiSZnFOSPlodfrodVqpY5B1USlUt33sjAiIjJfFlFGtFotzp8/D71eL3UUqiZyuRx169aFSqWSOgoREVUxsy8jQghcu3YNCoUCAQEB/N+zBbo96d21a9cQGBjIie+IiCyM2ZeR0tJSFBQUwNfXF/b29lLHoWri6emJq1evorS0FDY2NlLHISKiKmT2uxF0Oh0AcPe9hbv9+3v795uIiCyH2ZeR27jr3rLx95eIyHJZTBkhIiIi88QyYqJ27doFmUyG7OxsqaPcV2xsLJo2bfpQh0927NiBkJAQXg1Line truncated
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAHHCAYAAABtF1i4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABPfElEQVR4nO3dd1hT5+MF8JMEkrD3VBBXXSggS7StHbRqLa0d7oG22lpx1W+H1r3tz2qtu+7W3aG2tWpraW21dQHuPVBcgIAShhBI3t8faGoqKihwE3I+z5NHuHlvcnJRcry5970yIYQAERERkUTkUgcgIiIiy8YyQkRERJJiGSEiIiJJsYwQERGRpFhGiIiISFIsI0RERCQplhEiIiKSFMsIERERSYplhIiIiCTFMkJERgICAvDyyy9LHYOILAjLCJGJWLFiBWQyGWQyGXbt2nXP/UII+Pn5QSaTmXxZ+PHHH9G8eXOo1Wr4+/tj7NixKC4ufuh6J0+exEcffYTg4GA4ODjAx8cH7du3R0JCQqnjf/vtNzz77LNwd3eHs7MzIiIisHLlylLHLl26FI0aNYJarUb9+vUxZ86c++ZYv349oqKiYGdnB2dnZ7Rs2RK///77PePS0tLw7rvvokaNGlCr1QgICMDbb79tNObUqVN4//330bJlS6jVashkMly4cOGh24LIkrCMEJkYtVqNNWvW3LP8zz//xOXLl6FSqSRIVXZbt25Fhw4d4OzsjDlz5qBDhw6YNGkSBg0a9NB1lyxZgsWLFyMsLAwzZszAsGHDcOrUKbRo0QK//fab0dgff/wRL774IrRaLcaNG4fJkyfDxsYGvXr1wueff2409ssvv0Tfvn3RpEkTzJkzB1FRURg8eDA+/fTTezKMGzcOXbt2hZ+fH2bOnIlJkyahWbNmuHLlitG4S5cuITw8HFu3bkX//v0xf/589O3bF9evXzcat3v3bsyePRs5OTlo1KhRWTcjkWURRGQSli9fLgCI119/Xbi7u4uioiKj+/v16ydCQ0NFrVq1RPv27Sstx+M+fuPGjUVQUJBR/pEjRwqZTCZOnDjxwHUTEhJETk6O0bKMjAzh4eEhWrVqZbT8hRdeEL6+vqKgoMCwrKioSNStW1c0a9bMsCw/P1+4ubnd85q6d+8u7OzsRFZWlmHZ7t27hUwmEzNnznzo62zXrp2oXbu2yMjIeOC4zMxModFohBBCTJ8+XQAQycnJD318IkvCPSNEJqZr167IzMzE9u3bDcu0Wi2+++47dOvWrdR1PvvsM7Rs2RJubm6wsbFBaGgovvvuu1LHrlq1ChEREbC1tYWLiwuefvpp/Prrr/eM27VrFyIiIqBWq1GnTh18/fXXD81+/PhxHD9+HO+88w6srKwMywcMGAAhxH0z3REaGgp7e3ujZW5ubnjqqadw4sQJo+UajQYuLi5Ge4qsrKzg7u4OGxsbw7I//vgDmZmZGDBggNH6cXFxyMvLw88//2xYNmvWLHh7e2PIkCEQQiA3N7fUnCdPnsTWrVvx4Ycfws3NDQUFBSgqKip1rKurKxwcHB74uoksHcsIkYkJCAhAVFQU1q5da1i2detWZGdno0uXLqWu88UXXyAkJAQTJkzAlClTYGVlhY4dOxq90QLA+PHj0bNnT1hbW2PChAkYP348/Pz87jke4uzZs3jzzTfxwgsvYMaMGXBxcUHv3r1x7NixB2Y/cOAAACAsLMxoua+vL2rWrGm4v7xSU1Ph7u5utOyZZ57BsWPHMHr0aJw9exbnzp3DxIkTkZCQgI8++uihmUJDQyGXy40yxcfHIzw8HLNnz4aHh4fhuJW5c+carXvnIyMvLy88//zzsLGxgY2NDdq1a8fjQYgehdS7ZoioxJ2Pafbv3y/mzp0rHBwcRH5+vhBCiI4dO4pnn31WCFH6xyh3xt2h1WpFYGCgeO655wzLzpw5I+RyuXjttdeETqczGq/X6w1f16pVSwAQf/31l2FZenq6UKlU4n//+98DX8OdjyFSUlLuuS88PFy0aNHigeuX5q+//hIymUyMHj3aaHlubq7o1KmTkMlkAoAAIGxtbcWmTZuMxsXFxQmFQlHqY3t4eIguXboIIYTIysoSAISbm5uwt7cX06dPF+vXrxdt27YVAMTChQsN6w0ePNgwtm3btmL9+vVi+vTpwt7eXtStW1fk5eWV+nz8mIaodNwzQmSCOnXqhFu3bmHz5s3IycnB5s2b7/sRDQCjjyVu3LiB7OxsPPXUU0hKSjIs37RpE/R6PcaMGQO53PifvkwmM/q+cePGeOqppwzfe3h4oEGDBjh//vwDc9+6dQsASj3IVq1WG+4vq/T0dHTr1g21a9c22ttx5zmeeOIJvPnmm1i7di1WrVqFsLAw9OjRA3v27DHKpFQqS338uzPd+UgmMzMTS5YswQcffIBOnTrh559/RuPGjTFp0iTDenfGent74+eff0anTp3wwQcfYPHixTh37lypByAT0f2xjBCZIA8PD0RHR2PNmjXYsGEDdDod3nzzzfuO37x5M1q0aAG1Wg1XV1d4eHhgwYIFyM7ONow5d+4c5HI5Gjdu/NDn9/f3v2eZi4sLbty48cD17pSiwsLCe+4rKCgwKk0Pk5eXh5dffhk5OTn44Ycf7jmWZODAgfjpp5+wbt06dOnSBd27d8dvv/0GHx8fDBkyxCiTVqst9TnuznTnT2tra6NtLZfL0blzZ1y+fBkpKSlGYzt16mRU7Dp27AgrKyv8888/ZX6dRMQyQmSyunXrhq1bt2LhwoVo164dnJ2dSx23c+dOvPLKK1Cr1Zg/fz62bNmC7du3o1u3bhBCPNJzKxSKUpc/7PF8fHwAANeuXbvnvmvXrsHX17dMz6/VavH666/j8OHD+OGHHxAYGHjP/UuXLkX79u2NyoC1tTXatWuHhIQEQwHx8fGBTqdDenr6PY+RmZlpyOTq6gq1Wg03N7d7Xr+npycAGMrYnXW8vLyMxikUCri5uT20tBGRMZYRIhP12muvQS6XY8+ePQ/8iOb777+HWq3GL7/8grfeegvt2rVDdHT0PePq1q0LvV6P48ePV1rm4OBgALhnkrKrV6/i8uXLhvsfRK/Xo1evXoiPj8eaNWvQunXre8ZkZmaiuLgYOp3unvuKioqg1+sN990vU0JCAvR6veF+uVyO4OBgXL9+/Z49KVevXgVQsscKKDn4FcA9c49otVpkZGQYxhFR2bCMEJkoe3t7LFiwAOPGjUNMTMx9xykUCshkMqM35gsXLmDTpk1G4zp06AC5XI4JEyZAr9cb3feoe1D+q0mTJmjYsCEWLVpklGfBggWQyWRGH39kZ2fj5MmTRh8lAcCgQYOwfv16zJ8/H6+//nqpz+Pp6QlnZ2ds3LjRqDjk5ubip59+QsOGDQ0fpTz33HNwdXXFggULjB5jwYIFsLW1Rfv27Q3LOnfuDJ1Oh6+++sqwrKCgAKtXr0bjxo0Ne0SeeeYZeHp6YvXq1SgoKDCMXbFiBXQ6HV544YUybzMiAqwePoSIpBIbG/vQMe3bt8fMmTPRtm1bdOvWDenp6Zg3bx7q1auHw4cPG8bVq1cPI0eOxMSJE/HUU0/h9ddfh0qlwv79++Hr64upU6dWSObp06fjlVdewYsvvoguXbrg6NGjmDt3Lvr27Ws0A+nGjRvRp08fLF++HL179wZQMs/H/PnzERUVBVtbW6xatcrosV977TXY2dlBoVDggw8+wKhRo9CiRQv06tULOp0OS5cuxeXLl43Ws7GxwcSJExEXF4eOHTuiTZs22LlzJ1atWoXJkyfD1dXVMPbdd9/FkiVLEBcXh9OnT8Pf3x8rV67ExYsX8dNPPxnGqVQqTJ8+HbGxsXj66afRs2dPpKSk4IsvvjBs2zuys7MNU8///fffAIC5c+fC2dkZzs7OGDhwYIVsdyKzJvHZPER0292n9j5Iaaf2Ll26VNSvX1+oVCrRsGFDsXz5cjF27FhR2j/xZcuWiZCQEKFSqYSLi4to3bq12L59+wMfXwghWrduLVq3bl2m17Jx40YRHBwsVCqVqFmzphg1apTQarWlvt7ly5cblsXGxhpO0y3t9t9TYlevXi0iIiKEs7OzsLGxEZGRkeK7774rNdOiRYtEgwYNhFKpFHXr1hWff/650SnNd6SlpYnY2Fjh6uoqVCqViIyMFNu2bSv1MdeuXSuCgoKESqUSXl5eYuDAgYbZVu9ITk6+7+upVavWwzcmkQWQCVFB+2eJiIiIHgGPGSEiIiJJsYwQERGRpFhGiIiISFIsI0RERCQplhEiIiKSFMsIERERScosJj3T6/W4evUqHBwc7rm6KBEREZkmIQRycnLg6+t7z9XC72YWZeTq1avw8/OTOgYRERE9gkuXLqFmzZr3vd8syoiDgwOAkhfj6OgocRoiIiIqC41GAz8/P8P7+P2YRRm589GMo6MjywgREZGZedghFjyAlYiIiCTFMkJERESSYhkhIiIiSZnFMSNlodfrodVqpY5BlUSpVD7wtDAiIjJf1aKMaLVaJCcnQ6/XSx2FKolcLkft2rWhVCqljkJERBXM7MuIEALXrl2DQqGAn58f//dcDd2Z9O7atWvw9/fnxHdERNWM2ZeR4uJi5Ofnw9fXF7a2tlLHoUri4eGBq1evori4GNbW1lLHISKiCmT2uxF0Oh0AcPd9NXfn53vn501ERNWH2ZeRO7jrvnrjz5eIqPqqNmWEiIiIzBPLiInasWMHZDIZbt68KXWUB4qPj0ejRo0e6+OTbdu2ITg4mGdDERFZKJYRCchksgfexo0bJ2m+gIAAzJo1q0xjP/roI4waNQoKhQIAcODAAYSEhMDe3h4xMTHIysoyjC0uLkZoaCj27dtn9Bht27aFtbU1Vq9eXWGvgYiIzAfLiASuXbtmuM2aNQuOjo5Gyz744INHetyqnvRt165dOHfuHN544w3Dsr59++K5555DUlISsrOzMWXKFMN9M2bLine truncated
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAHHCAYAAABtF1i4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABMu0lEQVR4nO3deVhU1f8H8PfMwMywg6yiLO67rIJoZgulZZSWprihpbaoWfwqNU3TcqssNXczU3Mrt2+ZaUaLWW7gLu6KG4sgwiDbwMz5/WFOoqAMAXeW9+t55lHunHvnM0dg3t577jkyIYQAERERkUTkUhdARERE1o1hhIiIiCTFMEJERESSYhghIiIiSTGMEBERkaQYRoiIiEhSDCNEREQkKYYRIiIikhTDCBEREUmKYYSIyggMDMQzzzwjdRlEZEUYRohMxNdffw2ZTAaZTIZdu3bd87wQAn5+fpDJZCYfFr7//nuEhoZCrVbD398fEydORGlp6QP3S01NRf/+/dGsWTM4OTnB1dUVERERWL58Oe5eueKDDz4w9NedD7Vafd/X2LVrl6FtVlZWlY+5YMEC9OrVC/7+/pDJZBg0aFCFr7ljxw489NBDsLe3h5ubG3r27ImUlJQH9geRtbCRugAiKkutVmP16tV46KGHymz/448/cOXKFahUKokqq5yffvoJ3bt3xyOPPIIvvvgCR48exUcffYRr165hwYIF9903KysLV65cQc+ePeHv74+SkhLs2LEDgwYNwqlTpzB16tR79lmwYAEcHR0NXysUigqPr9frMXLkSDg4OCA/P7/CdpU55owZM5CXl4eIiAikpaVVeKwtW7bgueeeQ2hoKKZPnw6NRoPZs2fjoYcewsGDB+Hp6VnhvkRWQxCRSVi2bJkAIJ5//nnh4eEhSkpKyjw/dOhQERYWJgICAkS3bt1qrI7/evyWLVuKoKCgMvWPGzdOyGQyceLEiSod85lnnhEODg6itLTUsG3ixIkCgMjMzKz0cRYsWCDc3d3FqFGjyt3XmGOmpKQIvV4vhBDCwcFBxMXFlduuZcuWonHjxqK4uNiw7dChQ0Iul4v4+PhK105kyXiZhsjExMbG4vr169ixY4dhm1arxfr169G3b99y9/n000/RoUMHuLu7w87ODmFhYVi/fn25bb/55htEREQYLhk8/PDD+Pnnn+9pt2vXLkRERECtVqNhw4ZYsWLFA2tPTk5GcnIyhg0bBhubf0+8vv766xBCVFjTgwQGBqKgoABarfae54QQ0Gg091zGuVt2djbGjx+PyZMnw9XV9b5tK3PMgIAAyGSyB75mcnIyevToAaVSadgeFBSEFi1aYO3atffdn8haMIwQmZjAwEBERUVhzZo1hm0//fQTcnNz0adPn3L3mT17NkJCQjB58mRMnToVNjY26NWrF3788ccy7SZNmoQBAwbA1tYWkydPxqRJk+Dn54dff/21TLuzZ8+iZ8+eeOKJJzBz5ky4ublh0KBBOH78+H1rP3jwIAAgPDy8zHZfX1/Ur1/f8PyDFBYWIisrCykpKVi+fDmWLVuGqKgo2NnZ3dO2YcOGcHFxgZOTE/r374+MjIxyj/n+++/Dx8cHr7zyygNfv7LHfJDi4mIAKLdue3t7pKamIj09vUrHJrIkHDNCZIL69u2LsWPHorCwEHZ2dli1ahU6d+4MX1/fctufPn26zAfeiBEjEBoais8++wzdunUDcCtgTJ48GT169MD69eshl//7f5G7zwCcOnUKO3fuRKdOnQAAL774Ivz8/LBs2TJ8+umnFdZ9e+xE3bp173mubt26SE1NrdT7nz17NsaOHWv4+vHHH8eyZcvKtHFzc8OIESMQFRUFlUqFP//8E/PmzcO+ffuQmJgIZ2dnQ9sjR45g0aJF2Lp1633HlBhzzMrw9vaGq6sr/vrrrzLbr1+/juTkZADA1atX4ePjY9RxiSwNwwiRCXrxxRfx5ptvYsuWLejatSu2bNmCOXPmVNj+ziBy48YN6HQ6dOrUqczZlc2bN0Ov12PChAllggiAey43tGzZ0hBEAMDT0xPNmjXD+fPn71t3YWEhAJQ7yFatVkOj0dx3/9tiY2MRHh6OzMxMbNmyBRkZGYZj3zZq1KgyX7/wwguIiIhAv379MH/+fIwZM8bw3BtvvIGnnnoKTz755H1f15hjVoZcLscrr7yCGTNmYOzYsXjppZeg0Wjw7rvvGi453f2+iKwRL9MQmSBPT09ER0dj9erV2LhxI3Q6HXr27Flh+y1btqB9+/ZQq9WoU6cOPD09sWDBAuTm5hranDt3DnK5HC1btnzg6/v7+9+zzc3NDTdu3LjvfrdD0e3LE3cqKioq93JFeQICAhAdHY3Y2FisWrUKDRs2RHR09AM/uPv27QsfHx/88ssvhm3r1q3D33//jZkzZ1bqtStzTGNMnjwZL7/8Mj7++GM0bdoU4eHhsLGxwcsvvwwAZe7aIbJWDCNEJqpv37746aefsHDhQjz11FMVDrr8888/8eyzz0KtVmP+/PnYunUrduzYgb59+z5wUGdFKrqU8aDj3b48U96trmlpaRVeZnqQnj174vLly9i5c+cD2/r5+SE7O9vw9TvvvINevXpBqVQiJSUFKSkpyMnJAQBcvny5UpeO7j6mMZRKJb788kukpqZi586dOHXqFLZv347c3FzI5XI0bty4SsclsiQMI0QmqkePHpDL5dizZ0+Fd9EAwIYNG6BWq7F9+3a89NJLeOqppxAdHX1Pu0aNGkGv1xvGKtSE4OBgAEBiYmKZ7ampqbhy5YrheWPdPiNy55me8gghkJKSUmbujsuXL2P16tVo0KCB4TF79mwAQGhoKJ5++mmjj1kV3t7e6NSpE5o2bQqdTofff/8dkZGRPDNCBI4ZITJZjo6OWLBgAVJSUhATE1NhO4VCAZlMBp1OZ9iWkpKCzZs3l2nXvXt3jB49GpMnTy53AOuDblOtjFatWqF58+ZYvHgxXnnlFcMZlgULFkAmk5W51JSbm4u0tDTUrVsXLi4uAIDMzMxyP/SXLl0KmUyG0NBQw7by2i5YsACZmZno2rWrYdumTZvuOd7atWuxbt06rFixAvXr1zf6mP/Vp59+irS0NHzxxRfVdkwic8YwQmTC4uLiHtimW7du+Oyzz9C1a1f07dsX165dw7x589C4cWMcOXLE0K5x48YYN24cPvzwQ3Tq1AnPP/88VCoV9u/fD19fX0ybNq1aav7kk0/w7LPP4sknn0SfPn1w7NgxzJ07F0OGDEGLFi0M7TZt2oTBgwdj2bJlhqnUp0yZgr/++gtdu3aFv78/srOzsWHDBuzfvx8jR44sc0kjICAAvXv3Rps2baBWq7Fr1y6sXbsWwcHBZW7f7d69+z01Hjp0CADw1FNPwcPDw+hjAsAPP/yAw4cPAwBKSkpw5MgRfPTRRwCAZ599Fm3btgVwa16XDRs24OGHH4ajoyN++eUXfPvttxgyZAheeOGFqnc0kSWRaLI1IrrL7RlY9+/ff9925c2QunTpUtGkSROhUqlE8+bNxbJlywyzid7tq6++EiEhIUKlUgk3NzfRuXNnsWPHjvseXwghOnfuLDp37lyp97Jp0yYRHBwsVCqVqF+/vhg/frzQarXlvt9ly5YZtv3888/imWeeEb6+vsLW1lY4OTmJjh07imXLlhlmO71tyJAhomXLlsLJyUnY2tqKxo0bi9GjRwuNRvPA+iqaadWYY8bFxQkA5T7ufE979+4VDz/8sHBzcxNqtVoEBQWJhQsX3vN+iKyZTIgqjnAjIiIiqgYcwEpERESSYhghIiIiSTGMEBERkaQYRoiIiEhSDCNEREQkKYYRIiIikpRZTHqm1+uRmpoKJyenapklkoiIiGqeEAJ5eXnw9fW9Z7XwO5lFGElNTYWfn5/UZRAREVEVXL58uczSC3czizDi5OQE4NabcXZ2lrgaIiIiqgyNRgM/Pz/D53hFzCKM3L404+zszDBCRERkZh40xIIDWImIiEhSDCNEREQkKYYRIiIikpRZjBmpDL1eD61WK3UZVEOUSuV9bwsjIiLzZRFhRKvV4sKFC9Dr9VKXQjVELpejQYMGUCqVUpdCRETVzOzDiBACaWlpUCgU8PPz4/+eLdDtSe/S0tLg7+/Pie+IiCyM2YeR0tJSFBQUwNfXF/b29lKXQzXE09MTqampKC0tha2trdTlEBFRNTL70wg6nQ4AePrewt3+9739701ERJbD7MPIbTx1b9n470tEZLksJowQERGReWIYMVG///47ZDIZcnJypC7lvhISEtCiRYv/dPlk27ZtCA4O5t1QRERWimFEAjKZ7L6PDz74QNL6AgMDMWvWrEq1fffddzF+/HgoFAoAwMGDBxESEgJHR0fExMQgOzvb0La0tBRhYWHYt29fmWN07doVtra2WLVqVbW9ByIiMh8MIxJIS0szPGbNmgVnZ+cy295+++0qHbe2J33btWsXzp07hxdeeMGwbciQIXjsscdw4MAB5ObmYurUqYbnZs6ciY4dOyIiIuKeYw0aNAhz5syplbqJiOhfRSU6HL6cA71eSFYDw4gEfHx8DA8XFxfIZLIy2xwdHQ1tk5KSEB4eDnt7e3To0AGnTp0yPPfBBx8gODgYX375JRo0aAC1Wg2g/DMbwcHBhjMuQgh88MEH8Pf3h0qlgq+vL9544w0AwCOPPIKLFy/irbfeMpypqcjatWvxxBNPGF4XAE6cOIGhQ4eiadOmiI2NxYkTJwAA58+fx9KlSzFlypRyjxUTE4PExEScO3eu8h1JRERGKSrR4eClG1i55yJGrz+Cp2f/idYTt+O5eX/hfFa+ZHWZ/TwjdxNCoLBEmts/7WwV1X7Xx7hx4zBz5kx4enri1VdfxUsvvYS//vrL8PzZs2exYcMGbNy40XCp5EE2bNiAzz//HGvXrkWrVq2Qnp6Ow4cPAwA2btyIoKAgDBs2DEOHDr3vcf7880/07du3zLaLine truncated
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAHHCAYAAABtF1i4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABOIUlEQVR4nO3deVhUZf8G8HsWGPZFQRAFwX0HBEG0spLSMkorcxeXTMty+7VoblluZZlvuaW575XLW25lmJW9KgLu4IqIG5sIwz4w8/z+UCdHFkGRM8v9ua65hOc858x3Dg5zc5bnkQkhBIiIiIgkIpe6ACIiIrJsDCNEREQkKYYRIiIikhTDCBEREUmKYYSIiIgkxTBCREREkmIYISIiIkkxjBAREZGkGEaIiIhIUgwjRGTA19cXL730ktRlEJEFYRghMhKrVq2CTCaDTCbDgQMHSi0XQsDb2xsymczow8LPP/+Mdu3awcbGBj4+Ppg2bRpKSkqqvJ3169dDJpPBwcHBoF2n02HVqlV4+eWX4e3tDXt7e7Ru3RozZsxAYWFhqe1kZ2fjww8/RJMmTWBra4sGDRpg2LBhSE5ONui3bds2dO3aFV5eXlCpVKhfvz5ef/11nDp1yqDf/v379T+rsh4zZ87U93366afL7WdlZVXlfUJkjpRSF0BEhmxsbLBhwwY88cQTBu1//vknrl69CpVKJVFllbN792706NEDTz/9NL799lucPHkSM2bMQFpaGhYvXlzp7eTm5uLDDz+Evb19qWX5+fkYMmQIOnTogJEjR6JOnTo4ePAgpk2bhqioKOzbtw8ymQzA7eDy3HPPIT4+Hu+88w6aNm2KCxcuYNGiRfj111+RkJAAR0dHAMDJkyfh6uqKMWPGwM3NDSkpKVixYgVCQkJw8OBB+Pv7AwBatGiBtWvXlqpr7dq1+O233/D888/r2yZNmoQ333zToF9eXh5Gjhxp0I/IogkiMgorV64UAMSrr74q3NzcRHFxscHy4cOHi6CgINGgQQPRvXv3x1bHo26/ZcuWwt/f36D+SZMmCZlMJhISEiq9nY8++kg0a9ZM9O/fX9jb2xssKyoqEv/880+pdaZPny4AiL179+rb/vnnHwFALFiwwKDvihUrBACxdevWCutISUkRSqVSjBgx4oE1N27cWDRp0uSB/dauXSsAiPXr1z+wL5El4GkaIiPTt29f3Lx5E3v37tW3aTQa/PTTT+jXr1+Z63z55Zfo2LEjateuDVtbWwQFBeGnn34qs++6desQEhICOzs7uLq64qmnnsJvv/1Wqt+BAwcQEhICGxsbNGzYEGvWrHlg7fHx8YiPj8dbb70FpfLfA6/vvPMOhBDl1nS/8+fP4+uvv8a8efMMtnOXtbU1OnbsWKq9Z8+eAICEhAR9m1qtBgB4eHgY9K1bty4AwNbWtsJa6tSpAzs7O2RlZVXYLzo6GhcuXED//v0r7AcAGzZsgL29PV555ZUH9iWyBAwjREbG19cXYWFh2Lhxo75t9+7dyM7ORp8+fcpc5z//+Q8CAwPx6aefYtasWVAqlejVqxd27txp0G/69OkYOHAgrKys8Omnn2L69Onw9vbGvn37DPpduHABr7/+Op577jl89dVXcHV1xeDBg3H69OkKaz969CgAIDg42KDdy8sL9evX1y9/kLFjx+KZZ57Biy++WKn+d6WkpAAA3Nzc9G3BwcGwt7fHlClTsG/fPly7dg1//vknPvzwQ7Rv3x7h4eGltpOVlYX09HScPHkSb775JtRqNbp06VLhc69fvx4AHhhG0tPTsXfvXvTo0aPMU1BEFknqQzNEdNvd0zRHjhwRCxYsEI6OjiI/P18IIUSvXr3EM888I4Qo+zTK3X53aTQa0bp1a/Hss8/q286fPy/kcrno2bOn0Gq1Bv11Op3+6wYNGggA4q+//tK3paWlCZVKJf7v//6vwtcwd+5cAUAkJyeXWta+fXvRoUOHCtcXQogdO3YIpVIpTp8+LYQQIjIystRpmvKEh4cLJycncevWrVLbrFu3rgCgf3Tt2lXk5OSUuZ1mzZrp+zk4OIjJkyeX2mf3KikpER4eHiIkJOSBNX777bcCgNi1a1elXhORJeCRESIj9MYbb6CgoAA7duxATk4OduzYUe4pGsDwVMOtW7eQnZ2NJ598EnFxcfr27du3Q6fTYerUqZDLDd/6dy/2vKtly5Z48skn9d+7u7ujWbNmSExMrLDugoICACjzIlsbGxv98vJoNBqMGzcOI0eORMuWLSvse79Zs2bh999/x5w5c+Di4mKwzN3dHYGBgZg5cya2b9+OTz75BH///TeGDBlS5rZWrlyJPXv2YNGiRWjRogUKCgqg1WrLfe6oqCikpqZW+hSNu7s7nnvuuSq9PiJzxrtpiIyQu7s7wsPDsWHDBuTn50Or1eL1118vt/+OHTswY8YMHDt2DEVFRfr2e0PGxYsXIZfLK/Uh7+PjU6rN1dUVt27dqnC9u6Ho3hruKiwsfOD1GV9//TUyMjIwffr0B9Z4r82bN2Py5MkYNmwY3n77bYNliYmJeOaZZ7BmzRq89tprAIBXXnkFvr6+GDx4MHbv3o0XXnjBYJ2wsDD913369EGLFi0A3L42pyzr16+HQqFA7969K6wzMTERBw8exLvvvlvmtTBElopHRoiMVL9+/bB7924sWbIEL7zwQqm/9u/6+++/8fLLL8PGxgaLFi3Crl27sHfvXvTr1w9CiId6boVCUWb7g7Z396LQGzdulFp248YNeHl5lbtudnY2ZsyYgeHDh0OtViMpKQlJSUnIzc2FEAJJSUlIS0srtd7evXsxaNAgdO/eHUuWLCm1fNWqVSgsLCw1NsvLL78MAPjnn38qfE2urq549tln9deE3K+goADbtm1DeHh4qYtk77dhwwYAD76uhMjSMIwQGamePXtCLpfj0KFDFZ6i2bJlC2xsbPDrr79i6NCheOGFF8q8KLNRo0bQ6XSIj49/bDUHBAQAAGJiYgzar1+/jqtXr+qXl+XWrVvIzc3FF198AT8/P/1jy5YtyM/Ph5+fH9566y2DdQ4fPoyePXsiODgYP/zwQ5lHG1JTUyGEKHWapbi4GAAqNRhbQUEBsrOzy1z2888/Iycnp9KnaBo1aoQOHTo8sC+RJWEYITJSDg4OWLx4MT755BNERESU20+hUEAmkxl82CYlJWH79u0G/Xr06AG5XI5PP/0UOp3OYNnDHkG5X6tWrdC8eXMsXbrUoJ7FixdDJpMZnGrKzs7GmTNn9B/yderUwbZt20o9nnnmGdjY2GDbtm2YOHGifv2EhAR0794dvr6+2LFjR7mngJo2bQohBH744QeD9rt3KwUGBurbyjrykpSUhKioqFJ3CN21YcMG2NnZ6W8rLs/Ro0eRkJBQYbAkslQ8aUlkxCIjIx/Yp3v37pg3bx66deuGfv36IS0tDQsXLkTjxo1x4sQJfb/GjRtj0qRJ+Oyzz/Dkk0/i1VdfhUqlwpEjR+Dl5YXZs2dXS81z587Fyy+/jOeffx59+vTBqVOnsGDBArz55pv6ay+A20OvDxkyBCtXrsTgwYNhZ2eHHj16lNre9u3bER0dbbAsJycHXbt2xa1bt/DBBx+UuoW5UaNG+us+Bg8ejC+//BIjRozA0aNH0apVK8TFxeH7779Hq1atDEJEmzZt0KVLFwQEBMDV1RXnz5/H8uXLUVxcjDlz5pSqLTMzE7t378Zrr71Wasj6+1X21l8iiyTlrTxE9K97b+2tSFm39i5fvlw0adJEqFQq0bx5c7Fy5Uoxbdo0UdZbfMWKFSIwMFCoVCrh6uoqOnfubDBiaXkjsHbu3Fl07ty5Uq9l27ZtIiAgQKhUKlG/fn0xefJkodFoyny9K1eurHBbZd3ae+nSJYPbdO9/REZGGvS/evWqGDp0qPDz8xPW1taibt26Yvjw4SI9Pd2g37Rp00RwcLBwdXUVSqVSeHl5iT59+ogTJ06UWduSJUsEAPHzzz9X+Bq0Wq2oV6+eaNeuXYX9iCyVTIhqOj5LRERE9BB4zQgRERFJimGEiIiIJMUwQkRERJJiGCEiIiJJMYwQERGRpBhGiIiISFImMeiZTqfD9evX4ejoWGp2USIiIjJOQgjk5OTAy8ur1Gzh9zKJMHL9+nV4e3tLXQYRERE9hCtXrqB+/frlLjeJMOLo6Ajg9otxcnKSuBoiIiKqDLVaDW9vb/3neHlMIozcPTXj5OTEMEJERGRiHnSJBS9gJSIiIkkxjBAREZGkGEaIiIhIUiZxzUhl6HQ6aDQaqcugx8Ta2rrC28KIiMh0mUUY0Wg0uHTpEnQ6ndSl0GMil8vh5+cHa2trqUshIqJqZvJhRAiBGzduQKFQwNvbm389m6G7g97duHEDPj4+HPiOiMjMmHwYKSkpQX5+Pry8vGBnZyd1OfSYuLu74/r16ygpKYGVlZXU5RARUTUy+cMIWq0WAHj43szd/fne/XkTEZH5MPkwchcP3Zs3/nyJiMyX2YQRIiIiMk0MI0Zq//79kMlkyMrKkrqUCkVFRaFFixaPdPpkz549CAgI4N1QREQWimFEAjKZrMLHJ598Iml9vr6+mD9/fqX6fvjhh5g8eTIUCgUA4OjRowgMDISDgwMiIiKQmZmp71tSUoKgoCBER0cbbKNbt26wsrLC+vXrq+01EBGR6WAYkcCNGzf0j/nz58PJycmg7f3333+o7db0oG8HDhzAxYsX8dprr+nb3nzzTTz77LOIi4tDdnY2Zs2apV/21VdfoVOnTggJCSm1rcGDB+Obb76pkbqJiOhfmhIdDl68KWkNDCMS8PT01D+cnZ0hk8kM2hwcHPR9Y2NjERwcDDs7O3Ts2BFnz57VL/vkk08QEBCA77//Hn5+frCxsQFQ9pGNgIAA/REXIQQ++eQT+Pj4QKVSwcvLC6NHjwYAPP3007h8+TLGjRunP1JTnk2bNuG5557TPy8AJCQkYPjw4WjatCnLine truncated
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAHHCAYAAABtF1i4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABPb0lEQVR4nO3deVxUVR8G8GdmYAaQTWQXEPddQBBCM1vINUzLcl+o9K3MJSvfyN1ULM3MPU1Nc01T8zWXjNS0TATFFXHDXUBE9mVg5rx/kJOjoKDgneX5fj7zKe6ce+c3B2Qe7j3nXJkQQoCIiIhIInKpCyAiIiLzxjBCREREkmIYISIiIkkxjBAREZGkGEaIiIhIUgwjREREJCmGESIiIpIUwwgRERFJimGEiIiIJMUwQkQ6vr6+eOWVV6Qug4jMDMMIkQH4/vvvIZPJIJPJcODAgQeeF0LA29sbMpnM4MPC1q1b0bJlS1hZWcHHxwcTJkxAcXHxI/e7dOmSrg/uf6xbt+6B9gkJCejYsSNsbW3h5OSE/v3749atW3ptzpw5g9GjR8Pf3x92dnbw8PBAly5dEBsb+8DxNm3ahJ49e6JOnTqwsbFBw4YN8dFHHyEjI+OBth9++CFatmwJJycn2NjYoHHjxpg4cSJycnL02u3du7fM9/T3338/sk+IzIWF1AUQ0b+srKywZs0aPPvss3rb9+3bh2vXrkGlUklUWfns2LED3bp1w/PPP4+5c+fixIkTmDJlClJTU7Fw4cJyHaN3797o3Lmz3rbQ0FC9r69du4bnnnsODg4OmDZtGnJycjBz5kycOHECMTExUCqVAIDvvvsOS5cuxeuvv473338fmZmZ+Pbbb/HMM89g586dCAsL0x1zyJAh8PT0RL9+/eDj44MTJ05g3rx52L59O44cOQJra2td28OHD6Nt27aIiIiAlZUVjh49iunTp+O3337DH3/8Ablc/++84cOHo1WrVnrb6tWrV67+IDILgogkt3z5cgFAvPbaa8LZ2VkUFRXpPT948GARGBgoatWqJbp06VJldTzp8Zs0aSL8/Pz06h8zZoyQyWQiISHhofsmJSUJAGLGjBmPfJ333ntPWFtbi8uXL+u27d69WwAQ3377rW5bbGysyM7O1ts3LS1NuLi4iDZt2uht37NnzwOvs2LFCgFALFmy5JE1zZw5UwAQBw8e1DsmALFhw4ZH7k9kzniZhsiA9O7dG7dv38bu3bt129RqNTZu3Ig+ffqUus/MmTPRunVr1KhRA9bW1ggMDMTGjRtLbbtq1SoEBwfDxsYG1atXx3PPPYdff/31gXYHDhxAcHAwrKysUKdOHaxcufKRtZ8+fRqnT5/GkCFDYGHx70nX999/H0KIMmsqTW5uLtRqdZnP//TTT3jllVfg4+Oj2xYWFoYGDRrgxx9/1G0LDAyEra2t3r41atRA27ZtkZCQoLf9+eeff+B1unfvDgAPtC2Nr68vAJR6WQcAsrOzy3W5isgcMYwQGRBfX1+EhoZi7dq1um07duxAZmYmevXqVeo+33zzDQICAjB58mRMmzYNFhYWeOONN/DLL7/otZs0aRL69+8PS0tLTJ48GZMmTYK3tzd+//13vXbnz59Hjx498PLLL+Orr75C9erVMWjQIJw6deqhtR89ehQAEBQUpLfd09MTXl5euucfZdKkSbC1tYWVlRVatWr1QFi6fv06UlNTH3gdAAgODi7X6yQnJ8PZ2blc7QCU2ra4uBhpaWm4ceMGfv31V4wdOxZ2dnYIDg5+oG1ERATs7e1hZWWFF154odQxK0TmjGNGiAxMnz59EBkZifz8fFhbW2P16tVo164dPD09S21/9uxZvfEMH3zwAVq2bIlZs2ahS5cuAEoCxuTJk9G9e3ds3LhRb0yDEELveImJifjjjz/Qtm1bAMCbb74Jb29vLF++HDNnziyz7ps3bwIAPDw8HnjOw8MDN27ceOj7lsvlaN++Pbp3746aNWvi4sWLmDVrFjp16oStW7fq3sujXic9PR2FhYVljq/Zv38/Dh48iLFjxz60HgD44osvoFAo0KNHjweei42N1RvL0rBhQ2zduhVOTk66bUqlEq+//jo6d+4MZ2dnnD59GjNnzkTbtm3x119/ISAg4JE1EJkDhhEiA/Pmm29i5MiR2LZtGzp27Iht27Zhzpw5Zba/N4jcuXMHGo0Gbdu21Tu7smXLFmi1WowfP/6BwZUymUzv6yZNmuiCCAC4uLigYcOGuHjx4kPrzs/PB4BSQ4CVlRWysrIeur+Pjw927dqlt61///5o0qQJPvroI10YedTr3G1T2vOpqano06cPateujdGjRz+0njVr1mDp0qUYPXo06tev/8DzTZo0we7du5Gbm4u//voLv/322wOzaVq3bo3WrVvrvu7atSt69OiBFi1aIDIyEjt37nxoDUTmgmGEyMC4uLggLCwMa9asQV5eHjQaTal/md+1bds2TJkyBfHx8SgsLNRtvzdkXLhwAXK5HE2aNHnk6987DuOu6tWr486dOw/d724oureGuwoKCvRCU3k5OTkhIiIC06dPx7Vr1+Dl5fXI17m3lnvl5ubilVdeQXZ2Ng4cOPDAWJJ77d+/H2+//TY6dOiAqVOnltrG3t5eNxvn1VdfxZo1a/Dqq6/iyJEj8PPzK/PY9erVw6uvvopNmzZBo9FAoVCU3QFEZoJjRogMUJ8+fbBjxw4sWrQInTp1gqOjY6nt9u/fj65du8LKygoLFizA9u3bsXv3bvTp0+eByy/lVdaH46OOd/eyyd3LKPe6efNmmZeZHsXb2xsAkJ6eXq7XcXJyeuCsiFqtxmuvvYbjx4/j559/RrNmzcp8vWPHjqFr165o1qwZNm7cqDcY92Fee+01ACh1TZTS3pNarUZubm65jk1k6hhGiAxQ9+7dIZfL8ffff5c5iwYomVViZWWFXbt24a233kKnTp301s64q27dutBqtTh9+nSV1ezv7w8ADwzOvHHjBq5du6Z7vqLuXh5ycXEBANSsWRMuLi6lDgKNiYl54HW0Wi0GDBiA6OhorFmzBu3atSvztS5cuICOHTvC1dUV27dvf+jZk/sVFhZCq9UiMzOzXO/JysqqQscnMmUMI0QGyNbWFgsXLsTEiRMRHh5eZjuFQgGZTAaNRqPbdunSJWzZskWvXbdu3SCXyzF58mRotVq95x73DMr9mjZtikaNGmHx4sV69SxcuBAymUzvUlNmZibOnDmj98F9/+qpQMnMmWXLlqFFixZ6A1Zff/11bNu2DVevXtVti46OxtmzZ/HGG2/oHWPYsGFYv349FixYoDt7UZrk5GS0b98ecrkcu3bt0oWf+2VkZKCoqOiB7d999x0A/dlEpb2nY8eOYevWrbrXIiKOGSEyWAMHDnxkmy5dumDWrFno2LEj+vTpg9TUVMyfPx/16tXD8ePHde3q1auHMWPG4PPPP0fbtm3x2muvQaVS4fDhw/D09ERUVFSl1Dxjxgx07doV7du3R69evXDy5EnMmzcP77zzDho3bqxrt3nzZkRERGD58uUYNGgQAGD06NG4cOECXnrpJXh6euLSpUv49ttvkZubi2+++UbvdT777DNs2LABL7zwAkaMGIGcnBzMmDEDzZs3R0REhK7d7NmzsWDBAoSGhsLGxgarVq3SO0737t1RrVo1AEDHjh1x8eJFjB49GgcOHNBblt/NzQ0vv/wygJIl3ocPH44ePXqgfv36UKvV2L9/PzZt2oSgoCD069dPt1/Pnj1hbW2N1q1bw9XVFadPn8bixYthY2OD6dOnV0qfE5kESZdcIyIhxL8rsB4+fPih7UpbIXXp0qWifv36QqVSiUaNGonly5eLCRMmiNL+eS9btkwEBAQIlUolqlevLtq1ayd279790OMLIUS7du1Eu3btyvVeNm/eLPz9/YVKpRJeXl5i7NixQq1Wl/p+ly9frtu2Zs0a8dxzzwkXFxdhYWEhnJ2dRffu3UVcXFypr3Py5EnRvn17YWNjIxwdHUXfvn1FcnKyXpuBAwcKAGU+kpKSdG0f1u7e937+/HkxYMAAUadOHWFtbS2srKxE06ZNxYQJE0ROTo7e63/zzTciODhYODk5CQsLC+Hh4SH69esnzp07V66+JDIXMiEq6RwtERER0WPgBUsiIiKSFMMIERERSYphhIiIiCTFMEJERESSYhghIiIiSTGMEBERkaSMYtEzrVaLGzduwM7O7oE7jBIREZFhEkIgOzsbnp6eD11x2CjCyI0bN3Q3yyIiIiLjcvXqVXh5eZX5vFGEETs7OwAlb8be3l7iaoiIiKg8srKy4O3trfscL4tRhJG7l2bs7e0ZRoiIiIzMo4ZYcAArERERSYphhIiIiCTFMEJERESSMooxI+Wh1WqhVqulLoOqiFKpfOi0MCIiMl4mEUbUajWSkpKg1WqlLoWqiFwuR+3ataFUKqUuhYiIKpnRhxEhBG7evAmFQgFvb2/+9WyC7i56d/PmTfj4+HDhOyIiE2P0YaS4uBh5eXnw9PSEjY2N1OVQFXFxccGNGzdQXFwMS0tLqcshIqJKZPSnETQaDQDw9L2Ju/v9vfv9JiIi02H0YeQunro3bfz+EhGZLpMJI0RERGScGEYM1N69eyGTyZCRkSF1KQ8VHR2Nxo0bP9Hlk507d8Lf35+zoYiIzBTDiARkMtlDHxMnTpS0Pl9fX8yePbtcbUePHo2xY8dCoVAAAI4ePYqAgADY2toiPDwc6enpurbFxcUIDAxETEyM3jE6duwIS0tLrF69utLeAxERGQ+GEQncvHlT95g9ezbs7e31tn388cePddynvejbgQMHcOHCBbz++uu6be+88w5efPFFHDlyBJmZmZg2bZruua+++gpt2rRBcHDwA8caNGgQ5syZ81TqJiKifxUWa/Dn+TRJa2AYkYC7u7vu4eDgAJlMprfN1tZW1zYuLg5BQUGwsbFB69atkZiYqHtu4sSJ8PfLine truncated
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAHHCAYAAABtF1i4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABO2klEQVR4nO3deVhUZcMG8HsGmGHfZBMEwX0XZAstbaHU1NIWN3Jf3nLJok2yNC3FsnzNvSwptzRz+crSMlzSNxVBcQV3xIVFRBj2gZnn+8OcHAUFBc4s9++6uJQz58zcc1C4Oec5z5EJIQSIiIiIJCKXOgARERGZN5YRIiIikhTLCBEREUmKZYSIiIgkxTJCREREkmIZISIiIkmxjBAREZGkWEaIiIhIUiwjREREJCmWESLS4+/vj969e0sdg4jMCMsIkYH47rvvIJPJIJPJsHfv3rseF0LA19cXMpnM4MvCzz//jE6dOsHa2hp+fn6YNm0aKioq7rtdWlqabh/c+bF27Vq9dataTyaT4emnn77ruc+dO4fBgwfDw8MDNjY2aN68OaZMmXLXeikpKejRowfs7e3h6uqKIUOG4Nq1a3etd/bsWbz00ktwcXGBra0tHn30UezcufOu9ZYtW4Zu3brB09MTSqUSAQEBGDFiBNLS0u67P4jMhaXUAYhIn7W1NdasWYNHH31Ub/nu3btx+fJlKJVKiZJVz9atW9G3b188/vjjWLBgAY4dO4ZPPvkE2dnZWLJkSbWeY9CgQXj22Wf1lkVEROh9vnLlyru2S0xMxJdffolnnnlGb3lycjIef/xx+Pj44K233kKDBg2Qnp6OS5cu6a13+fJldO3aFU5OTpg1axYKCwvx+eef49ixY0hISIBCoQAAXLp0CREREbCwsMA777wDOzs7xMXF4ZlnnkF8fDy6du2qe87Dhw8jICAAzz33HFxcXHDhwgUsW7YMW7ZswZEjR+Dt7V2tfUJk0gQRGYS4uDgBQLzwwgvCzc1NlJeX6z0+ZswYERwcLBo3bix69epVZzke9vnbtGkjOnbsqJd/ypQpQiaTiZSUlHtue+HCBQFAzJkz54Fee9SoUUImk4lLly7plmk0GtGuXTsRHh4uiouL77n9a6+9JmxsbMTFixd1y7Zv3y4AiK+++kq3bNy4ccLS0lKkpqbqlhUVFQlfX1/RqVOn++ZMTEwUAERsbGxN3h6RyeJpGiIDM2jQIFy/fh3bt2/XLVOr1fjpp58wePDgSrf5/PPP0blzZzRo0AA2NjYIDg7GTz/9VOm6q1atQlhYGGxtbeHi4oKuXbvijz/+uGu9vXv3IiwsDNbW1mjSpAlWrFhx3+wnT57EyZMnMXbsWFha/nvgddy4cRBCVJmpMkVFRVCr1dVev6ysDBs2bEC3bt3QqFEj3fI//vgDx48fx7Rp02BjY4Pi4mJoNJpKn2PDhg3o3bs3/Pz8dMsiIyPRokUL/Pjjj7ple/bsQVBQEFq2bKlbZmtri+eeew6HDh3CmTNn7pnV398fAJCXl1ft90dkylhGiAyMv78/IiIi8MMPP+iWbd26Ffn5+Rg4cGCl23z55ZcICgrCjBkzMGvWLFhaWuLll1/Gr7/+qrfe9OnTMWTIEFhZWWHGjBmYPn06fH19sWPHDr31bo2HePrpp/HFF1/AxcUFw4cPx4kTJ+6Z/fDhwwCAkJAQveXe3t5o1KiR7vH7mT59Ouzt7WFtbY3Q0NBKy9KdfvvtN+Tl5SEqKkpv+Z9//gkAUCqVCAkJgZ2dHWxtbTFw4EDk5ubq1rty5Qqys7Pvyg4AYWFhetnLyspgY2Nz13q2trYAgKSkpLseu379OrKzs5GYmIgRI0YAAJ566qn7vi8ic8AxI0QGaPDgwYiJiUFJSQlsbGywevVqdOvWrcrxBadPn9b74ThhwgR06tQJc+fORa9evQDcLBgzZsxAv3798NNPP0Eu//d3ESGE3vOdOnUKf/31Fx577DEAQP/+/eHr64u4uDh8/vnnVebOyMgAADRs2PCuxxo2bIirV6/e833L5XI888wz6NevH3x8fHD+/HnMnTsXPXv2xM8//6x7L5VZvXo1lEolXnrpJb3lt45S9O/fHz169EBMTAyOHDmC2NhYXLp0CXv37oVMJrtv9tzcXJSVlUGpVKJly5bYs2cPCgoK4ODgoFvv1sDjK1eu3PUcPj4+KCsrAwA0aNAA8+fPr3SgLZE5YhkhMkD9+/fHG2+8gS1btqBHjx7YsmUL5s+fX+X6txeRGzduQKPR4LHHHtM7urJ582ZotVpMnTpVr4gAN69MuV2bNm10RQQA3N3d0bJlS5w/f/6euUtKSgCg0kG21tbWUKlU99zez88Pv//+u96yIUOGoE2bNnjrrbeqLCMqlQq//vornn32WTg7O+s9VlhYCAAIDQ3FqlWrAAAvvvgibG1tERMTg/j4eERGRt43+633p1Qq8dprr+GXX37BgAEDMHPmTNjZ2WHx4sVITEzU2w+327p1K0pLS5GSkoJVq1ahqKjonvuCyJzwNA2RAXJ3d0dkZCTWrFmDjRs3QqPR3PUb/+22bNmCRx55BNbW1nB1dYW7uzuWLFmC/Px83Trnzp2DXC5HmzZt7vv6t4+ZuMXFxQU3bty453a3StGtIwC3Ky0trfTUxv24urpixIgROHXqFC5fvlzpOhs2bEBpaeldp2huzzRo0CC95bfG3/z999/Vyn77Oj179sSCBQvw119/oVOnTmjZsiV+/fVXzJw5EwBgb29/13M88cQT6NmzJ6Kjo7F+/XpMnz4dCxcuvP8OIDIDLCNEBmrw4MHYunUrli5dip49e971G/8te/bswXPPPQdra2ssXrwYv/32G7Zv347BgwffdfqluiwsLCpdfr/nu3WK49Ypj9tlZGQ88GWsvr6+AKA3xuN2q1evhpOTU6Xzr9x6TU9PT73lHh4eAKArWPfL7urqqnfUZMKECcjKysLff/+NxMREpKamwsnJCQDQokWLe76fpk2bIigoCKtXr77nekTmgmWEyED169cPcrkc+/fvr/IqGuDmUQFra2v8/vvvGDlyJHr27InIyMi71mvatCm0Wi1OnjxZZ5kDAwMBQHe64parV6/i8uXLusdr6tbpIXd397sey8jIwM6dO/Hiiy9WeoolODgYwN3jOG6NX7n1nD4+PnB3d78rOwAkJCRUmt3Ozg4REREIDg6GhYUF/vzzT9jY2KBLly73fU8lJSV6R66IzBnLCJGBsre3x5IlS/DRRx+hT58+Va5nYWEBmUymd7lqWloaNm/erLde3759IZfLMWPGDGi1Wr3HHvQIyp3atm2LVq1a4euvv9bLs2TJEshkMr1TTfn5+UhNTdX7gVzZTKdXrlzB8uXL0aFDh0oHl65duxZarbbSUzQA8Pzzz0OpVCIuLk7vfX/zzTcAoDeI9MUXX8SWLVv0JkOLj4/H6dOn8fLLL9/zvf/999/YuHEjRo0apTtCUlFRUemprYSEBBw7dqzSK3eIzBEHsBIZsGHDht13nV69emHu3Lno0aMHBg8ejOzsbCxatAjNmjXD0aNHdes1a9YMU6ZMwccff4zHHnsML7zwApRKJQ4ePAhvb2/ExsbWSuY5c+bgueeewzPPPIOBAwfi+PHjWLhwIUaPHo3WrVvr1tu0aRNGjBiBuLg4DB8+HADw7rvv4ty5c3jqqafg7e2NtLQ0fPXVVygqKsKXX35Z6eutXr0a3t7eePzxxyt93MvLC1OmTMHUqVPRo0cP9O3bF0eOHMGyZcswaNAghIaG6tZ9//33sX79ejzxxBOYNGkSCgsLMWfOHLRv3153OS4AXLx4Ef3798dzzz0HLy8vnDhxAkuXLkWHDh0wa9Ys3XqFhYXw9fXFgAED0LZtW9jZ2eHYsWOIi4uDk5MTPvzww4fY00QmRNIp14hI59YMrAcPHrznepXNkPrtt9+K5s2bC6VSKVq1aiXi4uLEtGnTRGX/xZcvXy6CgoKEUqkULi4uolu3bmL79u33fH4hhOjWrZvo1q1btd7Lpk2bRGBgoFAqlaJRo0bigw8+EGq1utL3GxcXp1u2Zs0a0bVrV+Hu7i4sLS2Fm5ub6Nevn0hKSqr0dVJTUwUAER0dfc88Wq1WLFiwQLRo0UJYWVkJX1/fSjMJIcTx48fFM888I2xtbYWzs7OIiooSmZmZeuvk5uaK559/Xnh5eQmFQiECAgLEe++9J1Qqld56ZWVlYtKkSaJDhw7C0dFRWFlZicaNG4tRo0aJCxcu3DMzkTmRCVFLx2eJiIiIHgDHjBAREZGkWEaIiIhIUiwjREREJCmWESIiIpIUywgRERFJimWEiIiIJGUUk55ptVpcvXoVDg4Od91dlIiIiAyTEAIFBQXw9va+627htzOKMnL16lXdjbKIiIjIuFy6dAmNGjWq8nGjKCMODg4Abr4ZR0dHidMQERFRdahUKvj6+up+jlfFKMrIrVMzjo6OLCNERERG5n5DLDiAlYiIiCTFMkJERESSYhkhIiIiSRnFmJHq0Gq1UKvVUsegOqJQKO55WRgRERkvkygjarUaFy5cgFarlToK1RG5XI6AgAAoFAqpoxARUS0z+jIihEBGRgYsLCzg6+vL355N0K1J7zIyMuDn58eJ74iITIzRl5GKigoUFxfD29sbtra2UsehOuLu7o6rV6+ioqICVlZWUschIqJaZPSHETQaDQDw8L2Ju/X1vfX1JiIi02H0ZeQWHro3bfz6EhGZLpMpI0RERGScWEYM1K5duyCTyZCXlyd1lHuKj49H69atH+r0ybZt2xAYGMiroYiIzBTLiARkMtk9Pz766CNJ8/n7+2PevHnVWvfdd9/FBx98AAsLCwDA4cOHERQUBHt7e/Tp0we5ubm6dSsqKhAcHIyEhAS95+jRowesrKywevXqWnsPRERkPFhGJJCRkaH7mDdvHhwdHfWWvf322w/0vPU96dvevXtx7tw5vPjii7plo0eLine truncated
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAHHCAYAAABtF1i4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABWv0lEQVR4nO3deVxUVf8H8M/MAMO+iYAgCG4oLqyCaGYLpWWWrbgjpT25pfE8lWRpWoql+bPcd829RW23jLQ0UQTEXVwRXNhEGARhYOb8/iBHJ0EZBO4An/frdV/EmXPvfOeY8PHec8+VCSEEiIiIiCQil7oAIiIiatoYRoiIiEhSDCNEREQkKYYRIiIikhTDCBEREUmKYYSIiIgkxTBCREREkmIYISIiIkkxjBAREZGkGEaISI+XlxeeeeYZqcsgoiaEYYTISKxZswYymQwymQx79+6963UhBDw8PCCTyYw+LHz//fcIDAyEubk5PD09MXXqVJSXl1d7/3PnzmHw4MFwdnaGhYUF2rVrh8mTJ1fZv6ysDL6+vpDJZJgzZ85dr589exYvvfQSHBwcYGlpiYceegi7du2q9FgLFixAx44doVQq4e7ujujoaBQVFen1SUtL0/1Z/XvbvHmzXt+EhASMGTMGQUFBMDU1hUwmq/Y4EDUVJlIXQET6zM3NsXHjRjz00EN67X/++ScuXboEpVIpUWXV88svv2DAgAF45JFHMH/+fBw9ehQff/wxsrOzsXjx4vvun5KSgkceeQTu7u7473//i2bNmiE9PR0ZGRlV7jN//nykp6dX+lpGRgbCwsKgUCjw9ttvw8rKCqtXr8aTTz6JuLg4PPzww7q+7777Lj799FO89NJLmDBhAk6cOIH58+fj+PHj+PXXX+869qBBg/D000/rtYWFhel9//PPP2PFihXo2rUrWrdujdOnT993DIiaHEFERmH16tUCgHjhhReEk5OTKCsr03t91KhRIigoSLRq1Ur069evzup40OP7+voKPz8/vfonT54sZDKZOHny5D331Wg0onPnziI0NFQUFxdX6/2ysrKEnZ2dmD59ugAgZs+erff6mDFjhImJiTh16pSuraioSHh4eIjAwEBd25UrV4SJiYkYNmyY3v7z588XAMT333+va7tw4UKl71WZzMxM3WcZO3as4I9dorvxMg2RkRk0aBCuXbuGnTt36trUajW++eYbDB48uNJ95syZgx49eqBZs2awsLBAUFAQvvnmm0r7rl+/HiEhIbC0tISDgwMefvhh/Pbbb3f127t3L0JCQmBubo7WrVvjyy+/vG/tJ06cwIkTJ/D666/DxOT2idcxY8ZACFFlTbf89ttvOHbsGKZOnQoLCwsUFxdDo9Hcc59JkybBx8cHQ4cOrfT1PXv2ICAgAD4+Pro2S0tLPPvss0hOTsaZM2cAAPHx8SgvL8fAgQP19r/1/b8vv9xSVFQEtVpdZX0uLi6wsLC452cgauoYRoiMjJeXF8LCwrBp0yZd2y+//IKCgoK7flHe8vnnnyMgIADTp0/HzJkzYWJigpdffhk//fSTXr9p06Zh2LBhMDU1xfTp0zFt2jR4eHjgjz/+0Ot3a47FE088gc8++wwODg4YMWIEjh8/fs/aDx06BAAIDg7Wa3dzc0PLli11r1fl999/BwAolUoEBwfDysoKlpaWGDhwIPLy8u7qn5CQgLVr12LevHlVzsUoLS2tNAxYWloCAJKSknT9ANzV99/97jRt2jRYW1vD3Nwc3bp1qzTUEdH9cc4IkREaPHgwYmJicPPmTVhYWGDDhg3o3bs33NzcKu1/+vRpvV+i48aNQ2BgIObOnYt+/foBqAgY06dPx/PPP49vvvkGcvntf4sIIfSOl5qair/++gu9evUCALzyyivw8PDA6tWrK50gesvVq1cBAC1atLjrtRYtWuDKlSv3/Ny3zlK88sor6Nu3L2JiYnD48GHExsYiIyMDe/fu1YUOIQTGjx+PiIgIhIWFIS0trdJj+vj4YM+ePSgsLISNjY2u/dYk4cuXL+v6AcDff/+NRx99VNdvz549ev0AQC6X48knn8Tzzz8Pd3d3nD9/HnPnzsVTTz2F77//XjfmRFQ9DCNERuiVV17BxIkT8eOPP6Jv37748ccf8cUXX1TZ/84gcv36dWg0GvTq1Uvv7Mr27duh1WoxZcoUvSAC4K6zCr6+vrogAgDNmzeHj48Pzp8/f8+6b968CQCVTrI1NzeHSqW65/43btwAAHTr1g3r168HALz44ouwtLRETEwM4uLiEB4eDqDi7qOjR4/e99LP6NGj8cMPPyAiIgIzZsyAlZUVFi1ahMTERL2aAwMDERoaik8++QTu7u549NFHcfLkSYwePRqmpqa6fgDg6el514TWYcOGwdfXF//9738ZRogMxMs0REaoefPmCA8Px8aNG7F161ZoNBq89NJLVfb/8ccf0b17d5ibm8PR0RHNmzfH4sWLUVBQoOtz7tw5yOVy+Pr63vf9PT0972pzcHDA9evX77nfrVB065LHnUpKSu47d+LW64MGDdJrvzVXZt++fQAAlUqFmJgYvP322/Dw8LjnMZ966inMnz8ff/31FwIDA+Hj44OffvoJM2bMAABYW1vr+n777bfw8/PDq6++Cm9vb/Tv3x+vvPIKAgIC9PpVxtHREVFRUUhNTcWlS5fu2ZeI9DGMEBmpwYMH45dffsGSJUvw1FNPwd7evtJ+e/bswbPPPgtzc3MsWrQIP//8M3bu3InBgwffdfmluhQKRaXt9zvercszty7X3Onq1atVXma65dbrLi4ueu3Ozs4AoAtDc+bMgVqtRkREBNLS0pCWlqYLANevX0daWprepNJx48YhKysL+/btQ2JiIk6dOgU7OzsAQPv27XX93N3dsXfvXpw+fRp//fUXLl26hE8//RQZGRl6/apyKxhVNr+FiKrGMEJkpJ5//nnI5XLs37+/yrtogIp/zZubm+PXX3/Fq6++iqeeekp3KeNObdq0gVarxYkTJ+qsZn9/fwDQXQK55cqVK7h06ZLu9aoEBQUB0J+fcWt/oOKMEQCkp6fj+vXr6NSpE7y9veHt7a27rDRz5kx4e3vf9TmtrKwQFhaGoKAgKBQK/P7777CwsEDPnj3vqqNdu3bo1asXXF1dceLECVy9erXSMf23W5exbtVJRNXDMEJkpKytrbF48WJ8+OGH6N+/f5X9FAoFZDKZ3i2waWlp2L59u16/AQMGQC6XY/r06dBqtXqv1fQMyr916tQJHTp0wLJly/TqWbx4MWQymd6lpoKCApw6dUrvUtJzzz0HpVKJ1atX69W4YsUKAMATTzwBAHjzzTexbds2vW3p0qUAgBEjRmDbtm3w9vauss59+/Zh69ateO2113RnSCqj1WrxzjvvwNLSEm+88YauPScn566+ly9fxqpVq9C1a9dKJ/ASUdU4gZXIiEVGRt63T79+/TB37lz07dsXgwcPRnZ2NhYuXIi2bdviyJEjun5t27bF5MmT8dFHH6FXr1544YUXoFQqcfDgQbi5uSE2NrZWap49ezaeffZZPPnkkxg4cCCOHTuGBQsWYOTIkejYsaOu37Zt2xAVFYXVq1djxIgRAABXV1dMnjwZU6ZMQd++fTFgwAAcPnwYy5cvx6BBg9CtWzcAFZNNAwMD9d731t00nTp1woABA3TtFy9exCuvvIJnn30Wrq6uOH78OJYsWYKuXbti5syZeseYMGECSkpK4O/vj7KyMmzcuFF3+/Cd82jeeecdnDt3Do8//jjc3NyQlpaGpUuXoqioCJ9//rneMS9evIh169YBuH3G6OOPPwYAtGrVCsOGDavhSBM1IlKuuEZEt91agfXgwYP37FfZCqkrV64U7dq1E0qlUnTo0EGsXr1aTJ06tdLVPletWiUCAgKEUqkUDg4Oonfv3mLnzp33PL4QQvTu3Vv07t27Wp9l27Ztwt/fXyiVStGyZUvx/vvvC7VaXennXb16tV67VqsV8+fPF+3btxempqbCw8Oj0v3/rapVUfPy8sRzzz0nXF1dhZmZmfD29hbvvvuuUKlUdx1j9erVws/PT1hZWQkbGxvx+OOPiz/++OOufhs3bhQPP/ywaN68uTAxMRFOTk7i+eefF0lJSXf13bVrlwBQ6Vbd8SRq7GRC1NL5WSIiIqIa4JwRIiIikhTDCBEREUmKYYSIiIgkxTBCREREkmIYISIiIkkxjBAREZGkGsSiZ1qtFleuXIGNjc1dTxclIiIi4ySEQGFhIdzc3O56WvidGkQYuXLlyn2fzElERETGKSMjAy1btqzy9QYRRmxsbABUfBhbW1uJqyEiIqLqUKlU8PDw0P0er0qDCCO3Ls3Y2toyjBARETUw95tiwQmsREREJCmGESIiIpIUwwgRERFJqkHMGakOrVYLtVotdRlUR8zMzO55WxgRETVcjSKMqNVqXLhwAVqtVupSqI7I5XJ4e3vDzMxM6lKIiKiWNfgwIoTA1atXoVAo4OHhwX89N0K3Fr27evUqPD09ufAdEVEj0+DDSHl5OYqLi+Hm5gZLS0upy6E60rx5c1y5cgXl5eUwNTWVuhwiIqpFDf40gkajAQCevm/kbv353vrzJiKixqPBh5FbeOq+ceOfLxFR49VowggRERE1TAwjRmr37t2QyWTIz8+XupR7iouLQ8eOHR/o8smOHTvg7+/Pu6GIiJqoGoWRhQsXwsvLC+bm5ggNDUVCQsI9+8+bNw8+Pj6wsLCAh4cH3nrrLZSUlNSo4MZAJpPdc/vwww8lrc/Lywvz5s2rVt933nkH77//PhQKBQDg0KFDCAgIgLW1Nfr374+8vDxd3/LycgQFBd31/0vfvn1hamqKDRs21NpnICKihsPgMLJlyxZER0dj6tSpSE5Ohp+fH/r06YPs7OxK+2/cuBGTJk3C1KlTcfLkSaxcuRJbtmzBe++998DFN1RXr17VbfPmzYOtra1e2//+978aHbeLine truncated
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAHHCAYAAABtF1i4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABVXklEQVR4nO3deVxUVf8H8M/MwAz7oigIgrjlLiBbaGkLPVhmqWnuIG5ZZha/FnezTOxxeSw1yd3cs9TKXDLKyiIREHdxF7dhEWGQbWDm/P4wR0lQBoE7A5/363VfxZlz73znuszHe889RyaEECAiIiKSiFzqAoiIiKhuYxghIiIiSTGMEBERkaQYRoiIiEhSDCNEREQkKYYRIiIikhTDCBEREUmKYYSIiIgkxTBCREREkmIYIaKH8vb2xosvvih1GURUSzGMEJmw1atXQyaTQSaTYf/+/fe9LoSAp6cnZDKZyYeF77//Hp06dYKVlRW8vLwwffp0lJSUPHS/Dz/80HAOytr+/PNPAIBer8fq1avx0ksvwdPTE7a2tmjfvj1mzpyJwsLCUscsKCjAiBEj0L59ezg6OsLOzg4+Pj747LPPUFxcfF8NiYmJePHFF+Hm5gY7Ozt07NgRn3/+OXQ6Xal+mzdvxpAhQ9CyZUvIZDI89dRTlT9hRHWIhdQFENHDWVlZYcOGDXjiiSdKtf/222+4cuUKVCqVRJVVzK5du9CrVy889dRTWLhwIY4ePYqZM2ciPT0dS5YseeC+ffr0QYsWLe5rnzRpEm7duoXAwEAAQH5+PiIjI/H4449jzJgxaNiwIeLi4jB9+nTExsbil19+gUwmA3A7jBw/fhwvvPACvL29IZfL8ddff+Gdd97BgQMHsGHDBsP7JCYmonPnzmjZsiU++OAD2NjYYNeuXRg/fjzOnTuHzz77zNB3yZIlSExMRGBgIG7cuFEVp46obhBEZLJWrVolAIg+ffoIFxcXUVxcXOr1UaNGCX9/f9GkSRPRo0ePaqvjUY/ftm1b4ePjU6r+yZMnC5lMJk6ePGn08VJTU4VMJhOjRo0ytBUVFYk///zzvr4zZswQAMTevXsfetw333xTABDXr183tI0aNUoolUpx48aNUn27du0qHBwc7qtLp9MJIYRo166d6NatmzEfi6jO4m0aIjMwcOBA3LhxA3v37jW0abVafPPNNxg0aFCZ+8ydOxedO3dG/fr1YW1tDX9/f3zzzTdl9l23bh2CgoJgY2MDZ2dndO3aFT/99NN9/fbv34+goCBYWVmhWbNm+Oqrrx5a+4kTJ3DixAmMHj0aFhZ3L8a+8cYbEEKUW9ODbNy4EUIIDB482NCmVCrRuXPn+/r27t0bAHDy5MmHHtfb2xsAkJ2dbWjTaDSwsrKCk5NTqb6NGjWCtbV1qTZPT0/I5fxrlchY/FNDZAa8vb0REhKCjRs3Gtp27dqFnJwcDBgwoMx9PvvsM/j5+eGjjz7CrFmzYGFhgX79+uHHH38s1W/GjBkYOnQoLC0t8dFHH2HGjBnw9PTEL7/8Uqrf2bNn0bdvXzz33HOYN28enJ2dMWzYMBw/fvyBtR86dAgAEBAQUKrd3d0djRs3NrxujPXr18PT0xNdu3Z9aF+1Wg0AcHFxue81rVaLzMxMXL58Gdu2bcPcuXPRpEmTUreFnnrqKWg0Grz22ms4efIkLl26hJiYGGzduhUTJ040unYiKoPUl2aIqHx3btMcPHhQLFq0SNjb24v8/HwhhBD9+vUTTz/9tBCi7Nsod/rdodVqRfv27cUzzzxjaDtz5oyQy+Wid+/ehtsLd+j1esP/N2nSRAAQv//+u6EtPT1dqFQq8X//938P/Axz5swRAERqaup9rwUGBorHH3/8gfv/27FjxwQA8f7771eof2hoqHBwcBA3b96877WNGzcKAIYtICBAHDlypFSfkpIS8eabbwpLS0tDP4VCIZYsWfLA9+VtGqKK45URIjPx6quvoqCgADt27EBubi527NhR7i0aAKVuIdy8eRM5OTl48sknkZSUZGjfvn079Ho9pk2bdt/thTuDPe9o27YtnnzyScPPDRo0QKtWrXD+/PkH1l1QUAAAZQ6ytbKyMrxeUevXrweAUrdoyjNr1iz8/PPPmD179n23WQDg6aefxt69e7FlyxaMGTMGlpaWyMvLK9VHoVCgefPmCAsLw5o1a7B582b07NkT48aNw/bt242qnYjKxqdpiMxEgwYNEBoaig0bNiA/Px86nQ59+/Ytt/+OHTswc+ZMJCcno6ioyNB+b8g4d+4c5HI52rZt+9D39/Lyuq/N2dkZN2/efOB+d0LRvTXcUVhYeN+4iwcRQmDDhg1o3749Onbs+MC+mzdvxpQpUzBixAi8/vrrZfZxdXWFq6srAKBv376YNWsWnnvuOZw5cwZubm4AgNmzZ+Ozzz7DmTNnYGdnB+B2MHz66acxduxYvPjii6XGwhCR8XhlhMiMDBo0CLt27UJMTAyef/75Mv+1DwB//PEHXnrpJVhZWeGLL77Azp07sXfvXgwaNAhCiEq9t0KhKLP9Ycdr1KgRAOD69ev3vXb9+nW4u7tXuIY///wTly5deuhVkb179yI8PBw9evRATExMhY/ft29f3Lp1C999952h7YsvvsAzzzxjCCJ3vPTSS7h27RouXrxY4eMTUdkYRojMSO/evSGXy/H3338/8BbNt99+CysrK+zZswfDhw/H888/j9DQ0Pv6NW/eHHq9HidOnKi2mn19fQEACQkJpdqvXbuGK1euGF6viPXr10Mmkz3wsx84cAC9e/dGQEAAvv76a6OuWty5ZZSTk2NoS0tLu29yMwCGydEqMnEbET0YwwiRGbGzs8OSJUvw4YcfomfPnuX2UygUkMlkpb5EL168eN8Yh169ekEul+Ojjz6CXq8v9Vplr6D8W7t27dC6dWssXbq0VD1LliyBTCYrdaspJycHp06dKhUG7iguLsaWLVvwxBNPlHnLCLj9+G6PHj3g7e2NHTt2lHsLKDMzs8zPt3z5cgCln/x57LHHsHfv3lKTmOl0Onz99dewt7dH8+bNH3IGiOhheKOTyMxEREQ8tE+PHj0wf/58dO/eHYMGDUJ6ejoWL16MFi1a4MiRI4Z+LVq0wOTJk/Hxxx/jySefRJ8+faBSqXDw4EG4u7sjOjq6SmqeM2cOXnrpJfznP//BgAEDcOzYMSxatAgjR45EmzZtDP22bduGyMhIrFq1CsOGDSt1jD179uDGjRvl3qLJzc1FWFgYbt68iffee+++R5ibN2+OkJAQALfnVYmJiUGvXr3QrFkz5ObmYs+ePdi7dy969uyJZ555xrDfhAkTMGTIEAQHB2P06NGwtrbGxo0bkZiYiJkzZ8LS0tLQ9/fff8fvv/8OAMjIyEBeXh5mzpwJAOjatWuFHkUmqpMkfZaHiB7o3kd7H6SsR3tXrFghWrZsKVQqlWjdurVYtWqVmD59uijrj/3KlSuFn5+fUKlUwtnZWXTr1q3UjKXlzcDarVu3Cj++um3bNuHr6ytUKpVo3LixmDJlitBqtWV+3lWrVt23/4ABA4SlpeV9M6HeceHChVKP6f57i4iIMPQ9ePCg6Nevn/Dy8hIqlUrY2tqKTp06ifnz5983y60QQuzevVt069ZNuLi4CKVSKTp06CBiYmLu63fn/Ja1TZ8+vULniagukglRRddiiYiIiCqBY0aIiIhIUgwjREREJCmGESIiIpIUwwgRERFJimGEiIiIJMUwQkRERJIyi0nP9Ho9rl27Bnt7+/tWEiUiIiLTJIRAbm4u3N3d71sZ/F5mEUauXbsGT09PqcsgIiKiSrh8+TIaN25c7uuVCiOLFy/GnDlzoFar4ePjg4ULFyIoKKjc/gsWLMCSJUuQmpoKFxcX9O3bF9HR0bCysqrQ+9nb2wO4/WEcHBwqUzIRERHVMI1GA09PT8P3eHmMDiObN29GVFQUYmJiEBwcjAULFiAsLAwpKSlo2LDhff03bNiACRMmYOXKlejcuTNOnz6NYcOGQSaTYf78+RV6zzu3ZhwcHBhGiIiIzMzDhlgYPYB1/vz5GDVqFCIjI9G2bVvExMTAxsYGK1euLLP/X3/9hS5dumDQoEHw9vbGf/7zHwwcOBDx8fHGvjURERHVQkaFEa1Wi8TERISGht49gFyO0NBQxMXFlblP586dkZiYaAgf58+fx86dO/HCCy+U+z5FRUXQaDSlNiIiIqqdjLpNk5mZCZ1OB1dX11Ltrq6uOHXqVJn7DBo0CJmZmXjiiScghEBJSQnGjBmDSZMmlfs+0dHRmDFjhjGlERERkZmq9qdp9u3bh1mzZuGLL75AcHAwzp49i/Hjx+Pjjz/G1KlTy9xn4sSJiIqKMvx8ZwDMg+j1emi12iqtnUyHUql84GNhRERkvowKIy4uLlAoFEhLSyvVnpaWBjc3tzL3mTp1KoYOHYqRI0cCADp06IC8vDyMHj0akydPLvMLRqVSQaVSVbgurVaLCxcuQK/XG/FpyJzI5XI0bdoUSqVS6lKIiKiKGRVGlEol/P39ERsbi169egG4fUUiNjYWb775Zpn75Ofn3xc4FAoFgNuToTwqIQSuX78OhUIBT09P/uu5Froz6d3169fh5eXFie+IiGoZo2/TREVFISIiAgEBAQgKCsKCBQuQl5eHyMhIAEB4eDg8PDwQHR0NAOjZsyfmz58PPz8/w22aqVOnomfPnoZQ8ihKSkqQn58Pd3d32NjYPPLxyDQ1aNAA165dQ0lJCSwtLaUuh4iIqpDRYaR///7IyMjAtGnToFar4evri927dxsGtaamppa6OjFlyhTIZDJMmTIFV69eRYMGDdCzZ0988sknVfIBdDodAPDyfS1359dXp9MxjBAR1TIyURX3SqqZRqOBo6MjcnJy7pv0rLCwEBcuXEDTpk0rPKMrmR/+OhMRmZ8HfX/fiwMsiIiISFIMIyZq3759kMlkyM7OlrqUB4qLine truncated
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAHHCAYAAABtF1i4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABULklEQVR4nO3deVhU1f8H8PcMMMO+CAKCILjvgGyimVaU9jP3zF3EssUljb6VZKlZimWZuaSZW7mklUu2WYZZWioK7rhvoLKICMM+w8z5/WFOTqIyCNyZ4f16nnmUM+fe+cxFmLf3nnuOTAghQERERCQRudQFEBERUd3GMEJERESSYhghIiIiSTGMEBERkaQYRoiIiEhSDCNEREQkKYYRIiIikhTDCBEREUmKYYSIiIgkxTBCRAYCAgLw1FNPSV0GEdUhDCNEJmLVqlWQyWSQyWTYvXv3Hc8LIeDn5weZTGbyYWHr1q3o0KEDbG1t4e/vj2nTpqG8vPy+202fPl1/DCp6/PXXXwb9Fy5ciFatWkGpVMLX1xdxcXEoKip6oH3qdDosXrwYwcHBsLOzg7u7Ox599FEcPnz4jnrPnTuHoUOHwtPTE3Z2dmjWrBmmTJlyR7+vv/4aHTt2hKurK9zd3dG1a1f8+OOPlTmURHWCtdQFEJEhW1tbrFu3Dg899JBB+x9//IHLly9DqVRKVFnl/Pzzz+jbty+6deuGBQsW4OjRo3jvvfeQnZ2NxYsX33Pb/v37o2nTpne0v/nmmygsLER4eLi+7Y033sAHH3yAp59+GhMnTkRqaioWLFiA48eP45dffqnSPgFg9OjRWLt2LUaOHInx48ejqKgIBw8eRHZ2tkG/Q4cOoVu3bvD19cWrr74Kd3d3pKWlIT093aDfggUL8PLLL6Nnz56YPXs2SktLsWrVKjz11FPYuHEj+vfvf89jQlQnCCIyCStXrhQARP/+/YWHh4fQaDQGz48ZM0aEhoaKRo0aiZ49e9ZYHQ+6/9atW4ugoCCD+qdMmSJkMpk4ceKE0ftLS0sTMplMjBkzRt929epVYW1tLUaMGGHQd8GCBQKA2Lp1q9H7FEKIDRs2CABi06ZN99xeq9WKtm3bisjISFFcXHzPvs2aNRPh4eFCp9Pp2/Lz84Wjo6Po3bv3Pbclqit4mYbIxAwZMgTXr1/H9u3b9W1qtRrffvsthg4dWuE2H374ITp16gR3d3fY2dkhNDQU3377bYV916xZg4iICNjb28PNzQ0PP/wwfv311zv67d69GxEREbC1tUXjxo3x5Zdf3rf21NRUpKam4vnnn4e19b8nXseOHQshxF1rupevvvoKQggMGzZM37Znzx6Ul5dj8ODBBn1vfb1+/Xqj9wkAc+fORUREBPr16wedTnfHJZ9bfv31Vxw7dgzTpk2DnZ0diouLodVqK+yrUqng6ekJmUymb3N2doajoyPs7OzuWSdRXcEwQmRiAgICEBUVha+++krf9vPPPyM/P/+OD99bPvnkE4SEhGDGjBmYNWsWrK2tMXDgwDvGJbzzzjsYMWIEbGxsMGPGDLzzzjvw8/PDjh07DPqdPXsWTz/9NB5//HF89NFHcHNzw6hRo3D8+PF71n7w4EEAQFhYmEG7j48PGjZsqH/eGGvXroWfnx8efvhhfVtZWRkA3PFhbm9vDwBITk42ep8qlQpJSUkIDw/Hm2++CRcXFzg6OqJx48b4+uuvDbb/7bffAABKpRJhYWFwcHCAvb09Bg8ejNzcXIO+3bp1w7Zt27BgwQJcvHgRJ0+exLhx45Cfn4+JEycaeTSILJTEZ2aI6B+3LtPs379fLFy4UDg5OekvAQwcOFA88sgjQoiKL6P891KBWq0Wbdu2FY8++qi+7cyZM0Iul4t+/foJrVZr0P/2SwiNGjUSAMSff/6pb8vOzhZKpVK8+uqr93wPc+bMEQBEWlraHc+Fh4eLjh073nP7/zp27JgAIF5//XWD9uTkZAFAvPvuuwbt27ZtEwCEo6Oj0ftMSUkRAIS7u7vw8vISn376qVi7dq2IiIgQMplM/Pzzz/q+vXv31vcdNmyY+Pbbb8Xbb78trK2tRadOnQyOZ1ZWlnjssccEAP3Dw8ND/P3330YdCyJLxgGsRCbomWeewaRJk/DDDz+gR48e+OGHHzB//vy79r/9DMGNGzeg1WrRpUsXg7MrW7ZsgU6nw9SpUyGXG54Uvf0SAgC0bt0aXbp00X9dv359tGjRAufPn79n3SUlJQBQ4SBbW1tbqFSqe27/X2vXrgWAOy6ndOjQAZGRkXj//ffh6+uLRx55BCdOnMBLL70EGxsbfR3G7LOwsBAAcP36dezduxeRkZEAgN69eyMwMBDvvfceevToYdA3PDwca9asAQAMGDAA9vb2iI+PR2JiIqKjowHcPFvTokULNGzYEE899RQKCgrw8ccfo3///ti1a1eFg2uJ6hpepiEyQfXr10d0dDTWrVuHTZs2QavV4umnn75r/x9++AEdO3aEra0t6tWrh/r162Px4sXIz8/X9zl37hzkcjlat25939f39/e/o83NzQ03bty453a3QtGtyyi3Ky0tNWqMhBAC69atQ9u2bdG+ffs7nt+4cSOCgoIwevRoBAYGolevXnjmmWcQEhICR0dHo/d5q7bAwEB9EAEAR0dH9OrVC0lJSfrbk2/1HTJkiME+bo3p+fvvv/VtAwcORFpaGlatWoWnn34asbGx2LlzJ9RqdYW3ARPVRQwjRCZq6NCh+Pnnn7FkyRI8+eSTcHV1rbDfrl270Lt3b9ja2uLTTz/FTz/9hO3bt2Po0KEQQlTpta2srCpsv9/+GjRoAADIyMi447mMjAz4+PhUuoa//voLly5duuMMxi2+vr7YvXs3Tp8+jT///BOXL1/GBx98gPT0dDRv3tzofd6qzcvL647nPD09odFo9ANa79bX09MTAPSh7fz589i2bRt69+5t0K9evXp46KGH7pjjhKiuYhghMlH9+vWDXC7H3r1773oXDXDzDIGtrS1++eUXjB49Gk8++aT+EsHtmjRpAp1Oh9TU1BqrOTg4GABw4MABg/arV6/i8uXL+ucrY+3atZDJZPd87wDQrFkzdOnSBd7e3khNTUVGRkaF7/9++/Tx8YG3tzeuXLlyx3NXr16Fra0tnJycAAChoaEAcEffq1evArh5ZgsAsrKyAKDCO200Gk2lJoIjqgsYRohMlKOjIxYvXozp06ejV69ed+1nZWUFmUxm8IF38eJFbNmyxaBf3759IZfLMWPGDOh0OoPnqnoG5b/atGmDli1bYunSpQb1LF68GDKZzOBSU35+Pk6ePGlwKekWjUaDb775Bg899FCFl4wqotPp8Prrr8Pe3h4vvvhilfY5aNAgpKenG9xWnZOTg++++w6PPvqofqxNnz59oFQqsXLlSoNjuWzZMgDA448/DgBo2rQp5HI5NmzYYHCML1++jF27diEkJKRS743I0nEAK5EJi4mJuW+fnj17Yu7cuejRoweGDh2K7OxsLFq0CE2bNsWRI0f0/Zo2bYopU6bg3XffRZcuXdC/f38olUrs378fPj4+SEhIqJaa58yZg969e+OJJ57A4MGDcezYMSxcuBDPPfccWrVqpe+3efNmxMbGYuXKlRg1apTBPn755Rdcv379rpdoAGDixIkoLS1FcHAwNBoN1q1bh6SkJHzxxRcVho3K7DM+Ph5ff/01BgwYgLi4OLi4uGDJkiXQaDSYNWuWvp+3tzemTJmCqVOnokePHujbty8OHz6Mzz//HEOGDNHP6lq/fn2MHj0ay5Ytw2OPPYb+/fujoKAAn376KUpKShAfH1/Zw0pk2aS8lYeI/nX7rb33UtGtvcuXLxfNmjUTSqVStGzZUqxcuVJMmzZNVPQjvmLFChESEiKUSqVwc3MTXbt2Fdu3b7/n/oUQomvXrqJr166Vei+bN28WwcHBQqlUioYNG4q33npLqNXqCt/vypUr79h+8ODBwsbGRly/fv2ur7Fy5UoRFBQkHBwchJOTk3jsscfEjh077tq/MvsUQohz586Jfv36CWdnZ2FnZyceffRRkZSUdEc/nU4nFixYIJo3by5sbGyEn59fhe9To9GIBQsWiODgYOHo6CgcHR3FI488cs9aieoamRDVdH6WiIiIqAo4ZoSIiIgkxTBCREREkmIYISIiIkkxjBAREZGkGEaIiIhIUgwjREREJCmzmPRMp9Ph6tWrcHJyumN1USIiIjJNQggUFBTAx8fnjtXCb2cWYeTq1avw8/OTugwiIiKqgvT0dDRs2PCuz5tFGLm1OFV6ejqcnZ0lroaIiIgqQ6VSwc/PT/85fjdmEUZuXZpxdnZmGCEiIjIz9xtiwQGsREREJCmGESIiIpIUwwgRERFJyizGjFSGTqeDWq2WugyqIQqF4p63hRERkfmyiDCiVqtx4cIF6HQ6qUuhGiKXyxEYGAiFQiF1KUREVM3MPowIIZCRkQErKyv4+fnxf88W6NakdxkZGfD39+fEd0REFsbsw0h5eTmKi4vh4+MDe3t7qcuhGlK/fn1cvXoV5eXlsLGxkbocIiKqRmZ/GkGr1QIAT99buFvf31vfbyIishxmH0Zu4al7y8bvLxGR5bKYMEJERETmiWHERO3cuRMymQx5eXlSl3JPiYmJaNWq1QNdPtm2bRuCg4N5NxQRUR3FMCIBmUx2z8f06dMlrS8gIADz5s2rVN/XX38db731FqysrAAABw8eREhICBwdHdGrVy/k5ubq+5aXlyM0NBRJSUkG++jRowdsbGywdu3aansPRERkPhhGJJCRkaF/zJs3D87OzgZt//vf/6q039qe9G337t04d+4cBgwYoG977rnn8OijjyIlJQX5+fmYNWuW/rmPPvoInTt3RkRExB37GjVqFObPn18rdRMR0b8Ky8rx97kcSWtgGJGAt7e3/uHi4gKZTGbQ5ujoqO+bnJyMsLAw2Nvbo1OnTjh16pT+uenLine truncated
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "# Generate and plot\n",
+    "\n",
+    "LUTs = {}\n",
+    "count = 0\n",
+    "for i in range(0, 3025, 55):\n",
+    "  if count != 0 and count < 10:\n",
+    "    count += 1\n",
+    "    continue\n",
+    "  if count > 0 and (count + 1) % 5 != 0:\n",
+    "    count +=1\n",
+    "    continue\n",
+    "\n",
+    "  # Data\n",
+    "  MACH = dataNp[:, 0][i:i + 1][0]\n",
+    "  N1 = dataNp[:, 1][i:i + 55]\n",
+    "  EPR = dataNp[:, 2][i:i+55]\n",
+    "\n",
+    "  # Interpolation\n",
+    "  interp = np.interp(n1s, N1, thrust(EPR, MACH) / max_thrust)\n",
+    "\n",
+    "  # Calculate the LUT: N1 vs. thrust using linear interpolation\n",
+    "  LUT = pd.DataFrame({\n",
+    "      'N1': n1s,\n",
+    "      'Thrust (%)': interp\n",
+    "  })\n",
+    "  LUT.loc[0] = [0,0]\n",
+    "\n",
+    "  LUTs[MACH] = LUT\n",
+    "\n",
+    "  ax = LUT.plot.line(x=\"N1\", y=\"Thrust (%)\")\n",
+    "  ax.set_title(f\"Mach {MACH}\")\n",
+    "\n",
+    "  count += 1\n",
+    "\n",
+    "plt.show()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 62,
+   "id": "c2d4f3af",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "0.000000:0.000000:0.206803:0.280661:0.354519:0.428377:0.502235:0.576093:0.649951:0.723810:0.797668,0.000000:0.000000:0.000000:0.000000:0.000000:0.000000:0.000000:0.000000:0.000000:0.000000:0.000000,20.000000:0.020334:0.018926:0.018365:0.017768:0.017149:0.016515:0.015863:0.015208:0.014534:0.013860,25.000000:0.029434:0.027422:0.026632:0.025797:0.024940:0.024059:0.023152:0.022235:0.021306:0.020375,30.000000:0.033986:0.030672:0.028919:0.026653:0.025928:0.025184:0.025802:0.021374:0.018702:0.023291,35.000000:0.070581:0.064962:0.062424:0.059648:0.056670:0.053536:0.050308:0.046924:0.043444:0.039910,40.000000:0.115706:0.107129:0.103436:0.099453:0.095210:0.090734:0.086072:0.081361:0.076653:0.071841,45.000000:0.164968:0.152931:0.147828:0.142341:0.136513:0.130440:0.124498:0.118283:0.111870:0.105388,50.000000:0.218546:0.202746:0.196107:0.188982:0.181423:0.174048:0.166338:0.158287:0.149967:0.142050,55.000000:0.276676:0.256828:0.248552:0.239686:0.230419:0.221263:0.211681:0.201686:0.191443:0.181712,60.000000:0.339274:0.315126:0.305110:0.294428:0.283152:0.272084:0.260539:0.248522:0.236179:0.224430,65.000000:0.406072:0.377362:0.365531:0.352959:0.339697:0.326347:0.312745:0.298645:0.284078:0.270104,70.000000:0.476658:0.443151:0.429436:0.414899:0.399618:0.383785:0.368032:0.351772:0.335030:0.318542,75.000000:0.550568:0.512033:0.496350:0.479778:0.462414:0.444307:0.426075:0.407561:0.388596:0.369500,80.000000:0.627228:0.583453:0.565725:0.547049:0.527537:0.507273:0.486477:0.465611:0.444349:0.422704,85.000000:0.705886:0.656705:0.636874:0.616041:0.594344:0.571897:0.548783:0.525433:0.501810:0.477902,90.000000:0.785094:0.730712:0.708789:0.685821:0.661969:0.637381:0.612132:0.586376:0.560384:0.534209,95.000000:0.863065:0.803680:0.779839:0.754920:0.729110:0.702557:0.675366:0.647478:0.619241:0.590907,100.000000:0.936071:0.872277:0.847042:0.820755:0.793551:0.765559:0.736891:0.707114:0.677066:0.646923,105.000000:0.994083:0.928996:0.903986:0.877967:0.850934:0.822172:0.792494:0.762245:0.731543:0.700495,110.000000:0.974676:0.921881:0.905331:0.888737:0.871094:0.851457:0.828909:0.802711:0.772336:0.737555,\n"
+     ]
+    }
+   ],
+   "source": [
+    "# Output for MSFS\n",
+    "\n",
+    "outputDict = {}\n",
+    "\n",
+    "for mach in LUTs:\n",
+    "  lut = LUTs[mach]\n",
+    "  _n1s = lut[\"N1\"].values\n",
+    "  scalars = lut[\"Thrust (%)\"].values\n",
+    "\n",
+    "\n",
+    "  \n",
+    "  for i in range(len(_n1s)):\n",
+    "    n1 = _n1s[i]\n",
+    "\n",
+    "    key = \"{0:6f}\".format(n1)\n",
+    "    if key not in outputDict:\n",
+    "      outputDict[key] = {}\n",
+    "\n",
+    "    _key = \"{0:6f}\".format(mach)\n",
+    "    if _key not in outputDict[key]:\n",
+    "      outputDict[key][_key] = []\n",
+    "    outputDict[key][_key] = \"{0:6f}\".format(scalars[i])\n",
+    "\n",
+    "header = True\n",
+    "output = \"0.000000\"\n",
+    "for key in outputDict:\n",
+    "  row = outputDict[key]\n",
+    "  # Machs\n",
+    "  if header:\n",
+    "    header = False\n",
+    "    for _key in row:\n",
+    "      output += \":\" + _key\n",
+    "    output += \",\"\n",
+    "  # Entries\n",
+    "  output += key\n",
+    "  for _key in row:\n",
+    "    output += \":\" + row[_key]\n",
+    "  output += \",\"\n",
+    "\n",
+    "print(output)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "5ac65794",
+   "metadata": {},
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.11.9"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}