diff --git a/04_Machine_Asynchrone_1.ipynb b/04_Machine_Asynchrone_1.ipynb
new file mode 100644
index 0000000..5f92c3a
--- /dev/null
+++ b/04_Machine_Asynchrone_1.ipynb
@@ -0,0 +1,194 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "id": "eea6c81b-cb1e-4f0e-83bf-3d7f71fa04d0",
+ "metadata": {},
+ "source": [
+ "
MACHINE ASYNCHRONE 1
"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "id": "1108e91d-d204-4a7e-98e1-c49e972c6efb",
+ "metadata": {},
+ "source": [
+ "1. Donnée
\n",
+ "\n",
+ "Une machine asynchrone triphasée couplée en étoile a les caractéristiques suivantes :\n",
+ "\n",
+ "\n",
+ " - Tension de ligne : Uligne = 380 V
\n",
+ " - Fréquence : f = 50 Hz
\n",
+ " - Couple utile : Tutile = 95 Nm
\n",
+ " - Nombre de paires pôles : p = 3
\n",
+ " - Facteur de puissance : cosφ = 0.75
\n",
+ " - Vitesse : N = 950 tr/min
\n",
+ " - Pertes par frottement et ventilation : Pfv = 400 W
\n",
+ " - Pertes fer : Pfer = 200 W
\n",
+ "
\n",
+ "\n",
+ "Sachant que les pertes joule statoriques sont égales aux pertes joule rotoriques (Pjs = Pjr), déterminer, pour le fonctionnement sous les conditions données ci-dessus :\n",
+ "\n",
+ "\n",
+ " - Le rendement de la machine
\n",
+ " - Le courant de ligne
\n",
+ "
"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "id": "b43109b2-aa0b-46bd-8a8c-df13481de04a",
+ "metadata": {},
+ "source": [
+ "2. Préambule
\n",
+ "\n",
+ "Le but de cet exercice est de comprendre les différentes puissances intervenant dans l’étude de la machine asynchrone et en particulier la puissance d’entrefer Pδ, la puissances mécanique Pmec et la puissance utile Putile. "
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "id": "4784c9a0-6ec1-4500-8882-5ed47ee1b015",
+ "metadata": {
+ "tags": []
+ },
+ "source": [
+ "3. Aide
\n",
+ "1) Bilan de puissance - Moteur
\n",
+ "\n",
+ "2) Puissance d'entrefer
\n",
+ "\n",
+ "3) Couple électromagnétique
\n",
+ "\n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 3,
+ "id": "ddf9e91e-253c-4ca4-a36b-e7ac24b8b697",
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Unligne = 380\n",
+ "f = 50\n",
+ "Tu = 95\n",
+ "p = 3\n",
+ "cosphi = 0.75000\n",
+ "Nn = 950\n",
+ "Pfv = 400\n",
+ "Pfer = 200\n",
+ "Nm = 950\n",
+ "Omega_m = 99.484\n",
+ "Pu = 9451.0\n",
+ "Pmec = 9851.0\n",
+ "Omega_s = 104.72\n",
+ "s = 0.050000\n",
+ "Ns = 1000\n",
+ "s = 0.050000\n",
+ "Pdelta = 10369.42937\n",
+ "Pjr = 518.47\n",
+ "Pjs = 518.47\n",
+ "Pel = 11087.90084\n",
+ "rend = 0.85237\n"
+ ]
+ }
+ ],
+ "source": [
+ "Unligne = 380 % V\n",
+ "f = 50 % Hz\n",
+ "Tu = 95 % Nm\n",
+ "p = 3 % nombre de paires de poles\n",
+ "cosphi = 0.75 % facteur de puissance\n",
+ "Nn = 950 % tr/min\n",
+ "Pfv = 400 % W\n",
+ "Pfer = 200 % W\n",
+ "\n",
+ "% 1. Le rendement de la machine\n",
+ "Nm = \n",
+ "Omega_m = \n",
+ "Pu = \n",
+ "Pmec = \n",
+ "\n",
+ "% s (calcul avec Omega_s)\n",
+ "Omega_s = \n",
+ "s = \n",
+ "\n",
+ "% s (calcul avec Ns)\n",
+ "Ns = \n",
+ "s = \n",
+ "\n",
+ "Pdelta = \n",
+ "Pjr = \n",
+ "Pjs = \n",
+ "Pel = \n",
+ "\n",
+ "rend = "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "id": "3c96b928-7555-48cf-a029-5ab23dbfbcb3",
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Uph = 219.39\n",
+ "Iph = 22.462\n",
+ "Iligne = 22.462\n"
+ ]
+ }
+ ],
+ "source": [
+ "% 2. Le courant de ligne\n",
+ "% Pel = 3*Uph*Iph*cosphi\n",
+ "Uph = \n",
+ "Iph = \n",
+ "\n",
+ "Iligne = "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "id": "a3bfaf9e-489f-4622-83a7-8433390981c6",
+ "metadata": {},
+ "outputs": [],
+ "source": []
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Octave",
+ "language": "octave",
+ "name": "octave"
+ },
+ "language_info": {
+ "file_extension": ".m",
+ "help_links": [
+ {
+ "text": "GNU Octave",
+ "url": "https://www.gnu.org/software/octave/support.html"
+ },
+ {
+ "text": "Octave Kernel",
+ "url": "https://github.com/Calysto/octave_kernel"
+ },
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://metakernel.readthedocs.io/en/latest/source/README.html"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "octave",
+ "version": "5.2.0"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}
diff --git a/04_Machine_Asynchrone_2.ipynb b/04_Machine_Asynchrone_2.ipynb
new file mode 100644
index 0000000..a240bd4
--- /dev/null
+++ b/04_Machine_Asynchrone_2.ipynb
@@ -0,0 +1,249 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "id": "f9875f2f-8874-4eb1-818b-75d999200069",
+ "metadata": {},
+ "source": [
+ "MACHINE ASYNCHRONE 2
"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "id": "078a4de4-b208-4628-8c27-3f1fd77016a6",
+ "metadata": {},
+ "source": [
+ "1. Donnée
\n",
+ "\n",
+ "Un moteur asynchrone à cage a les caractéristiques suivantes :\n",
+ "\n",
+ "\n",
+ " - Sn = 34.5 kVA
\n",
+ " - Unligne = 400 V
\n",
+ " - Stator couplé en Y
\n",
+ " - fn = 50 Hz
\n",
+ " - p = 3 (nombre de paires de pôles)
\n",
+ " - rs = 0.05 pu
\n",
+ " - xos = 0.10 pu
\n",
+ " - xh = 10 pu
\n",
+ " - r'r = 0.04 pu
\n",
+ " - x'or = 0.15 pu
\n",
+ " - Pfer = 0 (Rfer est négligé dans le schéma équivalent)
\n",
+ "
\n",
+ "\n",
+ "La machine est alimentée sous tension et fréquence nominales.\n",
+ "\n",
+ "\n",
+ " - Pour un glissement s = 0.04, calculer :
\n",
+ " \n",
+ " - Le courant de ligne au stator
\n",
+ " - La puissance électrique active consommée au stator
\n",
+ "
\n",
+ " - Calculer le glissement critique sK et le couple de décrochage TK
\n",
+ "
"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "id": "90a8c70d-5cdc-49b7-a628-1ffdac157e1a",
+ "metadata": {},
+ "source": [
+ "2. Préambule
\n",
+ "\n",
+ "Le but de cet exercice est multiple :\n",
+ "\n",
+ "\n",
+ " - Le passage des pu aux vraies grandeurs pour les paramètres du schéma équivalent.
\n",
+ " - L’application de l’équivalent de Thévenin pour le calcul du couple et du glissement critique (ceci pourrait évidemment être étendu au calcul du couple).
\n",
+ " - Pour finir, l’étude du cas de freinage permet de valider la compréhension du glissement et des 3 modes de fonctionnement (moteur, générateur et frein).
\n",
+ "
"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "id": "5840eecd-fc14-4337-a20d-f83ed43a5e32",
+ "metadata": {},
+ "source": [
+ "3. Aide
\n",
+ "1) Impédence équivalente
\n",
+ "\n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "id": "dfe7e73f-61af-485b-b849-b76c95103e39",
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Sn = 34500\n",
+ "Unligne = 400\n",
+ "Fn = 50\n",
+ "p = 3\n",
+ "rs_pu = 0.050000\n",
+ "xss_pu = 0.10000\n",
+ "xh_pu = 10\n",
+ "rrp_pu = 0.040000\n",
+ "xsrp_pu = 0.15000\n",
+ "s = 0.040000\n",
+ "Un = 230.94\n",
+ "In = 49.796\n",
+ "Zn = 4.6377\n",
+ "Rs = 0.23188\n",
+ "Xss = 0.46377\n",
+ "Xh = 46.377\n",
+ "Rrp = 0.18551\n",
+ "Xsrp = 0.69565\n"
+ ]
+ }
+ ],
+ "source": [
+ "Sn = 34.5e3\n",
+ "Unligne = 400\n",
+ "Fn = 50\n",
+ "p = 3\n",
+ "rs_pu = 0.05\n",
+ "xss_pu = 0.1\n",
+ "xh_pu = 10\n",
+ "rrp_pu = 0.04\n",
+ "xsrp_pu = 0.15\n",
+ "\n",
+ "s = 0.04\n",
+ "\n",
+ "% Vraies grandeurs\n",
+ "Un = \n",
+ "In = \n",
+ "Zn = \n",
+ "\n",
+ "% Calcul des vraies grandeurs\n",
+ "Rs = \n",
+ "Xss = \n",
+ "Xh = \n",
+ "Rrp = \n",
+ "Xsrp = "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 13,
+ "id": "154c72ab-fbb3-4ee8-8469-a88f1acd34a3",
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Zeq_cmpx = 4.6902 + 1.5884i\n",
+ "Zeq = 4.9519\n",
+ "Uph = 230.94\n",
+ "Iph = 46.637\n",
+ "cosphi = 0.94716\n",
+ "Pel_a = 30603.56849\n",
+ "Uph_cmpx = 230.94\n",
+ "Iph_cmpx = 44.172 - 14.959i\n",
+ "Iph = 46.637\n",
+ "S = 30603.56849 + 10364.15709i\n",
+ "Pel_b = 30603.56849\n"
+ ]
+ }
+ ],
+ "source": [
+ "% 1\n",
+ "Zeq_cmpx = \n",
+ "Zeq = \n",
+ "\n",
+ "% 1.a \n",
+ "Uph = \n",
+ "Iph = \n",
+ "cosphi = \n",
+ "\n",
+ "Pel_a = \n",
+ "\n",
+ "% 1.b \n",
+ "Uph_cmpx = \n",
+ "Iph_cmpx = \n",
+ "Iph = \n",
+ "S = \n",
+ "\n",
+ "Pel_b = "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 14,
+ "id": "0a6fea5c-5193-42cf-9b92-aa2b38371347",
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Uph_cmpx = 230.94\n",
+ "Ue_cmpx = 228.6480 + 1.1319i\n",
+ "Ue = 228.65\n",
+ "Ze_cmpx = 0.22731 + 0.46030i\n",
+ "Re = 0.22731\n",
+ "Xe = 0.46030\n",
+ "sk = 0.15746\n",
+ "Omega_s = 104.72\n",
+ "Tk = 532.85\n"
+ ]
+ }
+ ],
+ "source": [
+ "% 2\n",
+ "Uph_cmpx = \n",
+ "Ue_cmpx = \n",
+ "Ue = \n",
+ "Ze_cmpx = \n",
+ "Re = \n",
+ "Xe = \n",
+ "\n",
+ "sk = \n",
+ "Omega_s = \n",
+ "\n",
+ "Tk = "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "id": "1a8e6a4b-79be-4139-bccb-21211e34b479",
+ "metadata": {},
+ "outputs": [],
+ "source": []
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Octave",
+ "language": "octave",
+ "name": "octave"
+ },
+ "language_info": {
+ "file_extension": ".m",
+ "help_links": [
+ {
+ "text": "GNU Octave",
+ "url": "https://www.gnu.org/software/octave/support.html"
+ },
+ {
+ "text": "Octave Kernel",
+ "url": "https://github.com/Calysto/octave_kernel"
+ },
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://metakernel.readthedocs.io/en/latest/source/README.html"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "octave",
+ "version": "5.2.0"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}
diff --git a/04_Machine_Asynchrone_3.ipynb b/04_Machine_Asynchrone_3.ipynb
new file mode 100644
index 0000000..1ad12da
--- /dev/null
+++ b/04_Machine_Asynchrone_3.ipynb
@@ -0,0 +1,331 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "id": "f9875f2f-8874-4eb1-818b-75d999200069",
+ "metadata": {},
+ "source": [
+ "MACHINE ASYNCHRONE 3
"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "id": "078a4de4-b208-4628-8c27-3f1fd77016a6",
+ "metadata": {},
+ "source": [
+ "1. Donnée
\n",
+ "\n",
+ "Un moteur à cage a les caractéristiques suivantes :\n",
+ "\n",
+ "\n",
+ " - 4 pôles
\n",
+ " - Unligne = 3 kV
\n",
+ " - fn = 50 Hz
\n",
+ " - Putile = 1 MW
\n",
+ " - Stator couplé en triangle
\n",
+ " - Rs = 0.2 Ω
\n",
+ " - Xos = 2.3 Ω
\n",
+ " - Xh = 80 Ω
\n",
+ " - R'r = 0.5 Ω
\n",
+ " - X'or = 2.4 Ω
\n",
+ " - Pfv = 6kW
\n",
+ "
\n",
+ "\n",
+ "\n",
+ " - Déterminer le couple électromagnétique et le glissement quand le moteur fonctionne sous conditions nominales.
\n",
+ " \n",
+ " - A partir de l'équation du couple (complète).
\n",
+ " - A partir de l'équation du couple simplifiée.
\n",
+ "
\n",
+ " - Quels sont, dans ce cas, le courant absorbé par le moteur ainsi que le facteur de puissance ?
\n",
+ "
"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "id": "90a8c70d-5cdc-49b7-a628-1ffdac157e1a",
+ "metadata": {},
+ "source": [
+ "2. Préambule
\n",
+ "\n",
+ "Cet exercice met en oeuvre la comparaison de l’équation de couple complète comparée à l’équation simplifiée. L’exercice passe par un peu de mathématique où il faudra repartir des équations de couple et les exprimer à partir de la puissance électromagnétique pour en ressortir un polynôme en s.\n",
+ "\n",
+ "Finalement, la deuxième partie est une mise en pratique de notions vues dans les exercices MAS1 et MAS2."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "id": "5840eecd-fc14-4337-a20d-f83ed43a5e32",
+ "metadata": {},
+ "source": [
+ "3. Aide
\n",
+ "1) Bilan de puissance - Moteur
\n",
+ "\n",
+ "2) Equivalent de Thévenin
\n",
+ "\n",
+ "3) Couple électromagnétique
\n",
+ "\n",
+ "4) Caractéristique de couple
\n",
+ "\n",
+ "5) Couple et glissement critiques
\n",
+ "\n",
+ "6) Impédence équivalente
\n",
+ ""
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "id": "dfe7e73f-61af-485b-b849-b76c95103e39",
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "p = 2\n",
+ "Unligne = 3000\n",
+ "f = 50\n",
+ "Pu = 1000000\n",
+ "Rs = 0.20000\n",
+ "Xss = 2.3000\n",
+ "Xh = 80\n",
+ "Rrp = 0.50000\n",
+ "Xsrp = 2.4000\n",
+ "Pfv = 6000\n"
+ ]
+ }
+ ],
+ "source": [
+ "p = 2 % 4 pôles -> 2 paires\n",
+ "Unligne = 3e3 % V\n",
+ "f = 50 % Hz\n",
+ "Pu = 1e6\n",
+ "Rs = 0.2 % Ohm\n",
+ "Xss = 2.3 % Ohm\n",
+ "Xh = 80 % Ohm\n",
+ "Rrp = 0.5 % Ohm\n",
+ "Xsrp = 2.4 % Ohm\n",
+ "Pfv = 6e3 % W"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 5,
+ "id": "154c72ab-fbb3-4ee8-8469-a88f1acd34a3",
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> Uph = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> Ue_cmpx =\n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> Ue = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> Ze_cmpx = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> Re =\n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> Xe = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> Pmec = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> a = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> b =\n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> c = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> s1a = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> s2a = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> sk = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> Omega_m = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> Nm = \n",
+ " ^\n",
+ "\n",
+ "parse error:\n",
+ "\n",
+ " syntax error\n",
+ "\n",
+ ">>> Tem1 = \n",
+ " ^\n",
+ "\n"
+ ]
+ }
+ ],
+ "source": [
+ "% 1. Couple éléctromagnétique et glissement\n",
+ "\n",
+ "%1.a\n",
+ "Uph = \n",
+ "Ue_cmpx =\n",
+ "Ue = \n",
+ "Ze_cmpx = \n",
+ "Re =\n",
+ "Xe = \n",
+ "\n",
+ "Pmec = \n",
+ "\n",
+ "a = \n",
+ "b =\n",
+ "c = \n",
+ "\n",
+ "s1a = \n",
+ "s2a = \n",
+ "\n",
+ "sk = \n",
+ "\n",
+ "Omega_m = \n",
+ "Nm = \n",
+ "Tem1 = \n",
+ "\n",
+ "%1.b\n",
+ "a = \n",
+ "b = \n",
+ "c = \n",
+ "\n",
+ "s1b = \n",
+ "s2b = \n",
+ "\n",
+ "Omega_m_b = \n",
+ "Nm_b = \n",
+ "Tem2 = "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 4,
+ "id": "0a6fea5c-5193-42cf-9b92-aa2b38371347",
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "% 2\n",
+ "\n",
+ "s =\n",
+ "Zeq_cmpx = \n",
+ "Zeq =\n",
+ "Iph =\n",
+ "Iligne = \n",
+ "cosphi = "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "id": "1a8e6a4b-79be-4139-bccb-21211e34b479",
+ "metadata": {},
+ "outputs": [],
+ "source": []
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Octave",
+ "language": "octave",
+ "name": "octave"
+ },
+ "language_info": {
+ "file_extension": ".m",
+ "help_links": [
+ {
+ "text": "GNU Octave",
+ "url": "https://www.gnu.org/software/octave/support.html"
+ },
+ {
+ "text": "Octave Kernel",
+ "url": "https://github.com/Calysto/octave_kernel"
+ },
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://metakernel.readthedocs.io/en/latest/source/README.html"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "octave",
+ "version": "5.2.0"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}
diff --git a/MAS1_1.png b/MAS1_1.png
new file mode 100644
index 0000000..bba8ebd
Binary files /dev/null and b/MAS1_1.png differ
diff --git a/MAS1_2.png b/MAS1_2.png
new file mode 100644
index 0000000..02081bd
Binary files /dev/null and b/MAS1_2.png differ
diff --git a/MAS1_3.png b/MAS1_3.png
new file mode 100644
index 0000000..d8ba0bf
Binary files /dev/null and b/MAS1_3.png differ
diff --git a/MAS2_1.png b/MAS2_1.png
new file mode 100644
index 0000000..185d7cb
Binary files /dev/null and b/MAS2_1.png differ
diff --git a/MAS3_1.png b/MAS3_1.png
new file mode 100644
index 0000000..33b80cc
Binary files /dev/null and b/MAS3_1.png differ
diff --git a/MAS3_2.png b/MAS3_2.png
new file mode 100644
index 0000000..71430ca
Binary files /dev/null and b/MAS3_2.png differ
diff --git a/MAS3_3.png b/MAS3_3.png
new file mode 100644
index 0000000..44a9dd1
Binary files /dev/null and b/MAS3_3.png differ
diff --git a/octave-workspace b/octave-workspace
new file mode 100644
index 0000000..6015177
Binary files /dev/null and b/octave-workspace differ