Calc2-Lab6-Parametric-Curves.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "a33db59f",
   "metadata": {},
   "source": [
    "# Calculus 2, Lab 6: Parametric curves\n",
    "\n",
    "Recall that a **parametric curve** is given by equations $x=f(t), y=g(t)$, for $t$ in some domain (often a closed interval $[a,b]$).\n",
    "\n",
    "Using parametric equations allows us to define a wide variety of curves that would be impossible to define using the graph of a function, or an implicitly-defined function.\n",
    "\n",
    "The main things we are interested in with parametric curves are:\n",
    "\n",
    "1. How to plot them (this is mostly a job for the computer)\n",
    "2. How to convert (when possible) between \"rectangular\" curves like $y=x^2$ or $\\frac{x^2}{4}+\\frac{y^2}{9}=1$ and parametric curves. (Converting from rectangular to parametric is called \"parametrization\", and there is no unique answer for this. Converting from parametric to rectangular is called \"eliminating the parameter\".)\n",
    "3. Calculating tangent and normal lines to parametric curves, and determining if a curve is \"smooth\"\n",
    "4. Calculating arc length (this is mostly a job for the computer)\n",
    "5. Calculating area bounded by closed curves (this is the topic of Assignment 9)\n",
    "\n",
    "In this notebook you will learn how to plot parametric curves using SymPy and GeoGebra, and you will calculate arc length symbolically (using SymPy) and numerically (using NumPy).\n",
    "\n",
    "## Plotting Parametric curves with SymPy\n",
    "\n",
    "The basic syntax for plotting a parametric curve in SymPy is `plot_parametric((expr_x,expr_y),range)`, where `expr_x` is an expression that defines $x$ in terms of our parametric (most likely `t`), `expr_y` is an expression that defines $y$ in terms of the parameter, and `range` defines the domain for the parameter, using the same syntax as for integrals. For example, for $t\\in [1,3]$ we would use `(t,1,3)` for the range.\n",
    "\n",
    "You can plot more than one curve on the same set of axes by including an ordered list of curves: for two curves, you would have something like `plot_parametric((expr_x1,expr_y1),(expr_x2,expr_y2),range)` if the range is the same, or `plot_parametric((expr_x1,expr_y1,range1),(expr_x2,expr_y2,range2))` for different ranges.\n",
    "\n",
    "Optional arguments include `line_color` for the colour of the curve, `label`, `xlabel`, and `ylabel` to label the plot, $x$ axis, and $y$ axis, respectively, and `xlim`, `ylim` to set the bounds for the $x$ and $y$ axes (equal to an ordered pair).\n",
    "\n",
    "Here is an example, which plots the ellipse with parametric equations $x=4\\cos(t), y=3\\sin(t)$, for $t\\in [0,2\\pi]$."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "95add614",
   "metadata": {},
   "outputs": [],
   "source": [
    "import sympy as sy\n",
    "sy.init_printing()\n",
    "t = sy.symbols('t')\n",
    "sy.plot_parametric((4*sy.cos(t),3*sy.sin(t)),(t,0,2*sy.pi))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3741c91e",
   "metadata": {},
   "source": [
    "You can also plot more than one curve on the same axes: just include multiple ordered pairs. For two curves with the same parameter domain, the syntax looks like\n",
    "\n",
    "```\n",
    "sy.plot_parametric((f1(t),g1(t)),(f2(t),g2(t)),(t,a,b))\n",
    "```\n",
    "For example, the following will plot the curves $x=t,y=\\cos(t)$ (cosine curve) and $x=t, y=\\sin(t)$ (sine curve), for $t\\in [-2\\pi,2\\pi]$:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d74db8d2",
   "metadata": {},
   "outputs": [],
   "source": [
    "sy.plot_parametric((t,sy.cos(t)),(t,sy.sin(t)),(t,-2*sy.pi,2*sy.pi))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "66eb1f21",
   "metadata": {},
   "source": [
    "One catch: if you wanted to specify a different colour, you can't do this all in one go, because the colour gets applied to both curves."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "65d24b5c",
   "metadata": {},
   "outputs": [],
   "source": [
    "sy.plot_parametric((t,sy.cos(t)),(t,sy.sin(t)),(t,-2*sy.pi,2*sy.pi),line_color='red')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "92aa5d0d",
   "metadata": {},
   "source": [
    "But we can use the following to override the colour for one of the curves:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0f4d326e",
   "metadata": {},
   "outputs": [],
   "source": [
    "p=sy.plot_parametric((t,sy.cos(t)),(t,sy.sin(t)),(t,-2*sy.pi,2*sy.pi),line_color='red',show=False)\n",
    "p[0].line_color='blue'\n",
    "p.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f96470f9",
   "metadata": {},
   "source": [
    "### Problem 1\n",
    "\n",
    "Plot the following parametric curves. (You probably want to plot each on its own set of axes.)\n",
    "\n",
    "If you don't want to type out `sy.plot_parametric` every time, you could add a line like `pp=sy.plot_parametric`.\n",
    "\n",
    "(a) $x=t^3-2t, y=t^2-t, t\\in [-1,2]$\n",
    "\n",
    "(b) $x=t+\\sin(2t), y=t+\\sin(3t), t\\in [-2\\pi, 2\\pi]$\n",
    "\n",
    "(c) $x=\\sin(3t), y=\\sin(4t), t\\in [0,2\\pi]$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "98f78fd1",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "bc77594f",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "218d9d55",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "fc94e5cc",
   "metadata": {},
   "source": [
    "## Plotting parametric curves with GeoGebra\n",
    "\n",
    "Sometimes you want something a little simpler to use than a Jupyter notebook. Fortunately, there are some great free software programs out there, such as Desmos, and [GeoGebra](https://geogebra.org). GeoGebra can do derivatives, integrals, geometry, statistics, and more. And the interface involves a little more pointing and clicking, if you like that sort of thing."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "333ba6f0",
   "metadata": {},
   "source": [
    "We will use [GeoGebra classic](https://www.geogebra.org/classic). Choose graphing mode to create your plot.\n",
    "\n",
    "You can create plots using the input box, which is next to the $+$ sign at the top-left.\n",
    "\n",
    "The syntax for a parametric curve is: `Curve(Expression,Expression,Parameter Variable,Start value,End value)`. For example, you could plot the unit circle using `Curve(cos(t),sin(t),t,0,2*pi)`.\n",
    "\n",
    "### Problem 2\n",
    "\n",
    "Use GeoGebra to plot the curve $x=\\cos(t+\\cos(2t)),y=\\sin(t+\\sin(4t)), t\\in[0,2\\pi]$.\n",
    "\n",
    "When you are happy with your plot, click on the menu button in the top-right corner, and choose Export Image.\n",
    "\n",
    "Then, double-click the markdown cell below, and add your image. You should be able to just copy it to the clipboard, and then use `Ctrl-V` to paste the image into the cell.\n",
    "\n",
    "Or, to play it safe, download the image from GeoGebra, upload it to Jupyter in the same folder as this notebook, and include it using the syntax `![my amazing image](myimage.png)`, assuming that you saved the file as `myimage.png`."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c1eec21a",
   "metadata": {},
   "source": [
    "**Double-click this text** to add your image, and hit `Ctrl-Enter` when you're done."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "45d4dacc",
   "metadata": {},
   "source": [
    "\n",
    "### Problem 3\n",
    "\n",
    "Use parametric curves to generate an emoji of your choice.\n",
    "\n",
    "**Note**: you will need more than one curve. Unless you are feeling ambitious, I recommend sticking to the classic smiley face."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "cf059911",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "62df27d3",
   "metadata": {},
   "source": [
    "## Arc length for parametric curves\n",
    "\n",
    "The formula for arc length for parametric curves is a fairly straightforward extension of the one we learned earlier. For a graph $y=f(x)$, we started with $\\Delta s = \\sqrt{\\Delta x^2+\\Delta y^2}$ and factored out the $\\Delta x$ to get $\\Delta s = \\sqrt{1+\\left(\\frac{\\Delta y}{\\Delta x}\\right)^2}\\Delta x$. Taking the limit gave us $ds = \\sqrt{1+f'(x)^2}dx$.\n",
    "\n",
    "Now, given $x(t)$ and $y(t)$, we have $dx = x'(t)\\,dt$ and $dy = y'(t)\\,dt$. The same sort of argument gives us\n",
    "\n",
    "$$ds = \\sqrt{x'(t)^2+y'(t)^2}\\,dt.$$\n",
    "\n",
    "### Problem 4\n",
    "\n",
    "For the curve given by $x=\\dfrac{t}{1+t}, y=\\ln(1+t), t\\in [0,2]$:\n",
    "\n",
    "(a) Set up the integral using SymPy, and display the integral. (To display an integral without evaluating it, use `Integral` instead of `integrate`.) You will need to define the expressions for $x$ and $y$ and compute their derivatives before setting up the integral."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "103c353e",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "c320f578",
   "metadata": {},
   "source": [
    "(b) Evaluate the integral using SymPy. \n",
    "\n",
    "**Note**: if SymPy fails to evaluate the integral, try using the change of variables $u=t+1$, and see if it succeeds with the resulting integral. Don't forget to adjust the bounds!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "7d27eeec",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "4fbba08f",
   "metadata": {},
   "source": [
    "### Problem 5\n",
    "\n",
    "The equations\n",
    "$$\\begin{aligned}x&=11\\cos(t)-4\\cos(11t/2)\\\\y&=11\\sin(t)-4\\sin(11t/2)\\end{aligned}$$\n",
    "define a curve called the **epitrochoid**.\n",
    "\n",
    "(a) Plot the curve using SymPy, and experiment to determine what interval gives the complete length of the curve.\n",
    "\n",
    "(b) Approximate the length of the curve using Simpson's Rule with 50 subintervals."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "46a4aa6b",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.11.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}