Tag Info

New answers tagged

0

The issue at least on my machine (Linux, M10.1) is that the plots take a long time to generate and so creating a smooth animation with on-the-fly generated plots is, well, impossible. You could generate the plots beforehand, though. Here's some code to show you the progress of plot generation as well. Assuming that you have defined a function f[t] that ...


1

I guess there is some floating-point-related issue here... This works: f[t_] := RegionPlot[ TransformedRegion[ Rectangle[{-1, -1}, {1, 1}], { Indexed[#1, {1}] (1 + t (Indexed[#1, {2}]^2 - 1)) + 2 t, Indexed[#1, {2}] (1 + t (Indexed[#1, {1}]^2 - 1)) + 2 t } & ], PlotRange -> {{-1 + 2 t, 1 + 2 t}, {-1 + 2 t, 1 + 2 t}} ] ...


6

An FEM element-meshing approach. The quality is controlled by the option "MaxCellMeasure" -> {"Length" -> 0.05}. Note that the VertexNormals -> -coords option causes the polygonal sphere to be smoothed out when displayed on the screen. Needs["NDSolve`FEM`"]; points = {{-0.9207, -0.3896, 0.0091}, {-0.8272, 0.5077, -0.2399}, {0.2544, -0.3511, ...


11

Let me add another answer. This code is much shorter and faster than my previous one, and the resulting mesh of each face is much cleaner. The procedure is simple. Triangles are first made from the given face vertices and discretised. Each mesh point is then pushed onto a 2-sphere while its angular positions are maintained. points = { {-0.9207, -0.3896, ...


1

bug fixed in 10.1 (windows) code ArcLength[Line[{{0, 0}, {1, 0}, {2, 0}}]] ArcLength[Line[{{0}, {1}, {2}}]] ArcLength[Line[{{0, 0}, {1, 0}, {2.0, 0}}]] ArcLength[Line[{{0}, {1}, {2.0}}]]


11

What about some 2D Geo functionality for this? points = {{-0.9207, -0.3896, 0.0091}, {-0.8272, 0.5077, -0.2399}, {0.2544, -0.3511, 0.9010}, {0.3510, 0.6527, 0.6712}, {0.5436, -0.6326, -0.5513}, {0.6016, 0.2317, -0.7643}}; edges = {{1, 2}, {1, 3}, {1, 5}, {2, 4}, {2, 6}, {3, 4}, {3, 5}, {4, 6}, {5, 6}}; Construct the geodesics as GeoPath objects: latlons ...


13

Using the same initialization code as Taiki: origin = {0, 0, 0}; points = {{-0.9207, -0.3896, 0.0091}, {-0.8272, 0.5077, -0.2399}, {0.2544, -0.3511, 0.901}, {0.351, 0.6527, 0.6712}, {0.5436, -0.6326, -0.5513}, {0.6016, 0.2317, -0.7643}}; fs = {{1, 3, 5}, {1, 2, 4, 3}, {1, 2, 6, 5}, {3, 4, 6, 5}, {2, 4, 6}}; faces = points[[#]] & /@ fs; Then ...


15

A crude attempt This is for Mathematica 10+ only. To construct each face, I use an intersection between a unit 3-ball centred at the origin and a pyramid whose base is at infinity and apex is at the origin. Each edge of the pyramid passes through each vertex of the spherical face. The pyramid is given by ConicHullRegion[{origin}, {vertices}]. The ...


2

Try this: Clear["Global`*"]; m = 1; \[HBar] = 1; k = 1; V = -k/Sqrt[1 + x^2 + y^2]; A = 8; \[CapitalDelta] = 10^-3; SE[Etr_] := -\[HBar]^2/(2 m) \!\( \*SubsuperscriptBox[\(\[Del]\), \({x, y}\), \(2\)]\(\[Psi][x, y]\)\) + V \[Psi][x, y] - Etr \[Psi][x, y] == 0 \[CapitalOmega] = Disk[{0, 0}, A]; BC = DirichletCondition[\[Psi][x, y] == ...



Top 50 recent answers are included