You cannot select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
This repo is archived. You can view files and clone it, but cannot push or open issues/pull-requests.

419 lines
24 KiB
TeX

This file contains ambiguous Unicode characters!

This file contains ambiguous Unicode characters that may be confused with others in your current locale. If your use case is intentional and legitimate, you can safely ignore this warning. Use the Escape button to highlight these characters.

\documentclass[a4paper,12pt]{article}
\usepackage{cmap}
\usepackage[utf8]{inputenc}
\usepackage[warn]{mathtext}
\usepackage{epsf,amsmath,amsfonts,amssymb,amsbsy}
\usepackage[mathscr]{eucal}
\usepackage[english, russian]{babel}
\author{Гришаев Григорий С01-119}
\title{Отчёт о выполнении лабораторной работы 2.1.1}
\usepackage[left=2cm,right=2cm,top=2cm,bottom=2cm]{geometry}
\usepackage{graphicx}
\usepackage{indentfirst}
\graphicspath{{images/}}
\DeclareGraphicsExtensions{.pdf,.png,.jpg}
\usepackage{pgfplots}
\begin{document}
\maketitle
\begin{center}
{\Large Измерение удельной теплоёмкости воздуха при постоянном давленини}
\end{center}
\paragraph*{Цель работы:} 1) измерение повышения температуры воздуха в результате подвода тепла при стационарном течении через стеклянную трубу; 2) вычисление по результатам измерений теплоёмкости воздуха при постоянном давлении.
\paragraph*{В работе используются:} теплоизолированная трубка; электронагреватель; источник питания постоянного тока; амперметр; вольтметр; термопара; компрессор; газовый счётчик;
секундомер.
\section{Теоретические сведения}
Определение теплоёмкости обычно производится в калориметрах. При этом регистрируется количество тепла {Q}, полученное телом, и изменение температуры этого тела $ {\Delta T} $.
Теплоёмкость тела определяется как их отношение: $$ C = \frac{Q}{\Delta T} \; (1).$$
Необходимо, чтобы количество тепла, затрачиваемое на нагревание исследуемого тела, существенно превосходило тепло, расходуемое на нагревание самого калориметра, а также на потери тепла из установки.
Для увеличения количества нагреваемого газа при неизменных размерах установки в нашей работе исследуемый газ (воздух) продувается через калориметр, внутри которого установлен нагреватель. При этом
измеряются мощность нагревателя, масса воздуха, протекающего в единицу
времени (расход), и приращение его температуры.
\subsection{Экспериментальная установка:}
Схема установки изображена на рис.1. Воздух, нагнетаемый компрессором, прокачивается через калориметр. Калориметр представляет собой стеклянную цилиндрическую трубку с двойными стенками, запаянными с торцов.
\begin{figure}[h]
\center{\includegraphics{1}}
\end{figure}
Нагреватель в виде намотанной на пенопласт нихромовой проволоки расположен внутри калориметра непосредственно в воздушном потоке. Нагрев проволоки производится от регулируемого источника постоянного тока (ИП).
Напряжение $U$ на нагревателе и ток $I$ через него регистрируются цифровыми мультиметрами. Таким образом, мощность нагрева равна
$$N= UI \; (3).$$
Для измерения разности температур $\Delta T$ служит медно-константановая
термопара. Один спай термопары расположен в струе воздуха, входящего в
калориметр, и находится при комнатной температуре, а второй — в струе выходящего нагретого воздуха. Константановая проволока термопары расположена внутри калориметра, а медные проводники подключены к цифровому вольтметру. Возникающая в термопаре ЭДС $\varepsilon$ пропорциональна разности температур $\Delta T$ спаев: $$\varepsilon =\beta \Delta T \; (4),$$ где $\beta = 42.3 \frac{мкВ}{^\circ C}$ — чувствительность медно-константановой термопары в рабочем диапазоне температур (2030 $^\circ C$ ). ЭДС регистрируется с помощью микровольтметра.
Вычислим работу, совершаемую при протекании газа через калориметр. Внешняя работа по перемещению моля газа в направлении течения в начале трубки равна $ A_{1} = P_{1}V_{1} $, а в конце трубки давление препятствует движению и внешняя работа над газом отрицательна: $ A_{2}=-P_{2}V_{2}$. Полная работа над газом равна $ A_{1}+A_{2}=P_{1}V_{1}-P_{2}V_{2}$, а работа самого газа равна этой же величине, но с обратным знаком: $$ A =P_{2}V_{2}-P_{1}V_{1}$$
Внутренняя энергия газа изменяется на величину $ \Delta U = U_{2} - U_{1} $.
Из первого начала термодинамики: $$ Q = U_{2}-U_{1}+P_{2}V_{2}-P_{1}V_{1}=H_{2}-H_{1}$$,
где $H=U+PV$ - энтальпия.
Для идеального газа $ H = C_{p}T$, поэтому $$ Q = C_{p}(T_{2}-T_{1})$$.
Следовательно, в данном эксперименте измеряется теплоемкость при постоянном давлении.
Расчет удельной теплоемкости воздуха: $$ c_{p} = \frac{Q}{m \Delta T}= \frac{IU-N}{m \Delta T}$$,
% Объём воздуха, прошедшего через калориметр, измеряется газовым счётчиком ГС. Для регулировки расхода служит кран К. Время $\Delta t$ прохождения некоторого объема $\Delta V$ воздуха измеряется секундомером. Объёмный расход равен $\frac{\Delta V}{\Delta t} $, массовый расход может быть найден как $$q = \rho_{0} \frac{\Delta V}{\Delta t} \; (5),$$ где $rho_{0}$ — плотность воздуха при комнатной температуре, которая в свою очередь может быть получена из уравнения Менделеева–Клапейрона: $\rho_{0}= \frac{\mu P_{0} }{R T_{0}},$ где $P_{0}$ — атмосферное давление, $T_{0}$ — комнатная температура (в Кельвинах), $\mu = 29,0 {г/моль}$ — средняя молярная масса (сухого) воздуха.
% Учитывая особенности устройства калориметра, следует ожидать, что мощность нагревателя расходуется не только на нагрев массы прокачиваемого воздуха, но и частично теряется за счет нагрева внутренних стенок термостата и рассеяния тепла через торцы термостата. Можно предположить, что при небольшом нагреве ($\Delta T \ll T_{0}$) мощность потерь тепла $N_{пот}$ прямо пропорциональна разности температур:$$ N_{пот} = \alpha \Delta T \; (6),$$ где $\alpha$ — некоторая константа. При этом условии основное соотношение (2) принимает вид $$N = (c_{P}q +\alpha)\Delta T \;(7)$$
% Следовательно, при фиксированном расходе воздуха ($q = const$ ) подводимая мощность и разность температур связаны прямой пропорциональностью($\Delta T(N)$ — линейная функция).
\subsection{Ход работы}
% \begin{enumerate}
% \item Подготовим к работе газовый счетчик: проверим, что он заполнен водой, установим счетчик по уровню.
% \item Охладим калориметр до комнатной температуры.
% \item Включим вольтметр, предназначенный для измерения ЭДС термопары.
% \item Запишем показания компнатной температуры и давления. $$T_{0} = 297.05 \; ^\circ C, P_{0} = 99325 \pm 13 \; {Па} $$
% \item С помощью газового счетчика и секундомера измерим максимальный расход воздуха $\frac{\Delta V}{\Delta T}$ (в л/с). Измерения представлены в таблице 1. По найденным значениям определим среднее значение расхода и массовый расход воздуха $q_{max}$ [г/с].
%
% $$q = \rho_0 \frac{\Delta V}{\Delta t} = \frac{\mu P_0}{RT_0} \frac{\Delta V}{\Delta t}.$$
%
% Относительная погрешность косвенных измерений может быть найдена по формуле $$\frac{\sigma_{q_{max}кос}}{q_{max}} = \sqrt{(\frac{\sigma_{T_0}}{T_{0}})^2+(\frac{\sigma_{P_0}}{P_{0}})^2+ (\frac{\sigma_t}{t})^2}$$
%
%
% \begin{table}
% \begin{center}
% \begin{tabular}{|c|c|c|c|c|}
% \hline
% $\Delta V, л$ & $\Delta t, c$ & $\frac{\Delta V}{\Delta t},\frac{л}{с}$ & $q_{max},\frac{г}{c} $ & $\sigma_{q_{max}кос}, \; \frac{г}{с} \cdot 10^{-3}$ \\
% \hline
% x & 0 & 0 & 0 & 0 \\
% \hline
% x & 0 & 0 & 0 & 0\\
% \hline
% x & 0 & 0 & 0 & 0\\
% \hline
% x & 0 & 0 & 0 &0 \\
% \hline
%
% \end{tabular}
% \end{center}
% \caption{Измерение расхода воздуха}
% \end{table}
% $$\overline{\frac{\Delta V}{\Delta t}} = 0.1951 \; \frac{л}{с} ,\; \overline{q_{max}} = 0.2276 \; \frac{г}{с} $$
%
% Случайная погрешность массового расхода может быть найдена по формуле: $$\sigma_{q_{max}сл} = \sqrt{\frac{\sum_{i=1}^{7} (q_{max,i}-\overline{q_{max}})^2}{6}} = 0.0003 \; \frac{г}{c}.$$
%
% Косвенная погрещность для среднего значения: $q_{max}$ $$\sigma_{\overline{q_{max}}кос} = \sqrt{\frac{\sum_{i=1}^{7} (\sigma_{q_{max}кос})^2}{7^2}} = 0.0003 \; \frac{г}{c}.$$
%
% Суммарная погрешность: $$\sigma_{\overline{q_{max}}} = \sqrt{(\sigma_{q_{max}сл})+(\sigma_{\overline{q_{max}}кос})^2} = 0.0004 \; \frac{г}{c}.$$
%
% Окончательное значение: $$q_{max} = 0.2276 \pm 0.0004 \; \frac{г}{с} $$
%
% \item Оценим величину тока нагревателя $I_{0}$, требуемого для нагрева воздуха на $\delta T = 1 {К}$.
%
% Определим теоретическое значение удельной теплоемкости воздуха при постоянном давлении $C_{теорp} \; \frac{Дж}{г\cdot K}$, считая воздух смесью двухатомных идеальных газов: $Cp = 3.5R\mu \approx 1 \; \frac{Дж}{г\cdot K}.$
%
% Оценим минимальную мощность $N_0$, необходимую для нагрева газа при максимальном расходе. $N_{0} = c_{p}q\Delta T \approx 0.227 {Вт}.$
%
% Учитывая, что сопротивление проволоки нагревателя составляет приблизительно $R_{н} \approx 35 {Ом}$ и в процессе опыта практически не меняется, искомое значение тока $I_{0} = q N_{0} R_{н} \approx 0.08 \; {А}.$
%
% \item Проведем измерение зависимости разности температур от мощности нагрева $\Delta T(N)$ при максимальном расходе воздуха $q_0 = q_{max}.$
% \begin{table}
%
% \begin{center}
% \begin{tabular}{|c|c|c|c|c|c|c|c|c|}
% \hline
% $I, мA$ & $U, B$ & $N, Вт$ & $R_н, Ом$ & $\varepsilon, \muВ$ & $ \Delta T, K$
% \\
% \hline
% 108.18 & 3.19 & 0.345 & 29.49 & 52 & 1.28
% \\
% \hline
% 141.62 & 4.18 & 0.592 & 29.52 & 88 & 2.16
% \\
% \hline
% 165.70 & 4.89 & 0.810 & 29.51 & 118 & 2.90
% \\
% \hline
% 181.06 & 5.34 & 0.967 & 29.49 & 141 & 3.46
% \\
% \hline
% 210.7 & 6.21 & 1.308 & 29.48 & 188 & 4.61
% \\
% \hline
% \end{tabular}
% \end{center}
% \caption{Измерение $\Delta T (N) \; {при} \; q_{max}$}
% \end{table}
% Следует отметить, что погрешность измерения тока: $\sigma_{I} = 0.01 \; мA$, а напряжения: $\sigma_{U}= 0.01 \; В$, $\sigma_{\varepsilon}= 1 \; \muВ$
%
% Завершив первую серию измерении, охладим калориметр до комнатнои температуры.
% Для этого отключим источник питания нагревателя, откроем кран К и продуем калориметр при максимальном расходе воздуха до тех пор, пока показания ЭДС не достигнут нуля.
%
% Данные представлены в таблице 2.
%
% \begin{table}
% \begin{center}
% \begin{tabular}{lr}
% \begin{tabular}{|c|c|c|c|c|}
% \hline
% $\Delta V, л$ & $\Delta t, c$ & $\frac{\Delta V}{\Delta t},\frac{л}{с}$ & $q_1, \frac{г}{c}$ & $Z, \; \frac{г}{с} \cdot 10^{-3}$ \\
% \hline
% 5 & 46.02 & 0.1086 & 0.1266 & 0.550 \\
% \hline
% 10 & 91.4 & 0.1094 & 0.1275 & 0.279 \\
% \hline
% 15 & 134.78 & 0.1113 & 0.1297& 0.192 \\
% \hline
% 20 & 178.96 & 0.1117 & 0.1302 & 0.146\\
% \hline
% 25 & 221.14 & 0.1130 & 0.1317 & 0.119 \\
% \hline
% 30 & 266.36 & 0.1126 & 0.1312 & 0.099\\
% \hline
% 35 & 311.83 & 0.1122 & 0.1308 & 0.084\\
% \hline
% \end{tabular}
%
% \begin{tabular}{|c|c|c|c|c|}
% \hline
% $\Delta V, л$ & $\Delta t, c$ & $\frac{\Delta V}{\Delta t},\frac{л}{с}$ & $q_2, \frac{г}{c}$ & $Z, \; \frac{г}{с} \cdot 10^{-3}$ \\
% \hline
% 1 & 13.40 & 0.0747 & 0.0870 & 1.30\\
% \hline
% 2 & 27.56 & 0.0726 & 0.0846 & 0.61\\
% \hline
% 3 & 41.31 & 0.0726 & 0.0846 & 0.41\\
% \hline
% 4 & 54.95 & 0.0733 & 0.0855 & 0.31 \\
% \hline
% 5 & 69.05 & 0.0724 & 0.0844 & 0.24\\
% \hline
% 6 & 82.24 & 0.0729 & 0.0850 & 0.21\\
% \hline
% 7 & 96.18 & 0.0728 & 0.0848 & 0.18 \\
% \hline
%\end{tabular}
%\end{tabular}
%\end{center}
%\caption{Измерения других расходов; $q_1 = 0.123 \pm 0.002 \; \frac{г}{с}; \; q_2 = 0.085 \pm 0.001 \; \frac{г}{с} $ }
%\end{table}
% Проведем аналогичные измерения для других значений расхода воздуха. Новая температура $T_{0} = 297.4 ^\circ C$.
% Данные представлены в таблице 3 и 4. ($Z \equiv \sigma_{q_{max}кос}$). Погрешности рассчитаны аналогично $q_{max}.$
%
%\begin{table}
% \begin{center}
% \begin{tabular}{lr}
% \begin{tabular}{|c|c|c|c|c|c|}
% \hline
% $I, A$ & $U, B$ & $N, Вт$ & $R_н, Ом$ & $\varepsilon, \muВ$ & $ \Delta T, K$\\
% \hline
% 91.34 & 2.69 & 0.246 & 29.45 & 50 & 1.23 \\
% \hline
% 120.30 & 3.55 & 0.427 & 29.51 & 89 & 2.19 \\
% \hline
% 145.11 & 4.28 & 0.621 & 29.49 & 128 & 3.14 \\
% \hline
% 176.45 & 5.20 & 0.918 & 29.47 & 184 & 4.52 \\
% \hline
% 212.20 & 6.23 & 1.322 & 29.36 & 261 & 6.41 \\
% \hline
%
% \end{tabular}
% \begin{tabular}{|c|c|c|c|c|c|}
% \hline
% $I, A$ & $U, B$ & $N, Вт$ & $R_н, Ом$ & $\varepsilon, \muВ$ & $ \Delta T, K$\\
% \hline
% 105.95 & 3.13 & 0.332 & 29.54 & 81 & 1.99 \\
% \hline
% 125.74 & 3.71 & 0.466 & 29.51 & 128 & 3.14 \\
% \hline
% 149.46 & 4.41 & 0.659 & 29.50 & 175 & 4.30 \\
% \hline
% 177.33 & 5.23 & 0.927 & 29.49 & 236 & 5.80 \\
% \hline
% 211.5 & 6.25 & 1.322 & 29.55 & 310 & 7.62 \\
% \hline
% \end{tabular}
% \end{tabular}
% \end{center}
% \caption{Измерение $\Delta T (N) \; {при} \; q_{1} \; и \; q_{2}$}
%\end{table}
%
% Погрешности будем считать по слудующим формулам: $$ \sigma_{\Delta T} = \Delta T \frac{\sigma_{\varepsilon}}{\varepsilon}, \sigma_{N}= N\sqrt{( \frac{\sigma_{I}}{I})^2 + (\frac{\sigma_{U}}{U})^2} \approx N \frac{\sigma_{U}}{U}$$
%
% \begin{table}
% \begin{tabular}{lr}
% \begin{tabular}{|l|l|l|l|l|l|l|l|l|l|l|}
% \hline
% q, $\frac{г}{с}$ & \multicolumn{4}{|c|}{0.2276} & \multicolumn{4}{|c|}{0.123}
% \\
% \hline
% k & $\Delta T, \; ^\circ C$& $\sigma_{\Delta T}\; ^\circ C$ & $N$, мВт & $\sigma_{N}$, мВт & $\Delta T, \; ^\circ C$& $\sigma_{\Delta T}\; ^\circ C$ & $N$, мВт & $\sigma_{N}$, мВт
% \\
%
% \hline
% 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0
% \\
% \hline
% 1 & 1.28 & 0.02 & 345 & 1.1 & 1.23 & 0.02 & 246 & 0.9
% \\
% \hline
% 2 & 2.16 & 0.02 & 592 & 1.4 & 2.19 & 0.02 & 427 & 1.2
% \\
% \hline
% 3 & 2.90 & 0.02 & 810 & 1.7 & 3.14 & 0.02 & 621 & 1.5
% \\
% \hline
% 4 & 3.46 & 0.02 & 967 & 1.8 & 4.52 & 0.02 & 918 & 1.8
% \\
% \hline
% 5 & 4.61 & 0.02 & 1308 & 2.1 & 6.41 & 0.02 & 1322 & 2.1
% \\
% \hline
% \end{tabular}
%
% \\
%
% \begin{tabular} {|l|l|l|l|l|l|l|l|l|l|l|}
% \hline
% q, $\frac{г}{с}$ & \multicolumn{4}{|c|}{0.085}
% \\
% \hline
% k & $\Delta T, \; ^\circ C$& $\sigma_{\Delta T}\; ^\circ C$ & $N$, мВт & $\sigma_{N}$, мВт
% \\
% \hline
% 0 & 0 & 0 & 0 & 0
% \\
% \hline
% 1 & 1.99 & 0.02 & 332 & 1.1
% \\
% \hline
% 2 & 3.14 & 0.02 & 466 & 1.3
% \\
% \hline
% 3 & 4.30 & 0.02 & 659 & 1.5
% \\
% \hline
% 4 & 5.80 & 0.02 & 927 & 1.8
% \\
% \hline
% 5 & 7.62 & 0.02 & 1322 & 2.1
% \\
% \hline
% \end{tabular}
%
%
%
% \end{tabular}
% \caption{ погрешности, $\Delta T \; и \; N$}
% \end{table}
%
%
% После завершения опытов выключим источник питания нагревателя и мультиметры. Кран К откроем для максимального продува воздуха через калориметр.
%
%\item Построим на одном графике зависимости $\Delta T (N)$ при разных значениях $q$.
% \begin{center}
% \begin{tikzpicture}[scale = 2]
% \begin{axis}[
% axis lines = left,
% legend style={at={(1,0.3)}},
% xlabel = {$N$, {мВт}},
% ylabel = {$\Delta T, {^\circ C} $
% },
% xmin=0, xmax=1350,
% ymin=0, ymax=8,
% ymajorgrids = true,
% xmajorgrids = true,
% minor tick num = 4
% ]
% \addplot+[only marks ] plot[error bars/.cd, y dir=both, y explicit]
% coordinates {
% (0,0)
% (345,1.28)
% (592, 2.16)
% (810, 2.90)
% (967 , 3.46)
% (1308, 4.61)
% };
% \addplot+[only marks ] plot[error bars/.cd, y dir=both, y explicit] coordinates {
% (0,0)
% (246,1.23)
% (427, 2.19)
% (621, 3.14)
% (918 , 4.52)
% (1322, 6.41)
% };
% \addplot+[only marks ] plot[error bars/.cd, y dir=both, y explicit]
% coordinates {
% (0,0)
% (332,1.99)
% (466, 3.14)
% (659, 4.30)
% (927 , 5.80)
% (1322, 7.62)
% };
%
%
% \addplot[blue, domain=0:1350]{0.042+0.0035*x};
% \addplot[red, domain=0:1350]{0.0619+0.0048*x};
% \addplot[brown, domain=0:1350]{0.0210+0.0063*x};
% \legend{
% $q_{max} = 0.2276 \frac{г}{с}$, $q_{1} = 0.123 \frac{г}{с};$ , $ q_{2} = 0.085 \frac{г}{с};$
% };
% \end{axis}
% \end{tikzpicture}
% \end{center}
% ${y_{max} = (3.5 \pm 0.3 )10^{-3}x+0.042 \pm 0.027}$,\\
% $\; {y_{1} = (4.8 \pm 0.1)10^{-3}x+0.062 \pm 0.044}$,\\
% $\; {y_{2} = (6.3 \pm 0.2)10^{-3}x+0.021 \pm 0.12}.$
% $Чтобы \; найти \; k_{max},\; k_1, \; k_2$, необходимо перевернуть соотвествующие коэффициенты $k'$ наклонов прямых, а погрешность считать по формуле: $$\sigma_k = k^2 \sigma_{k'}.$$
% $k_{max} = 0.286 \pm 0.025 \; \frac{Вт}{К}, \; k_{1} = 0.208 \pm 0.004 \; \frac{Вт}{К}, \; k_{2} = 0.159 \pm 0.005 \; \frac{Вт}{К} $
%
% Прямая для расхода $q_2$ не идёт из нуля, вероятно, из-за поспешного начала третей серии измерений (начальная температура была несколько больше комнатной). Однако в целом из вида наших прямых, можно сделать вывод о том, что тепловые потери пропорциональны разности температур.
%
%
% Построим график зависимости $k(q)$ и по его наклону определим теплоёмкость воздуха при постоянном давлении $c_p$
%
% \begin{center}
% \begin{tikzpicture}[scale = 1]
% \begin{axis}[
% axis lines = left,
% legend style={at={(1,1)}},
% xlabel = {$q, \; \frac{г}{с} \cdot 10^{-1}$},
% ylabel = {$k, \frac{Вт}{К}$
% },
% xmin=0, xmax=3,
% ymin=0, ymax=0.5,
% ymajorgrids = true,
% xmajorgrids = true,
% minor tick num = 4
% ]
% \addplot+[only marks ] plot[error bars/.cd, x dir=both, x explicit ,y dir=both, y explicit]
% coordinates {
% (2.276,0.286) +- (0.004,0.025)
% (1.23, 0.208 ) +- (0.02,0.004)
% ( 0.85,0.159) +- (0.01,0.005)
% };
% \addplot [red, domain=0:3]{0.0928+0.0860*x};
%
% \legend{
% $$y=0.0928+0.0860x$$
% };
% \end{axis}
% \end{tikzpicture}
% \end{center}
%
%Итак, из графика найдём $c_p = (0.86 \pm 0.1) \; \frac{Дж}{К \cdot г} \cdot 29 \frac{г}{моль} = 24.94 \pm 2.29 \; \frac{Дж}{К моль}.$ A также $\alpha = 0.093 \pm 0.018 \; \frac{Вт}{К}.$ К сожалению, теплоёмкость значительно отличается от теоретического значения $29.09 \; \frac{Дж}{К моль}.$ Измерения необходимо было проводить более тщательно, ждать установления термодинамического равновесия более длительное время. Однако не смотря на это, результаты, полученные при обработке данных совпадают (почти) с теоретическими с учётом погрешностей.
%
%Посчитаем доли тепловых потерь в опытах : $\frac{N_{пот}}{N} =\frac{\alpha}{k}$
%
% \begin{center}
% \begin{tabular}{|c|c|}
% \hline
% $q, \; \frac{г}{с}$ & $\frac{N_{пот}}{N}$ \\
% \hline
% 0.085 & $0.59 \pm 0.11$ \\
% \hline
% 0.123 & $0.45 \pm 0.09$\\
% \hline
% 0.2276 & $0.33 \pm 0.06$ \\
% \hline
% \end{tabular}
% \end{center}
%
%\end{enumerate}
\end{document}