[?] Plot the following function and check if it is right?

Questions
,
1

I need to create a log-log plot of the following function:

f1=((v^2/(k - 3/2) + 1)^-k Gamma[k + 1] v^2)/(k Sqrt[(k - 3/2) \[Pi]]
  Gamma[k - 0.5])

The plot given by mathematica is:

and the plot given by matlab is:

What could be reason for this? Please see the code of mathematica

ClearAll["Global`*"]

 distFunckappa = ((v^2/(k - 3/2) + 1)^-k Gamma[k + 1] v^2)/(
 k Sqrt[(k - 3/2) \[Pi]] Gamma[k - 0.5])



 k = 2;

LogLogPlot[distFunckappa, {v, 0, 10}, 
 PlotStyle -> Directive[Blue, Thick], 
 AxesStyle -> Directive[Black, Bold, FontSize -> 14], 
 TicksStyle -> Directive[Black, Bold], AxesLabel -> {"v", "DNF"}, 
 Axes -> True]

and code in matlab is:

k=2;

v=0:0.1:10;

kappaDistFunc=(gamma(k+1)./(k.*gamma(k-1/2).*sqrt(k-3/2).*sqrt(pi))).*(1+(v.^2./(k-3/2))).^(-k).*v.^2;

loglog(kappaDistFunc)

Please help me Thanks in advance

2

Mathematica gives a correct answer. Execute yours matlab code, give the same form as Mathematica.

3

My version of Mathematica is Mathematica 11.0. It gives me the first plot and matlab gives me the second plot and Matlab is correct. But people using Mathematica 11.3 is getting the same plot as matlab

4

My version of Mathematica is Mathematica 11.0. It gives me the first plot and matlab gives me the second plot and Matlab is correct. But people using Mathematica 11.3 is getting the same plot as matlab

5

Mathematica 10.2.

102.png
102.png660×425 8.44 KB

Mathematica 11.3.

113.png
113.png663×432 8.14 KB

6

See this: I have also done the same. What is the issue??

mathematicaa1.png
mathematicaa1.png1920×1080 141 KB

7

Probably a bug in LogLogPlot in Mathematica 11.0.0.

Try plot in Plot:

 Plot[(0.900316*v^2)/(1 + 2 v^2)^2, {v, 0, 10}]

And what Yours conclusions ?

8

Yes, the plot gives the correct answer. But I want to represent the data in log scale as loglog plot

9

Try:

distFunckappa = ((v^2/(k - 3/2) + 1)^-k Gamma[k + 1] v^2)/(k Sqrt[(k - 3/2) \[Pi]] Gamma[k - 0.5])
k = 2;

ListLogLogPlot[Table[{v, distFunckappa}, {v, 0.1, 10, 1/10}], 
PlotStyle -> Directive[Blue, Thick], 
AxesStyle -> Directive[Black, Bold, FontSize -> 14], 
TicksStyle -> Directive[Black, Bold], AxesLabel -> {"v", "DNF"}, Joined -> True]
10

Thank you, sir. This works correctly