Mathematica Asked on February 19, 2021
logLogRegionPlot[rplot_] :=
Module[{pts, pgon},
pts = Cases[Normal@rplot, Line[a__] :> a, Infinity];
pgon = {EdgeForm[],
Directive[RGBColor[0.368417, 0.506779, 0.709798],
AbsoluteThickness[1.6], Opacity[0.3]],
Cases[Normal@rplot, Polygon[_], Infinity]};
ListLogLogPlot[
pts,
Joined -> True, Frame -> True,
PlotRange -> All, AspectRatio -> 1,
Axes -> False, PlotStyle -> ColorData[1][1],
Epilog -> (pgon /. {x_, y_?NumericQ} :> Log@{x, y})
]
]
logLogRegionPlot@
RegionPlot[
{y > 8*(10^-10) (x)^(1/2)*HeavisideTheta[(x)^(-1) - (y)] &&
x > 6*(10^4) && x < 6*(10^10) &&
y < (8*(10^-10))^-1 x^(-5/2)*HeavisideTheta[-(x)^(-1) + (y)] &&
y > 0.6*x^(-3/2)},
{x, 10^2, 10^6}, {y, 10^-6, 10^-2},
PlotPoints -> 100
]
How can I produce a log region plot satisfying those inequalities?
I have tried to produce it with the above code. An exclusion region is coming, but I need a large plot range (xaxis€[10^2,10^14], yaxis€[10^-16,10^0]) for which it is giving a wrong plot.
This was already asked in Wolfram community
[Wolfram community asked question][1]
[1]: https://community.wolfram.com/groups/-/m/t/2116321?p_p_auth=mXPP5Nq0![My plot in MMA 12.1](https://i.stack.imgur.com/qwReG.png)
The change of the variables does the job.
RegionPlot[({y > 8*(10^-10) (x)^(1/2)*HeavisideTheta[(x)^(-1) - (y)] &&
x > 6*(10^4) && x < 6*(10^10) &&
y < (8*(10^-10))^-1 x^(-5/2)*HeavisideTheta[-(x)^(-1) + (y)] &&
y > 0.6*x^(-3/2)}) /. {x -> Exp[s], y -> Exp[t]}, {s, Log[10^2], Log[10^6]}, {t, Log[10^-6], Log[10^-2]}, PlotPoints -> 50]
Correct answer by user64494 on February 19, 2021
Get help from others!
Recent Questions
Recent Answers
© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP