84 lines
13 KiB
Plaintext
84 lines
13 KiB
Plaintext
|
{
|
||
|
|
"cells": [
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 1,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"text/plain": [
|
||
|
|
"[<matplotlib.lines.Line2D at 0x7fcea17426d0>]"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"execution_count": 1,
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "execute_result"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAD4CAYAAADhNOGaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAfhElEQVR4nO3deXxU9b3/8ddnlmxkIytLiCwJCAIiRkDEKnVD7C3VbuJWu/wsVbpv/n699z7aX/u7tre3rbVXa22v93axbtVWWsG1rbaiQkCWIAIBBcISCEsgQNb5/v6YgcYQyCTM5Mxk3s/HYx5z5sz3TN4e5zFvzjkz55hzDhERSV0+rwOIiIi3VAQiIilORSAikuJUBCIiKU5FICKS4gJe/eGioiI3cuRIr/68iEhSWrFiRYNzrjiWr+lZEYwcOZLq6mqv/ryISFIys62xfk3tGhIRSXEqAhGRFKciEBFJcSoCEZEUpyIQEUlxPRaBmT1oZnvMrOYUz5uZ3WNmtWa2xsymxj6miIjESzRbBP8DzDnN81cDlZHbbcBPzzyWiIj0lx6LwDn3MrD/NEPmAb9yYa8B+WY2NFYBu9pYf5jv/OlNmts64vUnRERSSiyOEQwHtnd6XBeZdxIzu83Mqs2seu/evX36Y3UHjvKLv7/Nyq0H+rS8iIi8WyyKwLqZ1+3VbpxzDzjnqpxzVcXFffuF9LRRhQR8xt9rG/q0vIiIvFssiqAOGNHpcRmwMwav263s9ADnlefziopARCQmYlEEi4BbIt8emgE0Oud2xeB1T+miiiLW7Gik8WhbPP+MiEhKiObrow8DrwLjzKzOzD5pZgvMbEFkyGJgC1AL/By4PW5pIy6qKMI5eHWLtgpERM5Uj2cfdc7N7+F5B9wRs0RRmDIin0Fpfl6p3ceciXH7gpKISEpIyl8WB/0+po8u1HECEZEYSMoigPDuoS0NR9hx8JjXUUREklrSFsGsiiIAbRWIiJyhpC2CsaXZFGWnqwhERM5Q0haBmTGrInycIHy8WkRE+iJpiwDCxwkamlrZUH/Y6ygiIkkr6YsA4O+btHtIRKSvkroIhuVnMrpoEEs37/M6iohI0krqIoDwVsFrW/bR1hHyOoqISFIaEEVwtLWDVdsPeh1FRCQpJX0RXDi6EJ/pOIGISF8lfRHkZQWZVKbTUouI9FXSFwHArIpC3th+kMPNOi21iEhvDYgiuGhMER0hx7K3T3dpZRER6c6AKIKpZw0mPeDT5StFRPpgQBRBRtDPtFEFLK3V7wlERHprQBQBhL9GuqH+MHsON3sdRUQkqQyYIjh+WmptFYiI9M6AKYIJQ3PJzwrqOIGISC8NmCLw+YyLxhTptNQiIr00YIoAwscJdjU2s6XhiNdRRESSxgArgkJAl68UEemNAVUE5QVZlA3OVBGIiPTCgCqC8OUri1i6eR8dIR0nEBGJxoAqAggfJzjc3M7aHY1eRxERSQoDrghmjtFxAhGR3hhwRVCYnc6Eobm6PoGISJQGXBEAzKosYsXWAxxr7fA6iohIwhuQRTBzTCGtHSGWv6PTUouI9GRAFsG0UQUE/abjBCIiURiQRZCVFmBq+WBe2awiEBHpSVRFYGZzzGyDmdWa2Z3dPJ9nZn80s9Vmts7MPh77qL0zq6KIdTsPsf9Iq9dRREQSWo9FYGZ+4F7gamACMN/MJnQZdgfwpnPuXOBS4AdmlhbjrL1yUWURzsGrm3VaahGR04lmi2AaUOuc2+KcawUeAeZ1GeOAHDMzIBvYD7THNGkvTR6eR056gL9t2utlDBGRhBdNEQwHtnd6XBeZ19l/AuOBncBa4PPOuVDXFzKz28ys2syq9+6N7wd0wO/jPeOKeWF9vU43ISJyGtEUgXUzr+sn61XAKmAYMAX4TzPLPWkh5x5wzlU556qKi4t7GbX3rpk0lIamVl5/W7uHREROJZoiqANGdHpcRvhf/p19HHjShdUCbwNnxyZi380eV0Jm0M/itbu8jiIikrCiKYLlQKWZjYocAL4eWNRlzDbgMgAzKwXGAVtiGbQvMtP8zD67mGdqtHtIRORUeiwC51w7sBB4FlgPPOacW2dmC8xsQWTYt4GZZrYWeBH4unMuIb7EP3fSUBqaWvQrYxGRUwhEM8g5txhY3GXe/Z2mdwJXxjZabLz37BIygj4Wr93FjNGFXscREUk4A/KXxZ1lpQWYPa6EJTW7tXtIRKQbA74IILx7aO/hFqq1e0hE5CQpUQTvPbuE9IBP3x4SEelGShTBoPQAl44rZknNbkLaPSQi8i4pUQQQ3j2053ALK7Yd8DqKiEhCSZkiuGx8KWkBH0+v0e4hEZHOUqYIstMDXDq2mCU1u7R7SESkk5QpAoBrJg+l/lALK7V7SETkhJQqgveeXRLePaRvD4mInJBSRZCTEeQ9lcUsWatvD4mIHJdSRQBwzeQh7D7UzBvbD3odRUQkIaRcEVw2vpQ0v35cJiJyXMoVQW5GkPeMLWLJWn17SEQEUrAIIPzjsp2NzayqO+h1FBERz6VkEVw2vpSg31isH5eJiKRmEeRlBrm4MnzuIee0e0hEUltKFgGEdw/tOHiM1XWNXkcREfFUyhbBFRMiu4f07SERSXEpWwR5mUFmVRTx9Jpd2j0kIiktZYsA/rF7aI12D4lICkvpIrhiQikBn3YPiUhqS+kiyM9K46KKIhbXaPeQiKSulC4CgGsmDWX7/mPU7DjkdRQREU+kfBFceU5495BOTS0iqSrliyA/K42ZFUUsXqvdQyKSmlK+CACumTSEbfuP6tTUIpKSVATANZOHMSjNz29e2+p1FBGRfqciIHxh++umlvGnNbvYf6TV6zgiIv1KRRBx04yzaG0P8Xj1dq+jiIj0KxVBxLghOUwbVcBvXt+qC9aISEqJqgjMbI6ZbTCzWjO78xRjLjWzVWa2zsxeim3M/nHzjLPYvv8YL23a63UUEZF+02MRmJkfuBe4GpgAzDezCV3G5AP3Ae93zp0DfDj2UePvqnOGUJSdzm9e1UFjEUkd0WwRTANqnXNbnHOtwCPAvC5jbgCedM5tA3DO7YltzP6RFvAxf9oI/rxhD9v3H/U6johIv4imCIYDnY+g1kXmdTYWGGxmfzWzFWZ2S3cvZGa3mVm1mVXv3ZuYu1/mTyvHgN8u2+Z1FBGRfhFNEVg387oeTQ0A5wPXAFcB/2JmY09ayLkHnHNVzrmq4uLiXoftD8PyM7l8fCmPLt9OS3uH13FEROIumiKoA0Z0elwG7OxmzDPOuSPOuQbgZeDc2ETsfzdfeBb7j7SyZO1ur6OIiMRdNEWwHKg0s1FmlgZcDyzqMuYp4GIzC5hZFjAdWB/bqP3nojFFjCoaxK/1S2MRSQE9FoFzrh1YCDxL+MP9MefcOjNbYGYLImPWA88Aa4BlwC+cczXxix1fPp9x4/RyVmw9wLqdunqZiAxs5tUZN6uqqlx1dbUnfzsajUfbmH7XC1x7Xhl3XTfJ6zgiIgCY2QrnXFUsX1O/LD6FvKwg7z93GH94YweHmtu8jiMiEjcqgtO4ecZIjrV18OSKOq+jiIjEjYrgNCaV5XHuiHx+/dpWXbRGRAYsFUEPbp5xFpv3HuHVLfu8jiIiEhcqgh68b/JQ8rOCumiNiAxYKoIeZAT9fKRqBM+uq6f+ULPXcUREYk5FEIUbp5fTEXI8rPMPicgApCKIwlmFg7hkbDEPL9tGW0fI6zgiIjGlIojSzTPOov5QCy+8We91FBGRmFIRRGn22SUMz8/U+YdEZMBREUTJ7zNumF7O0s37qN1z2Os4IiIxoyLohY9eMIKg3/jlUm0ViMjAoSLohaLsdK47r4xHl29n58FjXscREYkJFUEvffayChyOn/x5k9dRRERiQkXQS2WDs5g/rZzHqut4p+GI13FERM6YiqAPFs6uIOAz7nlRWwUikvxUBH1QkpvBx2aO5PerdrCpXt8gEpHkpiLoowWXjGFQWoAfvbDR6ygiImdERdBHBYPS+MSsUSxeu5uaHbqusYgkLxXBGfjUxaPIywzyw+e1VSAiyUtFcAZyM4J8+pLR/PmtPazYesDrOCIifaI
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 432x288 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {
|
||
|
|
"needs_background": "light"
|
||
|
|
},
|
||
|
|
"output_type": "display_data"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"import numpy as np #import appropriate libraries\n",
|
||
|
|
"import matplotlib.pyplot as plt\n",
|
||
|
|
"\n",
|
||
|
|
"start = 0 #start point of the domain over which the diffeq will be solved\n",
|
||
|
|
"end = 1 #end-point\n",
|
||
|
|
"steps = 30 # number of steps\n",
|
||
|
|
"t = np.linspace(start, end, steps) # t is a uniformly spaced vector spanning the start and endpoints\n",
|
||
|
|
"y = np.zeros(len(t)) # y is a vector initialized to 0 as long as t\n",
|
||
|
|
"h = (end - start)/steps # h is the width of each step\n",
|
||
|
|
"t = t*h # t is h times the number of t's gone before\n",
|
||
|
|
"\n",
|
||
|
|
"y[0] = 1 # initial value\n",
|
||
|
|
"\n",
|
||
|
|
"for i in range (0,len(t)-1): # initiate loop\n",
|
||
|
|
" y[i+1] = y[i] + h * (-6*y[i]) # euler method\n",
|
||
|
|
"plt.plot(t,y) # plot the graph"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": []
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"metadata": {
|
||
|
|
"interpreter": {
|
||
|
|
"hash": "a7b5a1112396ad393f5588006c146c6954be07cba002f50e558ad6b83bb9990c"
|
||
|
|
},
|
||
|
|
"kernelspec": {
|
||
|
|
"display_name": "Python 3.9.7 64-bit ('base': conda)",
|
||
|
|
"language": "python",
|
||
|
|
"name": "python3"
|
||
|
|
},
|
||
|
|
"language_info": {
|
||
|
|
"codemirror_mode": {
|
||
|
|
"name": "ipython",
|
||
|
|
"version": 3
|
||
|
|
},
|
||
|
|
"file_extension": ".py",
|
||
|
|
"mimetype": "text/x-python",
|
||
|
|
"name": "python",
|
||
|
|
"nbconvert_exporter": "python",
|
||
|
|
"pygments_lexer": "ipython3",
|
||
|
|
"version": "3.9.7"
|
||
|
|
},
|
||
|
|
"orig_nbformat": 4
|
||
|
|
},
|
||
|
|
"nbformat": 4,
|
||
|
|
"nbformat_minor": 2
|
||
|
|
}
|