Created
February 29, 2016 16:54
-
-
Save den-run-ai/37a7ab28c722e7fbc3b2 to your computer and use it in GitHub Desktop.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"cell_type": "code", | |
"execution_count": 1, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"from sympy import symbols, solve" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 2, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [], | |
"source": [ | |
"mm, mwaq, mwt, mg, c = symbols(\"mm, mwaq, mwt, mg, c\")" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 3, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"[(c*mg + c*mwt + mm - mwt - sqrt(c**2*mg**2 + 2*c**2*mg*mwt + c**2*mwt**2 + 2*c*mg*mm - 2*c*mg*mwt - 2*c*mm*mwt - 2*c*mwt**2 + mm**2 + 2*mm*mwt + mwt**2))/(2*(c - 1)),\n", | |
" (c*mg + c*mwt + mm - mwt + sqrt(c**2*mg**2 + 2*c**2*mg*mwt + c**2*mwt**2 + 2*c*mg*mm - 2*c*mg*mwt - 2*c*mm*mwt - 2*c*mwt**2 + mm**2 + 2*mm*mwt + mwt**2))/(2*(c - 1))]" | |
] | |
}, | |
"execution_count": 3, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"mwaqcalc = solve((1 - mm/(mm+mwaq)) * c - (mwt - mwaq)/(mwt - mwaq + mg),mwaq)\n", | |
"mwaqcalc" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 4, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [ | |
"from sympy import init_printing, latex\n", | |
"init_printing(latex_printer=lambda *args, **kwargs: latex(*args, mul_symbol='dot', **kwargs))" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 5, | |
"metadata": { | |
"collapsed": false | |
}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAABX4AAAApCAYAAACC/zGFAAAABHNCSVQICAgIfAhkiAAAEDpJREFU\neJzt3XmsHWUZx/HvvVfKbgsFoSyWALXspCigQKllMUiRTS1BDEgQQYJhExDigmGxICWAGEACHMIi\nIGhQNstiWSQKqIAYFRcICjQIsraICPjHM4c7nXvOubO8M+87M79PcnPvmTtn5pn3ed/3zHbeARER\nERERERERERGRGhvyHYCIiIiIiIiIiIhI2dp0InQz4BRgLvCO31BEpCbmASf6DkJEREREREREpIeV\ngMW+g/BtE+AeYBXfgYhIrVzvOwARERERERERkR62BxYAy/oOxKf3A48Cm/oORERqZTpwmu8gRERE\nRERERET6OA640HcQPl0GnOw7CBGpnWOBj/kOQkRERERERESkj2HgV8CnfAfiwzbAs8AKvgMRkdq5\nDutARURERERERERCtTPwF2BC8h9NP6nxbeB8YInvQESkVlbGBkfXgyBFREREREREJGR3AS8CB/kO\npErTgbeAtXwHIiK1sy/wOd9BiIiIiIiIiIikcDjwe99BVOl0YKHvIESkli4GJvsOQkREREREREQk\nhcnA28CH4xObPNTDnsCdvoMQkVpaFfuahIiIiIiIiIhI6F4EHgH2ik9s6onf1YHNgPt9ByIitTMD\n6yxFREREREREROriPuxBb+95n8OFrwxcDhwLPO1wuXlsHf1u1dgWIuLEHOCWAu/fFtgOmAhsD5wK\n3OsgrirUOfas2rSt0q58t2lbRUREmkCf3UblIG3nog08DhwKjGDDPjjzReAU4F1gPZcLzukE4AXf\nQYhILd1Q4L0rAfNir+cCS4C1C0VUjTrHnlWbtlXale82bauIiEgT6LPbqByk7Vy1gZnYudlpjuIa\nI5QTvxdgZ7lFRLJYDXuwW15bAO8AG0Sv34/1i3MLxlWFOseeVZu2VdqV7zZtq4iISBPos9uoHKTt\nXLWBadH7dh5vxrxCOfF7A7DQdxAiUjufB/Yt8P4h7OsZQ9HrTbF+cUbBuKoQQuzbAD8GbsOG6rmU\ncq7yh7CtUp025dv3tlbVhiUcynm7Kf8yiOpHOr4/u0MRQjmozopPrtrApOh9B7gLbWmhnPi9LfqR\ndtkAeBaY6jsQySSkvHWw8cpduRKY73B5Vao69q2ABdgHFdhXXe4Fnqf8z5U65yl0IbXvrqbkO03Z\nVrmtPtuw+KGc10cZfbHy3xyqH9Uar7ybsp9SVJuORWRUiMcOZSmrL1gOOzf7pZxxjSuUE793Azf5\nDsKTLXD70L46ORr4D1bRs2hzmYWgyrx9EBvkvJcR4PqMyxvkEOAsRq/Y1Ymr2LPk6BZgw8S0Gdjn\nyrUF4xikznmqg9D65Sble7yyrbod+2rD4o9yXh95++JBlP/mUP2o1qDy9rEPHqI2HYvI0kI7dihT\nWX3BMFZvv5I/tMFCOfHb5jt+O4SRAx9uBO7J8b4O7S2zEFSVt+WAh7GvT/SyPXBcjjh62QPrrLvr\nXc/RcqvgMvZOhve/DjwNfCAx/SXKe1hnnfNUFyH1y03L96Cy9dGOfbRh8Us5r4+8ffEgyn9zqH5U\nq195+9oHD02bjkVkrJCOHcpWVl8w5o7f4RzB1cFi3F6xlHrYAfs6htRLVXk7Dhujqd8g53OwK71F\nzQLWiJa1JrAbMMXBcqvgM/Yno3WvmJj+JrB8Ceurc57qJJR+uYn57le2vra16jYs/inn9VFGX6z8\nN4fqR7V6lXcT91PyaNOxiPQWyrFDFcrqC7r1dXF3gqtboQ/AggaYB9wPXOBo2XksAqanmG9z7GTQ\nq8AbWKM+A7vdOo29gF2ALYGDgFWBz2Jn17cDzsbuPD4GmIxdPZoAHAy81WN5E4CTsafw/Qv4E7Ar\ncAdwUcqY8gihHDbCbkWfClzF0l+nuAG7U3Ne9N4HsfKYB+wXLWeV6H+zgduBnwHfT1sAGbjKueu6\n42tb6pK3dbHOczq2zWf0mGdjrM2l1avdXAvcjI0JFTcxY7xxVfUzdwFnOo49i49i4ysvik1bC/vg\nW5hzmf36trVwm6cQPwva1L67qso3+CvfJxlctuuXsK1pldGGIYx9lK5Q2qbr7cqryn677fkuo78o\nqsn5r6p9qX5kE0r9cF3et+HvsxvCaTc+92HAbZ0NpUyrjrkuxw5N7wtWj34vGjhXAxyP7aANMhcb\nSHnL6PWOwGtYJUhjAnBu9PdD2Mnu4xgdg+NEbCDw+Yzemj0SreOgHstbHutQrotNOwn4H7BPypgg\n+y3uIZTDCHAxdiHiKOCx2PInAm/H5p2MVeDkUB6Hkf9qXId0ZeYq567rTh5tylvcgdjOYXJ96wDn\nZ1hO0XaTVqj9TFodin3lZh5Wj7bL8d4256iN7buqfEMY5VukbLPqkL8dF2nDEMY+SlcobdP1drnm\ns99uar7r1F80If9VtS/Vj2xCqR91Ke8O7Tu+7aVD9ccidSzTEPZtoV7ndELuC2ZhJ7GndSfUeaiH\nyfR/SNPjwGrYmfteZmBPyPsa8Gg0bRLwBPCblOuchVWSIezq1CKscrwb/f+taP3XAE9F096Oftbo\nsewzgU2AQ2PTfofl6L4BMRWRpxyS5e6iHOYAP8V2pj8J/Dm2/B2wMujeAv8icCxjx9mZjTXcN/pu\nbXGucu667vjalrrkLe4uYFlsPN+43Uk/Lnje/iONMtpXnI9+Jq8NgSOB7wAPDJiv12dBmTlKrjOE\nHLW5X4Zq8w1hlG/VfWceadswuGvHZfahobRN19vlUtX9dmj5TsbkKp669BdNyX+Z7Uv1I99nQkj1\noy7lnVYd2o0veetsHco0xH1bqNc5nZD7gg8BS4C/l7DsSm2LJeSmPv9fDUvcjn3+fwc2cHe/E8dp\n1jkFOzu/RbSumYn5f8jYDmL9aN45ienrRss+NzH9NOAPGWKEbFe1spZDr3J3UQ7rYldK1sUa1B6x\nec+MYozbjrFPKFwEnJ5yO5I6pCszVzl3WXfyalPekv6IXb2Nu4b044Ln6T/SKKt9dbnsZ9LqkC9H\ny2JfhzlnnPn6fRaUlaNe6/Sdo7b3y1BtviGM8i1Stll1yN6O07ZhcNeOy+xDQ2qbLpfjUtX9dmj5\n7hWTy3hC7y+alP+y2pfqR/7PhJDqRx3KG9p5fNtLh3KPRaC8tl1VX+RqXXU5dmhDX3AudnL7Pa7G\n+K3a68DL9B+z4gXsVutegyWvjj3c6UIsQXnX+Vz0eyfsLP2vE/N/HLg8MW037Gvm9ySmfxrLRfKO\nw1n0H9j6Cka/6hL3QWAb4L89/ncIo1dF85RDr3J3UQ7/iH4fALwC/Dw2746MLa+tE/Nsil19Sc6X\nVLTMXOXc1XI2B37A6NcSxvMIcLjDGOqSt6Q7WfoBbxOwK3JpxgfL23+kUVb76srTz6TlMkdD2Dbd\nDnxznPX2KrMyc9Rrnb5z1PZ+uep8g//yTVu2Wblqx1naMLhrx2X2oSG1TZfLKbIfEeej3w4t371i\nchVPWf2F8l9u+xov5tDrB7ipI0U/E0KrH6F9fjft+DYvX8ciUF7brqovcrWuuhw7NLUviNuBPg+t\nfzfgn7xOxb7WnbRttNyDCyw77qYe69k4WscnEtPvAH4S/T2V0auWF0XzT4rNuxxWufbPGE+HdFc6\nQiyHu4EfxeYbwcY9OTLx/ksSr4/Erlp1B8GehF15SatDtiuBLrbV5XKKaFPeuvbGdhxXiV7vChyR\n8r2u201aofUzaXXInqPTgG8kph2Y4f3KkduYQm/fvvIN/sq3aNlm1SFbOy7ahiG8fZQQ26bL5RQV\nWr/d1HyH2l80Lf+ul1NVvKof6TS9vLs6tPf4Nq5D9cciXXUs0zYcO3Q1tS9YFTvXsVV84nD0eyjg\nn7yuxG7bnpKY/nz0+7Ue75mKjc+R1jB2Nn9hYvpsLJG/jE1bNZp+dfT6REavWr4cxfNybP5Z2FcM\nyroKEGI5rMnS45Csid2NGZ82HfhrYj0zsSvMr0evj8LGWimDq211tZwi2pS3uF9gHfrs6PUc4NaU\n73XZbtKqcz+T1cHAO9iFu7gdMixDOXIbU+jt20e+wW/5+uo703DRhiG8fZQQ22YI+xEQXr/d1HxD\nmP1FE/PvcjlVxQuqH2k0vbzzqlu7KZOr/Zg6lmlbjh2g2X3BftiwVL+NTxzuPW8jPIF9pTt5deZJ\nYAFjx/+dCXyLbA84moGdnV+YmD4bG6h5cWzaetjZ/wXYbd7xMcJuAFYAPhC9nop9neVv2NNLyxBi\nOfwGG0S96wDg38A60esRbFDsCxLrGWF04O2tsVv2n6McrrbV1XKKaFPe4l4BHgZ2iV6vF4tjPC7b\nTVp17mey2An4LjaG0lWxn2uBlTMsRzlyG1Po7dtHvsFv+frqO8fjqg1DePsoIbbNEPYjQuy3m5pv\nCK+/aGr+XS6nqnhB9SONJpd3EXVrN2VxuR9TxzJty7EDNLsvOBg4b9AM2wLHAKdgtzL3ezCaa2Wu\n9yPYDtoKiemTsHE7LsUa93nYuB9ZT4TvCTyOndGPWwgclpg2gt0CfhG2rcl1HQrcDJwV/f9h4LKM\n8UC2W9xDK4fVsIdsXR7Fsjc2CPZd2ADV3wM26rH+LbGrMvOBE3LE3yF9mbnaVpd1J6825S3pdOzJ\nm9Owup+Fq3aTVoj9TFod0ufoJfoP+ZO86j4e5ahd7bvqfIPf8i1atll1SJcLl20YwttHCa1thrAf\nEWK/3dR8Q3j9RZPzX1X7Uv1IJ7T6AeGVd1KHdh7fJnXwcyxSxzJt07FDU/uCXbEbYJfpN8NKLP2U\n+7nAEmBtRwH4XO8lwNcdLq8KK2CDUn8hx3s75D8Z1lYdVGZ11CF/3mZjH+JnYwO4t1GRfiatDmpb\nRVSRo1B1UN0JRQflIqnNbbONlG8RqZsO+uwGlUObdGh3roeBB7FhLHv+E+x29hOADaLXtwPLA9uX\nHFwV6z0W+Az2VNIQrY5dcYjbHRvfeEGO5b2KPXxC0lOZ1VORvD2AfaXiQOB+ZxGFy3U/k5baVnq+\nchQq1Z1wtD0XapvtonyLSBO0/bO7S+XQHm3P9VexE7+3DJppCBtyYSh6vSl2N9yMcRY+jJ1YPRob\nS+LLGYPLu96sNsJu217F8XJduBa7i6A7HMUUbAyxY7xFJNIOC7Bx/NpA/Uz4lCORMKlttovyLSIi\nInUyE7uJNjl0xXtGYn8/E/v7HOzBaFck5t8Te9DB09Hri7En3s7HHnxwGDa2RRZp1lvUC9jZ7/nA\njdjJ5VAsgw0WPh17Kun+2PijV/oMSqQFpmCDtz/mO5AKqJ8Jn3IkEia1zXZRvkVERKQupgEnYfsr\nfe94Huox7RBsZ+dExp4gvRr4Z/S/rYFbsafVvQmsGP08DxyBDePQz0PAdRnWKyIiIiIiIiIiIiI5\n7YGdgAVYjsGDIx+PPd2u6vWKiIiIiIiIiIiIyADDsb9nAWtggwGvCeyGfRW6n2eAxYlp++SIIet6\nRURERERERERERGSA7lAP6wOPAisl/j8Rezpev/eeATwFLMEegnAzS4/ZO5486xURERERERERERER\nEREREREREREREREREREREREREREREREREREREREREREREREREREREREREXHu/8odd1MDBjc4AAAA\nAElFTkSuQmCC\n", | |
"text/latex": [ | |
"$$\\frac{1}{2 \\cdot c - 2} \\cdot \\left(c \\cdot mg + c \\cdot mwt + mm - mwt + \\sqrt{c^{2} \\cdot mg^{2} + 2 \\cdot c^{2} \\cdot mg \\cdot mwt + c^{2} \\cdot mwt^{2} + 2 \\cdot c \\cdot mg \\cdot mm - 2 \\cdot c \\cdot mg \\cdot mwt - 2 \\cdot c \\cdot mm \\cdot mwt - 2 \\cdot c \\cdot mwt^{2} + mm^{2} + 2 \\cdot mm \\cdot mwt + mwt^{2}}\\right)$$" | |
], | |
"text/plain": [ | |
" _________________________________________________\n", | |
" ╱ 2 2 2 2 2 \n", | |
"c⋅mg + c⋅mwt + mm - mwt + ╲╱ c ⋅mg + 2⋅c ⋅mg⋅mwt + c ⋅mwt + 2⋅c⋅mg⋅mm - 2⋅c\n", | |
"──────────────────────────────────────────────────────────────────────────────\n", | |
" 2⋅(c - 1) \n", | |
"\n", | |
"________________________________________________________\n", | |
" 2 2 2 \n", | |
"⋅mg⋅mwt - 2⋅c⋅mm⋅mwt - 2⋅c⋅mwt + mm + 2⋅mm⋅mwt + mwt \n", | |
"────────────────────────────────────────────────────────\n", | |
" " | |
] | |
}, | |
"execution_count": 5, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"mwaqcalc[1]" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": null, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"display_name": "Python 2", | |
"language": "python", | |
"name": "python2" | |
}, | |
"language_info": { | |
"codemirror_mode": { | |
"name": "ipython", | |
"version": 2 | |
}, | |
"file_extension": ".py", | |
"mimetype": "text/x-python", | |
"name": "python", | |
"nbconvert_exporter": "python", | |
"pygments_lexer": "ipython2", | |
"version": "2.7.9" | |
}, | |
"latex_envs": { | |
"bibliofile": "biblio.bib", | |
"cite_by": "apalike", | |
"current_citInitial": 1, | |
"eqLabelWithNumbers": true, | |
"eqNumInitial": 0 | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 0 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment