## Plot The volume of a cylinder between a cone and XY

I want to plot the volume of a cylinder between a cone and the plane XY The equation of the cone is zz == √x2+y2\sqrt{x^2 + y^2} and the cylinder equation y2−2y+x2=0y^2 – 2y + x^2 = 0 I am using the following code to plot the graphs: p3 = …

## Interpreting P-Value on ParameterTable

When using LinearModelFit[] we can generate a table of parameter values and statistics for the corresponding model using the property “ParameterTable”. Looking in the documentation I can’t find an explicit statement of the hypothesis test being performed for the p-values listed in this table. I would like to make sure …

## PlotLegends in ListPlot does not display different markers

I would like to plot a graph with the legends at its side. However, I have managed to plot the legends, but the markers do not match the one in plot (picture bellow): My code for that plot is: ListPlot[{Table[{(n – 1/2)*1/1000, VelocityU[[1, n]]}, {n, 1, NN}], Table[{(n – 1/2)*1/1000, …

## Manipulating List?

I have a huge list with 1000 elements. Manually entering and calculating SUM of probabilities it would take a few hours. Example: list = {1, 2, 3, 3, 4, 3} Calculating probabilities: Counts[list]/Length[list] 1/6, 2 -> 1/6, 3 -> 1/2, 4 -> 1/6|> Calculating the sum of the list elements …

## FindMinimum function does not produce minimum from near point

Using my simple program and I would like to find such matrix elements that minimize difference of their eigenvalues. H1 = {{1, 2}, {3, 4}}; H2 = {{2, 2}, {3, 2}}; f1[R_, S_] := R*H1 + S*IdentityMatrix[2] f2[R_, S_] := S*H2 + R*IdentityMatrix[2] Evl1 = N[Eigenvalues[f1[R, S]]]; Evl2 = N[Eigenvalues[f2[R, …

## Removing solutions which are particular cases of others

I have an algorithm which gives me a list of solutions (in matrix form, if it is of any importance) of an equation. Some of these solutions are particular cases of other. For example, in this list there can be solutions (10a0b1−a000)\begin{pmatrix} 1 & 0 & a \\ 0 & …

## Animation of a plane wave

I want to animate in 3D the plane wave sin(ωt+kx+ϕ)sin(ωt+kx+ϕ)\sin(\omega t+kx+\phi), but following the form from the documentation, I can only make kkk or ωω\omega vary, and I want ttt to change. How do I do that? ================= 1   Welcome! Please do add more info, especially your code. Consider …

## Table of ParametricPlots crashes Kernel

I can usually crash the Mathematica Kernel by using Table to generate a List of ParametricPlots, as follows: testpara[α_] := ParametricPlot[ { {Cos[θ], Sin[θ]}, {2 Cos[α] + Cos[θ], 2 Sin[α] + Sin[θ]}, {r, 0} }, {θ, 0, 2 π}, {r, 1, 2}, PlotRange -> 3, Frame -> False ] Table[testpara[α], …

## What particular situation do we have to use RuleDelayed after replacement? [closed]

On the answer of How to impose normalization condition (not initial condition) on DSolve or NDSolve? k = vars /. sol /. cs[[1]] /. C[_] :> 0 // N; I see that in the documentation of ReplaceAll and examples of ReplaceAll on the internet, we use Rule (->) instead of …

## How to Fit experimental data with 4D ParametricNDSolve?

Let me first explain the background problem. I’m trying to solve a problem with a projectile motion in 4D such that the position vector is: r(t)=r[x(t),y(t),z(t)]\boldsymbol{r}(t)=\boldsymbol{r}[x(t),y(t),z(t)] I have the differential equations with F(t)=mr″(t)F(t)=mr”(t) such that: \boldsymbol{F}_{drag}(t)=\frac{1}{2} \rho A C_{drag} R (\boldsymbol{r’}(t))^2; \\ \boldsymbol{F}_{grav}=m\boldsymbol{g};\\ \boldsymbol{F}_{total}(t)=\boldsymbol{F}_{drag}(t)+\boldsymbol{F}_{grav}\boldsymbol{F}_{drag}(t)=\frac{1}{2} \rho A C_{drag} R (\boldsymbol{r’}(t))^2; \\ …