summaryrefslogtreecommitdiff
path: root/graphics/asymptote/examples/impact.asy
diff options
context:
space:
mode:
authorNorbert Preining <norbert@preining.info>2019-09-02 13:46:59 +0900
committerNorbert Preining <norbert@preining.info>2019-09-02 13:46:59 +0900
commite0c6872cf40896c7be36b11dcc744620f10adf1d (patch)
tree60335e10d2f4354b0674ec22d7b53f0f8abee672 /graphics/asymptote/examples/impact.asy
Initial commit
Diffstat (limited to 'graphics/asymptote/examples/impact.asy')
-rw-r--r--graphics/asymptote/examples/impact.asy29
1 files changed, 29 insertions, 0 deletions
diff --git a/graphics/asymptote/examples/impact.asy b/graphics/asymptote/examples/impact.asy
new file mode 100644
index 0000000000..d903969777
--- /dev/null
+++ b/graphics/asymptote/examples/impact.asy
@@ -0,0 +1,29 @@
+// Contributed by Philippe Ivaldi.
+// http://www.piprime.fr/
+
+import graph3 ;
+import contour;
+size (6cm,0);
+currentprojection=orthographic(1,1,1) ;
+
+real rc=1, hc=2, c=rc/hc;
+draw(shift(hc*Z)*scale(rc,rc,-hc)*unitcone,blue);
+
+triple Os=(0.5,0.5,1);
+real r=0.5;
+draw(shift(Os)*scale3(r)*unitsphere,red);
+
+real a=1+1/c^2;
+real b=abs(Os)^2-r^2;
+
+real f(pair z)
+{
+ real x=z.x, y=z.y;
+ return a*x^2-2*Os.x*x+a*y^2-2*Os.y*y-2*Os.z*sqrt(x^2+y^2)/c+b;
+}
+
+real g(pair z){return (sqrt(z.x^2+z.y^2))/c;}
+
+draw(lift(g,contour(f,(-rc,-rc),(rc,rc),new real[]{0})),linewidth(2bp)+yellow);
+
+axes3("$x$","$y$","$z$",Arrow3);