Skip to content

Instantly share code, notes, and snippets.

@den-run-ai
Created February 29, 2016 16:54
Show Gist options
  • Save den-run-ai/37a7ab28c722e7fbc3b2 to your computer and use it in GitHub Desktop.
Save den-run-ai/37a7ab28c722e7fbc3b2 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"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