import numpy as np
import matplotlib.pyplot as plt
from matplotlib.ticker import FormatStrFormatter
from scipy.integrate import quad

t = [0.0000,  183.0000,  366.0000,  549.0000,  732.0000,  915.0000,  
1098.0000,  1281.0000,  1464.0000,  1647.0000,  1830.0000,  2013.0000,  
2196.0000,  2379.0000,  2562.0000,  2745.0000,  2928.0000,  3111.0000,  
3294.0000,  3477.0000,  3660.0000,  3843.0000,  4026.0000,  4209.0000,  
4392.0000,  4575.0000,  4758.0000,  4941.0000,  5124.0000,  5307.0000,  
5490.0000,  5673.0000,  5856.0000,  6039.0000,  6222.0000,  6405.0000,  
6588.0000,  6771.0000,  6954.0000,  7137.0000,  7320.0000,  7503.0000,  
7686.0000,  7869.0000,  8052.0000,  8235.0000,  8418.0000,  8601.0000,  
8784.0000,  8967.0000,  9150.0000,  9333.0000,  9516.0000,  9699.0000,  
9882.0000,  10065.0000,  10248.0000,  10431.0000,  10614.0000,  10797.0000,  10980.0000, 11163.0000]

T_t_f_c_1 = [20.0000,  20.1235,  20.2347,  20.3351,  20.4258,  20.5078,  
20.5820,  20.6490,  20.7096,  20.7644,  20.8140,  20.8588,  20.8993,  
20.9359,  20.9689,  20.9989,  21.0259,  21.0503,  21.0724,  21.0924,  
21.1105,  21.1268,  21.1415,  21.1549,  21.1669,  21.1778,  21.1877,  
21.1966,  21.2046,  21.2119,  21.2184,  21.2244,  21.2297,  21.2346,  
21.2390,  21.2429,  21.2465,  21.2497,  21.2526,  21.2553,  21.2576,  
21.2598,  21.2617,  21.2635,  21.2650,  21.2665,  21.2678,  21.2689,  
21.2700,  21.2709,  21.2717,  21.2725,  21.2732,  21.2738,  21.2744,  
21.2749,  21.2753,  21.2757,  21.2761,  21.2764,  21.2767,  21.2769]

T_t_f_c_2 = [20.0000,  21.2355,  22.3467,  23.3506,  24.2579,  25.0783,  
 25.8199,  26.4904,  27.0965,  27.6445,  28.1399,  28.5877,  28.9926, 
 29.3586,  29.6895,  29.9886,  30.2590,  30.5034,  30.7244,  30.9241,  
 31.1047,  31.2679,  31.4154,  31.5488,  31.6693,  31.7782,  31.8767,  
 31.9657,  32.0461,  32.1188,  32.1844,  32.2438,  32.2974,  32.3459,  
 32.3897,  32.4292,  32.4649,  32.4972,  32.5263,  32.5526,  32.5764,  
 32.5979,  32.6172,  32.6347,  32.6505,  32.6647,  32.6775,  32.6891,  
 32.6995,  32.7089,  32.7174,  32.7250,  32.7319,  32.7381,  32.7436,  
 32.7486,  32.7531,  32.7571,  32.7607,  32.7639,  32.7668,  32.7694]


T_t_f_c_3 = [20.0000,  44.7095,  66.9347,  87.0111,  105.1590,  
121.5654,  136.3978,  149.8071,  161.9297,  172.8892,  182.7970,  
191.7540,  199.8514,  207.1716,  213.7892,  219.7715,  225.1795,  
230.0682,  234.4875,  238.4823,  242.0934,  245.3576,  248.3081,  
250.9751,  253.3856,  255.5644,  257.5336,  259.3134,  260.9218,  
262.3754,  263.6890,  264.8759,  265.9485,  266.9175,  267.7931,  
268.5840,  269.2985,  269.9438,  270.5267,  271.0530,  271.5282,  
271.9573,  272.3446,  272.6941,  273.0095,  273.2941,  273.5507,  
273.7822,  273.9908,  274.1788,  274.3481,  274.5006,  274.6379,  
274.7614,  274.8724,  274.9722,  275.0618,  275.1422,  275.2143,  275.2788,  275.3366,  275.3882]



chart = plt.figure()    
#plt.plot(t, T_t_f_c_1,'r', linewidth=1)
#plt.plot(t, T_t_f_c_2,'g', linewidth=1)
plt.plot(t, T_t_f_c_3,'r', linewidth=1)
plt.xlabel(r'Tempo $(s)$')
plt.ylabel(r'Temperatura $(^oC)$')
plt.grid(True)

chart.savefig("/Users/Libotte/Desktop/chart_citt.pdf", bbox_inches='tight')

plt.show()
54ab168e7cd3c52438cca72decb086171c449fdf
<mask token>

<mask token>
plt.plot(t, T_t_f_c_3, 'r', linewidth=1)
plt.xlabel('Tempo $(s)$')
plt.ylabel('Temperatura $(^oC)$')
plt.grid(True)
chart.savefig('/Users/Libotte/Desktop/chart_citt.pdf', bbox_inches='tight')
plt.show()

<mask token>
t = [0.0, 183.0, 366.0, 549.0, 732.0, 915.0, 1098.0, 1281.0, 1464.0, 1647.0,
    1830.0, 2013.0, 2196.0, 2379.0, 2562.0, 2745.0, 2928.0, 3111.0, 3294.0,
    3477.0, 3660.0, 3843.0, 4026.0, 4209.0, 4392.0, 4575.0, 4758.0, 4941.0,
    5124.0, 5307.0, 5490.0, 5673.0, 5856.0, 6039.0, 6222.0, 6405.0, 6588.0,
    6771.0, 6954.0, 7137.0, 7320.0, 7503.0, 7686.0, 7869.0, 8052.0, 8235.0,
    8418.0, 8601.0, 8784.0, 8967.0, 9150.0, 9333.0, 9516.0, 9699.0, 9882.0,
    10065.0, 10248.0, 10431.0, 10614.0, 10797.0, 10980.0, 11163.0]
T_t_f_c_1 = [20.0, 20.1235, 20.2347, 20.3351, 20.4258, 20.5078, 20.582, 
    20.649, 20.7096, 20.7644, 20.814, 20.8588, 20.8993, 20.9359, 20.9689, 
    20.9989, 21.0259, 21.0503, 21.0724, 21.0924, 21.1105, 21.1268, 21.1415,
    21.1549, 21.1669, 21.1778, 21.1877, 21.1966, 21.2046, 21.2119, 21.2184,
    21.2244, 21.2297, 21.2346, 21.239, 21.2429, 21.2465, 21.2497, 21.2526, 
    21.2553, 21.2576, 21.2598, 21.2617, 21.2635, 21.265, 21.2665, 21.2678, 
    21.2689, 21.27, 21.2709, 21.2717, 21.2725, 21.2732, 21.2738, 21.2744, 
    21.2749, 21.2753, 21.2757, 21.2761, 21.2764, 21.2767, 21.2769]
T_t_f_c_2 = [20.0, 21.2355, 22.3467, 23.3506, 24.2579, 25.0783, 25.8199, 
    26.4904, 27.0965, 27.6445, 28.1399, 28.5877, 28.9926, 29.3586, 29.6895,
    29.9886, 30.259, 30.5034, 30.7244, 30.9241, 31.1047, 31.2679, 31.4154, 
    31.5488, 31.6693, 31.7782, 31.8767, 31.9657, 32.0461, 32.1188, 32.1844,
    32.2438, 32.2974, 32.3459, 32.3897, 32.4292, 32.4649, 32.4972, 32.5263,
    32.5526, 32.5764, 32.5979, 32.6172, 32.6347, 32.6505, 32.6647, 32.6775,
    32.6891, 32.6995, 32.7089, 32.7174, 32.725, 32.7319, 32.7381, 32.7436, 
    32.7486, 32.7531, 32.7571, 32.7607, 32.7639, 32.7668, 32.7694]
T_t_f_c_3 = [20.0, 44.7095, 66.9347, 87.0111, 105.159, 121.5654, 136.3978, 
    149.8071, 161.9297, 172.8892, 182.797, 191.754, 199.8514, 207.1716, 
    213.7892, 219.7715, 225.1795, 230.0682, 234.4875, 238.4823, 242.0934, 
    245.3576, 248.3081, 250.9751, 253.3856, 255.5644, 257.5336, 259.3134, 
    260.9218, 262.3754, 263.689, 264.8759, 265.9485, 266.9175, 267.7931, 
    268.584, 269.2985, 269.9438, 270.5267, 271.053, 271.5282, 271.9573, 
    272.3446, 272.6941, 273.0095, 273.2941, 273.5507, 273.7822, 273.9908, 
    274.1788, 274.3481, 274.5006, 274.6379, 274.7614, 274.8724, 274.9722, 
    275.0618, 275.1422, 275.2143, 275.2788, 275.3366, 275.3882]
chart = plt.figure()
plt.plot(t, T_t_f_c_3, 'r', linewidth=1)
plt.xlabel('Tempo $(s)$')
plt.ylabel('Temperatura $(^oC)$')
plt.grid(True)
chart.savefig('/Users/Libotte/Desktop/chart_citt.pdf', bbox_inches='tight')
plt.show()

import numpy as np
import matplotlib.pyplot as plt
from matplotlib.ticker import FormatStrFormatter
from scipy.integrate import quad
t = [0.0, 183.0, 366.0, 549.0, 732.0, 915.0, 1098.0, 1281.0, 1464.0, 1647.0,
    1830.0, 2013.0, 2196.0, 2379.0, 2562.0, 2745.0, 2928.0, 3111.0, 3294.0,
    3477.0, 3660.0, 3843.0, 4026.0, 4209.0, 4392.0, 4575.0, 4758.0, 4941.0,
    5124.0, 5307.0, 5490.0, 5673.0, 5856.0, 6039.0, 6222.0, 6405.0, 6588.0,
    6771.0, 6954.0, 7137.0, 7320.0, 7503.0, 7686.0, 7869.0, 8052.0, 8235.0,
    8418.0, 8601.0, 8784.0, 8967.0, 9150.0, 9333.0, 9516.0, 9699.0, 9882.0,
    10065.0, 10248.0, 10431.0, 10614.0, 10797.0, 10980.0, 11163.0]
T_t_f_c_1 = [20.0, 20.1235, 20.2347, 20.3351, 20.4258, 20.5078, 20.582, 
    20.649, 20.7096, 20.7644, 20.814, 20.8588, 20.8993, 20.9359, 20.9689, 
    20.9989, 21.0259, 21.0503, 21.0724, 21.0924, 21.1105, 21.1268, 21.1415,
    21.1549, 21.1669, 21.1778, 21.1877, 21.1966, 21.2046, 21.2119, 21.2184,
    21.2244, 21.2297, 21.2346, 21.239, 21.2429, 21.2465, 21.2497, 21.2526, 
    21.2553, 21.2576, 21.2598, 21.2617, 21.2635, 21.265, 21.2665, 21.2678, 
    21.2689, 21.27, 21.2709, 21.2717, 21.2725, 21.2732, 21.2738, 21.2744, 
    21.2749, 21.2753, 21.2757, 21.2761, 21.2764, 21.2767, 21.2769]
T_t_f_c_2 = [20.0, 21.2355, 22.3467, 23.3506, 24.2579, 25.0783, 25.8199, 
    26.4904, 27.0965, 27.6445, 28.1399, 28.5877, 28.9926, 29.3586, 29.6895,
    29.9886, 30.259, 30.5034, 30.7244, 30.9241, 31.1047, 31.2679, 31.4154, 
    31.5488, 31.6693, 31.7782, 31.8767, 31.9657, 32.0461, 32.1188, 32.1844,
    32.2438, 32.2974, 32.3459, 32.3897, 32.4292, 32.4649, 32.4972, 32.5263,
    32.5526, 32.5764, 32.5979, 32.6172, 32.6347, 32.6505, 32.6647, 32.6775,
    32.6891, 32.6995, 32.7089, 32.7174, 32.725, 32.7319, 32.7381, 32.7436, 
    32.7486, 32.7531, 32.7571, 32.7607, 32.7639, 32.7668, 32.7694]
T_t_f_c_3 = [20.0, 44.7095, 66.9347, 87.0111, 105.159, 121.5654, 136.3978, 
    149.8071, 161.9297, 172.8892, 182.797, 191.754, 199.8514, 207.1716, 
    213.7892, 219.7715, 225.1795, 230.0682, 234.4875, 238.4823, 242.0934, 
    245.3576, 248.3081, 250.9751, 253.3856, 255.5644, 257.5336, 259.3134, 
    260.9218, 262.3754, 263.689, 264.8759, 265.9485, 266.9175, 267.7931, 
    268.584, 269.2985, 269.9438, 270.5267, 271.053, 271.5282, 271.9573, 
    272.3446, 272.6941, 273.0095, 273.2941, 273.5507, 273.7822, 273.9908, 
    274.1788, 274.3481, 274.5006, 274.6379, 274.7614, 274.8724, 274.9722, 
    275.0618, 275.1422, 275.2143, 275.2788, 275.3366, 275.3882]
chart = plt.figure()
plt.plot(t, T_t_f_c_3, 'r', linewidth=1)
plt.xlabel('Tempo $(s)$')
plt.ylabel('Temperatura $(^oC)$')
plt.grid(True)
chart.savefig('/Users/Libotte/Desktop/chart_citt.pdf', bbox_inches='tight')
plt.show()

import numpy as np
import matplotlib.pyplot as plt
from matplotlib.ticker import FormatStrFormatter
from scipy.integrate import quad

t = [0.0000,  183.0000,  366.0000,  549.0000,  732.0000,  915.0000,  
1098.0000,  1281.0000,  1464.0000,  1647.0000,  1830.0000,  2013.0000,  
2196.0000,  2379.0000,  2562.0000,  2745.0000,  2928.0000,  3111.0000,  
3294.0000,  3477.0000,  3660.0000,  3843.0000,  4026.0000,  4209.0000,  
4392.0000,  4575.0000,  4758.0000,  4941.0000,  5124.0000,  5307.0000,  
5490.0000,  5673.0000,  5856.0000,  6039.0000,  6222.0000,  6405.0000,  
6588.0000,  6771.0000,  6954.0000,  7137.0000,  7320.0000,  7503.0000,  
7686.0000,  7869.0000,  8052.0000,  8235.0000,  8418.0000,  8601.0000,  
8784.0000,  8967.0000,  9150.0000,  9333.0000,  9516.0000,  9699.0000,  
9882.0000,  10065.0000,  10248.0000,  10431.0000,  10614.0000,  10797.0000,  10980.0000, 11163.0000]

T_t_f_c_1 = [20.0000,  20.1235,  20.2347,  20.3351,  20.4258,  20.5078,  
20.5820,  20.6490,  20.7096,  20.7644,  20.8140,  20.8588,  20.8993,  
20.9359,  20.9689,  20.9989,  21.0259,  21.0503,  21.0724,  21.0924,  
21.1105,  21.1268,  21.1415,  21.1549,  21.1669,  21.1778,  21.1877,  
21.1966,  21.2046,  21.2119,  21.2184,  21.2244,  21.2297,  21.2346,  
21.2390,  21.2429,  21.2465,  21.2497,  21.2526,  21.2553,  21.2576,  
21.2598,  21.2617,  21.2635,  21.2650,  21.2665,  21.2678,  21.2689,  
21.2700,  21.2709,  21.2717,  21.2725,  21.2732,  21.2738,  21.2744,  
21.2749,  21.2753,  21.2757,  21.2761,  21.2764,  21.2767,  21.2769]

T_t_f_c_2 = [20.0000,  21.2355,  22.3467,  23.3506,  24.2579,  25.0783,  
 25.8199,  26.4904,  27.0965,  27.6445,  28.1399,  28.5877,  28.9926, 
 29.3586,  29.6895,  29.9886,  30.2590,  30.5034,  30.7244,  30.9241,  
 31.1047,  31.2679,  31.4154,  31.5488,  31.6693,  31.7782,  31.8767,  
 31.9657,  32.0461,  32.1188,  32.1844,  32.2438,  32.2974,  32.3459,  
 32.3897,  32.4292,  32.4649,  32.4972,  32.5263,  32.5526,  32.5764,  
 32.5979,  32.6172,  32.6347,  32.6505,  32.6647,  32.6775,  32.6891,  
 32.6995,  32.7089,  32.7174,  32.7250,  32.7319,  32.7381,  32.7436,  
 32.7486,  32.7531,  32.7571,  32.7607,  32.7639,  32.7668,  32.7694]


T_t_f_c_3 = [20.0000,  44.7095,  66.9347,  87.0111,  105.1590,  
121.5654,  136.3978,  149.8071,  161.9297,  172.8892,  182.7970,  
191.7540,  199.8514,  207.1716,  213.7892,  219.7715,  225.1795,  
230.0682,  234.4875,  238.4823,  242.0934,  245.3576,  248.3081,  
250.9751,  253.3856,  255.5644,  257.5336,  259.3134,  260.9218,  
262.3754,  263.6890,  264.8759,  265.9485,  266.9175,  267.7931,  
268.5840,  269.2985,  269.9438,  270.5267,  271.0530,  271.5282,  
271.9573,  272.3446,  272.6941,  273.0095,  273.2941,  273.5507,  
273.7822,  273.9908,  274.1788,  274.3481,  274.5006,  274.6379,  
274.7614,  274.8724,  274.9722,  275.0618,  275.1422,  275.2143,  275.2788,  275.3366,  275.3882]



chart = plt.figure()    
#plt.plot(t, T_t_f_c_1,'r', linewidth=1)
#plt.plot(t, T_t_f_c_2,'g', linewidth=1)
plt.plot(t, T_t_f_c_3,'r', linewidth=1)
plt.xlabel(r'Tempo $(s)$')
plt.ylabel(r'Temperatura $(^oC)$')
plt.grid(True)

chart.savefig("/Users/Libotte/Desktop/chart_citt.pdf", bbox_inches='tight')

plt.show()