10

There is a typical typo here when in equation $r'=r^2\sqrt {…}$ only one branch is used, while there are two branches $r'=\pm r^2\sqrt{…}$. To take both branches into account, we must square the equation and differentiate by $\theta $ , then we obtain G = 6.672*(10^-11)(*Gravitational constant \ N·m\.b2/kg\.b2*); M = 5.965*10^24(*The mass property of ...


7

DateListPlot[Accumulate[RandomReal[{1, -0.5}, {93}]], {2011, 1}, FrameTicks -> {Automatic, {date[[;; ;; 6]], None}}, DateTicksFormat -> {"Month", "/", "Year"}, ImageSize -> 700] Thanks to Mr. Wizard: How to rotate TickMarks in DateListPlot? DateListPlot[Accumulate[RandomReal[{1, -0.5}, {93}]], {2011, 1}, FrameTicks -> {Automatic, {date[[;;...


5

I guess, the simplest way is following: SetDirectory["D:\\field data\\"]; fn=FileNames["f0*.txt"]; data=Table[Import[f,"Table"][[4;;10]],{f,fn}]; You can use the data as you wish..


5

Out-of-box solution (left figure): ListContourPlot[m[[All, 1 ;; 3]], ContourStyle -> Directive[Black, Dashed], ColorFunction -> "Rainbow", ImageSize -> Small, FrameLabel -> {"α", "δ"}, PlotTheme -> "Monochrome"] This is pure esthetics consideration that the contours are dashed, the image is small, the theme is monochrome. It is ...


4

How about: y1 = s[3, 1] AbsArgPlot[y1[x], {x, -5, 5}]


4

lpp3d1 = ListPointPlot3D[{#[[3 ;; 5]]} & /@ Sidata1, PlotStyle -> colorsSi, AxesLabel -> {"q", "s2", "α"}, ImageSize -> Large, LabelStyle -> {18}, PlotLabel -> Style["Si", 24], ViewPoint -> {-3, -2, 1}]; 1. Remove braces from {#[[3 ;; 5]]} & in the first argument of ListPointPlot3D and use replacement rule Point -> ...


3

The fix is something like this: PlotStyle -> (Directive[PointSize[0.01],#]&/@colorsSi) A few things to note: Only one PlotStyle can be used. If you use a few the 1st one will override the rest. Per each point specify all styles inside Directive Use (...) for proper operation order for complex PlotStyle and other options


3

See the Wolfram Function Repository function PhaseUnwrap: https://resources.wolframcloud.com/FunctionRepository/resources/PhaseUnwrap Note that this only works with numerical data, not with functions. Without any specifics on your data I can't make an example for you.


3

You are almost there. Try the option FrameTicksas follows: FrameTicks -> {{Automatic, None}, {Table[{10^i, Style[Superscript[10, i], 12, Italic, Blue] }, {i, 0, 10}], None}} Have fun! Edit: To address your question of how to remove the decimal points of ticks along the y axis. Instead of the Automatic you might insert this line: Table[...


3

Your ODE is stiff. There is a limit of 10,000 steps build in NDSolve and it reached that at Pi. You can see that from the output Notice the upper limit is 3.16 and not 2 Pi which is 6.28319 You can improve this by using Method -> "StiffnessSwitching" sol = NDSolve[{r'[θ] == r[θ]^2*Sqrt[(ε/p)^2 - (1/r[θ] - 1/p)^2], r[0] == r0}, r, {θ, 0, 2 Pi}, ...


3

I would do it this way. Show[ Plot[{1/x}, {x, 0, 26}, PlotRange -> {Automatic, {0, .6}}, PlotStyle -> AbsoluteThickness[3]], Plot[{1/x}, {x, 2, 6}, PlotRange -> {Automatic, {0, .6}}, PlotStyle -> Transparent, Filling -> Axis, FillingStyle -> Red], Plot[1/x, {x, 8, 24}, PlotRange -> {Automatic, {0, .6}},...


2

I want to use just single "Import" commond, because I have 30 data data file I do not see the need to use one import command. Why not just use a loop? See if this does what you want. (Not tested) fileName0 = "D:\\field data\\f0 ("; numOfFiles = 30; data = Table[0, {numOfFiles}]; Do[ fileName = fileName0 <> ToString[n] <> ").txt"; tmp =...


2

This problem can be vastly simplified. It involves BoxRatios, the fact that your data is essentially flat, and you are forcing it to a non-flat form factor which will scale the Tube to be very non-circular. Consider the following: With[{data = {{0, 0, 0}, {1, 1, 0.001}}}, Grid[{{ Graphics3D[Arrow[Tube[data]], BoxRatios -> {1, 1, 1}], Graphics3D[...


2

It misses it due to sampling. You can see this by asking to explicitly sample Pi/2 Plot[Floor[Sin[x]], {x, -Pi, Pi}, PlotPoints -> {50, {Pi/2}}, PlotLabel -> Floor[Sin[x]], Frame -> True, GridLines -> Automatic, GridLinesStyle -> Directive[LightGray, Dashed], FrameTicks -> {Automatic, {Table[x, {x, {-Pi, 0, Pi}}], None}}, ...


2

How it exactly should rotate is not clear to me (still image), but here is my interpretation of doing it: ClearAll[func,regpolygon] func[θ_]:={Cos[θ],Sin[θ]} regpolygon[n_][θ_]:=func[θ] Cos[Pi/n]/Cos[(2Pi)/n ((θ n)/(2Pi)-Floor[(θ n)/(2Pi)])-Pi/n] MakeScene[θ_,f_]:=Module[{plog,t,p1,pt0,pt1,pt2}, plog=ParametricPlot[f[t]-{1,0},{t,0,2Pi},PlotStyle->Red]; ...


2

star = Graphics@ WindingPolygon[CirclePoints[5][[{1, 3, 5, 2, 4}]]]; ContourPlot[{x^2 + y^2}, {x, -4, 4}, {y, -4, 4}, Contours -> {1, 9, 16}, ContourShading -> {Red, Green, Blue, None}, Epilog -> Inset[star, {0, 0}, Center, 1], ContourStyle -> {Red, Green, Blue}, FrameLabel -> {"a", "b"}] Use Epilog -> Inset[star,{0,0}, Center,...


2

The simplest way to find stability function is to use formula $$R(z)=\frac{Det[I-z\times A+z\times e\times b^T]}{Det[I-z\times A]}$$ where $A$ is Runge-Kutta matrix and $b$ is weights (Butcher tableau).R = Det[MI - z*A + z*ME.Transpose[B]]/Det[MI - z*A] In Wolfram Mathematica we could do that by using the function RungeKuttaLinearStabilityFunction For ...


1

Some progress can be made by integrating over y only (with the correction that {…} be replaced by (…)). Integrate[x DiracDelta[r x - y] Exp[1/g^2 (Cos[x2 - x] + Cos[x2] + Cos[y + x2])], {y, 0, 2 Pi}, Assumptions -> r > 0 && 0 < x < 2 Pi] (* E^((Cos[x - x2] + Cos[x2] + Cos[r x + x2])/g^2) x HeavisideTheta[2 Pi - r x] *) In contrast, ...


1

We can use ParametricNDSolveValue[], then we have f1[t_] := s*(y[t] - x[t]); f2[t_] := x[t]*(k - z[t]) - y[t]; f3[t_] := x[t]*y[t] - m*z[t]; system = {x'[t] == f1[t], y'[t] == f2[t], z'[t] == f3[t]}; ini = {x[0] == 1, y[0] == -1, z[0] == 10}; n = ParametricNDSolveValue[ Join[system, ini], {x[t], y[t], z[t]}, {t, 0, 100}, {s, m, k}]; Manipulate[ ...


1

Put a test of the argument to Sqrt to enforce it being positive before the inequality is checked: RegionPlot3D[(t2 + t4)^2 + 4 (t2 t4 - (1/2 (2*t1 + t2 + t4 - 1))^2) > 0 && 1/(1 + 2 (-(t2 + t4) + Sqrt[(t2 + t4)^2 + 4 (t2 t4 - (1/2 (2*t1 + t2 + t4 - 1))^2)])) >= 0, {t1, 0, 1}, {t2, 0, 1}, {t4, 0, 1}]


1

If I understand your question, I think the code below does the equivalent. You basically want to alter the x-values to account for the phase shift in rotFrameSin relative to labFrameSin, such that you can plot both on different x-values but the "shape" is the same? I've avoided ListAnimate as I am not all that familiar with that function. You can alter ...


1

ClearAll[st, sign, f, g] f[x_] := x^3 - x + 2; g[x_] := x^3 - 3*x + 1; sign[x_] := Sign[x] /. {-1 -> Negative, _ :> Positive} st = {StringTemplate["f(``) = <* f[#]*>."], StringTemplate["f'(``) is <* sign[f'[#]]*>."]}; pt1 = Plot[g[x], {x, -1, 3}, PlotRange -> {-2, 6}, ImageSize -> 200]; pt2 = Manipulate[st[[p]][x], {x, 0, 1}, ...


Only top voted, non community-wiki answers of a minimum length are eligible