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
