Differences between revisions 47 and 54 (spanning 7 versions)

Sage Interactions - Calculus

goto interact main page

Contents

Sage Interactions - Calculus

Root Finding Using Bisection

by William Stein

Newton's Method

Note that there is a more complicated Newton's method below.

by William Stein

http://sagenb.org/home/pub/2824/

A contour map and 3d plot of two inverse distance functions

by William Stein

http://sagenb.org/home/pub/2823/

A simple tangent line grapher

by Marshall Hampton

Numerical integrals with the midpoint rule

by Marshall Hampton

Numerical integrals with various rules

by Nick Alexander (based on the work of Marshall Hampton)

Some polar parametric curves

by Marshall Hampton. This is not very general, but could be modified to show other families of polar curves.

Function tool

Enter symbolic functions f, g, and a, a range, then click the appropriate button to compute and plot some combination of f, g, and a along with f and g. This is inspired by the Matlab funtool GUI.

Newton-Raphson Root Finding

by Neal Holtz

This allows user to display the Newton-Raphson procedure one step at a time. It uses the heuristic that, if any of the values of the controls change, then the procedure should be re-started, else it should be continued.

Coordinate Transformations

by Jason Grout

Taylor Series

by Harald Schilly

Illustration of the precise definition of a limit

by John Perry

I'll break tradition and put the image first. Apologies if this is Not A Good Thing.

A graphical illustration of sin(x)/x -> 1 as x-> 0

by Wai Yan Pong

Quadric Surface Plotter

by Marshall Hampton. This is pretty simple, so I encourage people to spruce it up. In particular, it isn't set up to show all possible types of quadrics.

The midpoint rule for numerically integrating a function of two variables

by Marshall Hampton

Gaussian (Legendre) quadrature

by Jason Grout

The output shows the points evaluated using Gaussian quadrature (using a weight of 1, so using Legendre polynomials). The vertical bars are shaded to represent the relative weights of the points (darker = more weight). The error in the trapezoid, Simpson, and quadrature methods is both printed out and compared through a bar graph. The "Real" error is the error returned from scipy on the definite integral.

Vector Calculus, 2-D Motion

By Rob Beezer

A fast_float() version is available in a worksheet

Vector Calculus, 3-D Motion

by Rob Beezer

Available as a worksheet

Multivariate Limits by Definition

by John Travis

http://www.sagenb.org/home/pub/2828/

Directional Derivatives

This interact displays graphically a tangent line to a function, illustrating a directional derivative (the slope of the tangent line).

directional derivative.png

3D graph with points and curves

By Robert Marik

This sagelet is handy when showing local, constrained and absolute maxima and minima in two variables. Available as a worksheet

Approximating function in two variables by differential

by Robert Marik

Taylor approximations in two variables

by John Palmieri

This displays the nth order Taylor approximation, for n from 1 to 10, of the function sin(x² + y²) cos(y) exp(-(x²+y²)/2).

Volumes over non-rectangular domains

by John Travis

http://www.sagenb.org/home/pub/2829/

Lateral Surface Area

by John Travis

http://www.sagenb.org/home/pub/2826/

Parametric surface example

by Marshall Hampton

Line Integrals in 3D Vector Field

by John Travis

http://www.sagenb.org/home/pub/2827/

-  ⇤ ← Revision 47 as of 2012-03-16 06:53:53 → 
  Size: 57910
  Editor: jason
  Comment:
+   ← Revision 54 as of 2012-04-26 15:50:40 → ⇥
  Size: 57403
  Editor: jason
  Comment:
-Deletions are marked like this.
+Additions are marked like this.
 Line 6:
-[[HTML(<b>hi</b><script>alert('hi');</script>)]]
{{{#!rhtml
<script type="text/javascript" src="http://aleph.sagemath.org/static/jquery.min.js"></script>
<script type="text/javascript" src="http://aleph.sagemath.org/embedded_sagecell.js"></script>
<script type="text/javascript">alert('hi');</script>
    <script>
$(function() {
    var makecells = function() {
        sagecell.makeSagecell({
            inputLocation: '#interact1',
            evalButtonText: 'Interact'});
    }
    sagecell.init(makecells);
})</script>

}}}
-Line 26:
+Line 10:
-{{{#!html
<div id="interact1"><script type="text/code">
+{{{#!sagecell
-Line 66:
+Line 49:
-</script></div>
-Line 77:
+Line 59:
-{{{
+{{{#!sagecell
-Line 119:
+Line 101:
-{{{
+{{{#!sagecell
-Line 135:
+Line 117:
-{{{
+{{{#!sagecell
-Line 155:
+Line 137:
-{{{
+{{{#!sagecell
-Line 183:
+Line 165:
-{{{
+{{{#!sagecell
-Line 256:
+Line 238:
-{{{
+{{{#!sagecell
-Line 269:
+Line 251:
-{{{
+{{{#!sagecell
-Line 367:
+Line 349:
-{{{
+{{{#!sagecell
-Line 473:
+Line 455:
-{{{
+{{{#!sagecell
-Line 478:
+Line 460:
-def trans(x=input_box(u^2-v^2, label="x=",type=SR), \
         y=input_box(u*v+cos(u*v), label="y=",type=SR), \
         t_val=slider(0,10,0.2,6, label="Length of curves"), \
+def trans(x=input_box(u^2-v^2, label="x",type=SR),
         y=input_box(u*v+cos(u*v), label="y",type=SR),
-Line 483:
+Line 464:
-         u_range=input_box(range(-5,5,1), label="u lines"),
         v_range=input_box(range(-5,5,1), label="v lines")):
+         t_val=slider(0,10,0.2,6, label="Length"),
         u_range=input_box('[-5..5]', label="u lines"),
         v_range=input_box('[-5..5]', label="v lines")):
-Line 486:
+Line 468:
-Line 490:
+Line 473:
-     g1=sum([parametric_plot((i,v), (v,t_min,t_max), rgbcolor=(1,0,0)) for i in u_range])
     g2=sum([parametric_plot((u,i), (u,t_min,t_max), rgbcolor=(0,0,1)) for i in v_range])
     vline_straight=parametric_plot((u,v_val), (u,t_min,t_max), rgbcolor=(0,0,1), linestyle='-',thickness=thickness)
     uline_straight=parametric_plot((u_val, v), (v,t_min,t_max),rgbcolor=(1,0,0), linestyle='-',thickness=thickness)
+     uvplot=sum([parametric_plot((i,v), (v,t_min,t_max), color='red',axes_labels=['u','v'],figsize=[5,5]) for i in u_range])
     uvplot+=sum([parametric_plot((u,i), (u,t_min,t_max), color='blue',axes_labels=['u','v']) for i in v_range])
     uvplot+=parametric_plot((u,v_val), (u,t_min,t_max), rgbcolor=(0,0,1), linestyle='-',thickness=thickness)
     uvplot+=parametric_plot((u_val, v), (v,t_min,t_max),rgbcolor=(1,0,0), linestyle='-',thickness=thickness)
-Line 495:
+Line 478:
-     (g1+g2+vline_straight+uline_straight).save("uv_coord.png",aspect_ratio=1, figsize=[5,5], axes_labels=['$u$','$v$'])
-Line 500:
+Line 482:
-     g3=sum([parametric_plot((partial(xuv,i),partial(yuv,i)), (v,t_min,t_max), rgbcolor=(1,0,0)) for i in u_range])
     g4=sum([parametric_plot((partial(xvu,i),partial(yvu,i)), (u,t_min,t_max), rgbcolor=(0,0,1)) for i in v_range])
     uline=parametric_plot((partial(xuv,u_val),partial(yuv,u_val)),(v,t_min,t_max),rgbcolor=(1,0,0), linestyle='-',thickness=thickness)
     vline=parametric_plot((partial(xvu,v_val),partial(yvu,v_val)), (u,t_min,t_max), rgbcolor=(0,0,1), linestyle='-',thickness=thickness)
     (g3+g4+vline+uline).save("xy_coord.png", aspect_ratio=1, figsize=[5,5], axes_labels=['$x$','$y$'])
     print jsmath("x=%s, \: y=%s"%(latex(x), latex(y)))
     print "<html><table><tr><td><img src='cell://uv_coord.png'/></td><td><img src='cell://xy_coord.png'/></td></tr></table></html>"
+     xyplot=sum([parametric_plot((partial(xuv,i),partial(yuv,i)), (v,t_min,t_max), color='red', axes_labels=['x','y'],figsize=[5,5]) for i in u_range])
     xyplot+=sum([parametric_plot((partial(xvu,i),partial(yvu,i)), (u,t_min,t_max), color='blue') for i in v_range])
     xyplot+=parametric_plot((partial(xuv,u_val),partial(yuv,u_val)),(v,t_min,t_max),color='red', linestyle='-',thickness=thickness)
     xyplot+=parametric_plot((partial(xvu,v_val),partial(yvu,v_val)), (u,t_min,t_max), color='blue', linestyle='-',thickness=thickness)
     html("$$x=%s, \: y=%s$$"%(latex(x), latex(y)))
     html.table([[uvplot,xyplot]])
-Line 515:
+Line 496:
-{{{
+{{{#!sagecell
-Line 538:
+Line 519:
-{{{
+{{{#!sagecell
-Line 565:
+Line 546:
-{{{
+{{{#!sagecell
-Line 591:
+Line 572:
-{{{
+{{{#!sagecell
-Line 606:
+Line 587:
-{{{
+{{{#!sagecell
-Line 652:
+Line 633:
-{{{
from scipy.special.orthogonal import p_roots
+{{{#!sagecell
import scipy
import numpy
from scipy.special.orthogonal import p_roots, t_roots, u_roots
-Line 662:
+Line 645:
-                'Chebyshev': {'w': 1/sqrt(1-x**2), 'xmin': -1, 'xmax': 1, 'func': t_roots},
                'Chebyshev2': {'w': sqrt(1-x**2), 'xmin': -1, 'xmax': 1, 'func': u_roots},
                'Trapezoid': {'w': 1, 'xmin': -1, 'xmax': 1, 'func': lambda n: (linspace(-1r,1,n), numpy.array([1.0r]+[2.0r]*(n-2)+[1.0r])*1.0r/n)},
                'Simpson': {'w': 1, 'xmin': -1, 'xmax': 1, 'func': lambda n: (linspace(-1r,1,n), numpy.array([1.0r]+[4.0r,2.0r]*int((n-3.0r)/2.0r)+[4.0r,1.0r])*2.0r/(3.0r*n))}}
+     'Chebyshev': {'w': 1/sqrt(1-x**2), 'xmin': -1, 'xmax': 1, 'func': t_roots},
     'Chebyshev2': {'w': sqrt(1-x**2), 'xmin': -1, 'xmax': 1, 'func': u_roots},
     'Trapezoid': {'w': 1, 'xmin': -1, 'xmax': 1,          'func': lambda n: (linspace(-1r,1,n), numpy.array([1.0r]+[2.0r]*(n-2)+[1.0r])*1.0r/n)},
     'Simpson': {'w': 1, 'xmin': -1, 'xmax': 1,          'func': lambda n: (linspace(-1r,1,n), 
            numpy.array([1.0r]+[4.0r,2.0r]*int((n-3.0r)/2.0r)+[4.0r,1.0r])*2.0r/(3.0r*n))}}
-Line 669:
+Line 655:
-    return polygon([(center-width2,0),(center+width2,0),(center+width2,height),(center-width2,height)],**kwds)
+    return polygon([(center-width2,0),
        (center+width2,0),(center+width2,height),(center-width2,height)],**kwds)
-Line 673:
+Line 660:
-def weights(n=slider(1,30,1,default=10),f=input_box(default=3*x+cos(10*x)),show_method=["Legendre", "Chebyshev", "Chebyshev2", "Trapezoid","Simpson"]):
+def weights(n=slider(1,30,1,default=10),f=input_box(default=3*x+cos(10*x),type=SR),
    show_method=["Legendre", "Chebyshev", "Chebyshev2", "Trapezoid","Simpson"]):
-Line 682:
+Line 670:
-    scaled_ff = fast_float(scaled_func)
+    scaled_ff = fast_float(scaled_func, 'x')
-Line 690:
+Line 678:
-    stems = sum(line([(x,0),(x,scaled_ff(x))],rgbcolor=(1-y,1-y,1-y),thickness=2,markersize=6,alpha=y) for x,y in coords_scaled)
    points = sum([point([(x,0),(x,scaled_ff(x))],rgbcolor='black',pointsize=30) for x,_ in coords])
+    stems = sum(line([(x,0),(x,scaled_ff(x))],rgbcolor=(1-y,1-y,1-y),
        thickness=2,markersize=6,alpha=y) for x,y in coords_scaled)
    points = sum([point([(x,0),
        (x,scaled_ff(x))],rgbcolor='black',pointsize=30) for x,_ in coords])
-Line 696:
+Line 686:
-    show(graph,xmin=plot_min,xmax=plot_max)
+    show(graph,xmin=plot_min,xmax=plot_max,aspect_ratio="auto")
-Line 704:
+Line 694:
-    html("$$\sum_{i=1}^{i=%s}w_i\left(%s\\right)= %s\\approx %s =\int_{-1}^{1}%s \,dx$$"%(n,latex(f.subs(x="x_i")), approximation, integral, latex(scaled_func)))
+    html("$$\sum_{i=1}^{i=%s}w_i\left(%s\\right)= %s\\approx %s =\int_{-1}^{1}%s \,dx$$"%(n,
        latex(f), approximation, integral, latex(scaled_func)))
-Line 717:
+Line 708:
-{{{
+{{{#!sagecell
-Line 829:
+Line 820:
-{{{
+{{{#!sagecell
-Line 866:
+Line 857:
-velocity = derivative( position(t) )
acceleration = derivative(velocity(t))
+velocity = derivative( position(t), t)
acceleration = derivative(velocity(t), t)
-Line 869:
+Line 860:
-speed_deriv = derivative(speed)
+speed_deriv = derivative(speed, t)
-Line 871:
+Line 862:
-dT = derivative(tangent(t))
+dT = derivative(tangent(t), t)
-Line 874:
+Line 865:
-## dB = derivative(binormal(t))
+## dB = derivative(binormal(t), t)
-Line 957:
+Line 948:
-{{{
+{{{#!sagecell
-Line 992:
+Line 983:
-{{{
+{{{#!sagecell
-Line 1029:
+Line 1020:
-{{{
+{{{#!sagecell
-Line 1082:
+Line 1073:
-{{{
%hide
%auto
+{{{#!sagecell
-Line 1143:
+Line 1132:
-{{{
+{{{#!sagecell
-Line 1181:
+Line 1170:
-{{{
+{{{#!sagecell
-Line 1221:
+Line 1210:
-{{{
+{{{#!sagecell
-Line 1328:
+Line 1317:
-{{{
+{{{#!sagecell
-Line 1385:
+Line 1374:
-{{{
+{{{#!sagecell
-Line 1406:
+Line 1395:
-{{{
+{{{#!sagecell

Diff for "interact/calculus"

Sage Interactions - Calculus

Root Finding Using Bisection

Newton's Method

A contour map and 3d plot of two inverse distance functions

A simple tangent line grapher

Numerical integrals with the midpoint rule

Numerical integrals with various rules

Some polar parametric curves

Function tool

Newton-Raphson Root Finding

Coordinate Transformations

Taylor Series

Illustration of the precise definition of a limit

A graphical illustration of sin(x)/x -> 1 as x-> 0

Quadric Surface Plotter

The midpoint rule for numerically integrating a function of two variables

Gaussian (Legendre) quadrature

Vector Calculus, 2-D Motion

Vector Calculus, 3-D Motion

Multivariate Limits by Definition

Directional Derivatives

3D graph with points and curves

Approximating function in two variables by differential

Taylor approximations in two variables

Volumes over non-rectangular domains

Lateral Surface Area

Parametric surface example

Line Integrals in 3D Vector Field