9 years, 40 days

## variable...

gamma[o] is just a variable. actually i want a the variation of function with respect to gamma[o]

capital GAMMA is the GAMMA function in maple

and summation is over p.

## I am having maple 9.5 I recently opted m...

I am having maple 9.5

I recently opted maple for my simulations as it can solve special functions easily. I am totaly unaware of its functionalities and just trying to explore. So, my posts may bother you all.

Thanks for all suggestions.

## @Carl Love  Thanks for your solutio...

Thanks for your solution. It worked. But I want my vertical axis in log axis.

 > restart:
 > k:=2.01: m:=2: L1:=1: L2:=3: L3:=5: with(plots):

Warning, the name changecoords has been redefined

 > a1:=(m*L1)+k-1: a2:=(m*L2)+k-1: a3:=(m*L3)+k-1:
 > b1:=k-(m*L1): b2:=k-(m*L2): b3:=k-(m*L3):
 > p1:=((((k*m)/snr)^((a1+1)/2))*MeijerG([[-a1/2,(1-a1)/2],[]],[[b1/2,-b1/2],[-(1+a1)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L1)*GAMMA(k)):
 > p2:=((((k*m)/snr)^((a2+1)/2))*MeijerG([[-a2/2,(1-a2)/2],[]],[[b2/2,-b2/2],[-(1+a2)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L2)*GAMMA(k)):
 > p3:=((((k*m)/snr)^((a3+1)/2))*MeijerG([[-a3/2,(1-a3)/2],[]],[[b3/2,-b3/2],[-(1+a3)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L3)*GAMMA(k)):
 > result1:=[seq([i,evalf(eval(p1,snr=i))],i=1..30)];

 > result2:=[seq([i,evalf(eval(p2,snr=i))],i=1..30)]:
 > result3:=[seq([i,evalf(eval(p3,snr=i))],i=1..30)]:
 > plots[pointplot](result1);

 >
 >

## @Markiyan Hirnyk  >&n...

 > restart:
 > k:=2.01: m:=2: L1:=1: L2:=3: L3:=5: with(plots):

Warning, the name changecoords has been redefined

 > a1:=(m*L1)+k-1: a2:=(m*L2)+k-1: a3:=(m*L3)+k-1:
 > b1:=k-(m*L1): b2:=k-(m*L2): b3:=k-(m*L3):
 > p1:=((((k*m)/snr)^((a1+1)/2))*MeijerG([[-a1/2,(1-a1)/2],[]],[[b1/2,-b1/2],[-(1+a1)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L1)*GAMMA(k)):
 > p2:=((((k*m)/snr)^((a2+1)/2))*MeijerG([[-a2/2,(1-a2)/2],[]],[[b2/2,-b2/2],[-(1+a2)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L2)*GAMMA(k)):
 > p3:=((((k*m)/snr)^((a3+1)/2))*MeijerG([[-a3/2,(1-a3)/2],[]],[[b3/2,-b3/2],[-(1+a3)/2]],((k*m)/snr)))/(2*sqrt(Pi)*GAMMA(m*L3)*GAMMA(k)):
 > result1:=[seq([i,evalf(eval(p1,snr=i))],i=1..30)];

 > result2:=[seq([i,evalf(eval(p2,snr=i))],i=1..30)]:
 > result3:=[seq([i,evalf(eval(p3,snr=i))],i=1..30)]:
 > plots[pointplot](result1,axis[2]=[mode=log]);

Error, (in plot/options2d) unknown or bad argument: axis[2] = [mode = log]

 >
 >

## @Carl Love  > plots[pointplot(re...

> plots[pointplot(result1)];
plots[INTERFACE_PLOT(POINTS([1., 0.1459286401], [2., 0.08869547820],

[3., 0.06236755560], [4., 0.04727672818], [5., 0.03756222199],

[6., 0.03083117361], [7., 0.02592152085], [8., 0.02220128203],

[9., 0.01929787603], [10., 0.01697787638], [11., 0.01508789788],

[12., 0.01352319762], [13., 0.01220994115], [14., 0.01109467362],

[15., 0.01013780300], [16., 0.009309422223], [17., 0.008586548553],

[18., 0.007951250986], [19., 0.007389350534], [20., 0.006889499432],

[21., 0.006442517024], [22., 0.006040902793], [23., 0.005678474074],

[24., 0.005350092865], [25., 0.005051457327], [26., 0.004778940890],

[27., 0.004529466759], [28., 0.004300409197], [29., 0.004089515096],

[30., 0.003894841233]))]

using plots[pointplot] returned me the above values.

## @Carl Love  How can i use logplot c...

How can i use logplot command to simultaneously plot some equations.

That means  I am having three eqs p1, p2 and p3 and I want a logplot of these eqs with respect to i from 0 to 30.

## thanks...

@Axel Vogt  I made certain changes to my function and got the same plot as u are having. Thanks for helping. Regards.

## >  k:=1: &g...

 > k:=1:
 > m:=1:
 > a:=k-m:
 > b:=k+m-1:
 > z:=(k*m)/snr:
 > x:=(z^((b+1)/2))/(0.693*GAMMA(k)*GAMMA(m));

 > c:=x*MeijerG([[-(b+1)/2],[(1-b)/2]],[[a/2,-a/2,-(b+1)/2,-(b+1)/2],[]],z);

 > with(plots):
 > plot(c,snr=0..10);

 > kernelopts(version);

 >

 > restart: k:=5:
 > m:=10:
 > a:=k-m:
 > b:=k+m-1:
 > z:=(k*m)/snr:
 > x:=(z^((b+1)/2))/(0.693*GAMMA(k)*GAMMA(m));

 > c:=x*MeijerG([[-(b+1)/2],[(1-b)/2]],[[a/2,-a/2,-(b+1)/2,-(b+1)/2],[]],z);

Error, (in MeijerG/normal/parameters) function does not exist: the parameters lists are inconsistent in MeijerG([[-15/2], [-13/2]],[[-15/2, -15/2, -5/2, 5/2], []],<...>)

 > with(plots):

Warning, the name changecoords has been redefined

 > plot(c,snr=0..10);

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

Error, empty plot

 > kernelopts(version);

 >

## worksheet attached...

code is working for low values of m but if I choose k=5 and m=10 then error is prompted.

## same plot...

plotting with c and exact c gave me the same plot.  however for certain values of k amd m I am getting error for MeijerG function saying "parameters are inconsistent".

## @Preben Alsholm  OK. Thanks for sho...

OK. Thanks for showing me the plot. Will you please check out the same plot with snr defined between the range -5 to 10 and show me the output.

## @Preben Alsholm  I am having maple ...

I am having maple 9.5

## @Axel Vogt  I used "with(RealDomain...

I used "with(RealDomain)" to solve my csc function.

However avoiding the use of with(RealDomain) also doesn't overcome the problem.

Thanks.

## @gkokovidis  I changed Digits:=15: ...

I changed Digits:=15:

but changing this doesn't proved to be a solution for my problem.

## @Carl Love  Thanks for helping. I a...

Thanks for helping.

I am using maple 9.5. I started working on a new worksheet and it is woking fine. However I encountered new error while plotting my function i.e.  "unable to evaluate the function to numeric values in the region";

which I think is due to negative arguments in MeijerG as evaluation of my function gives

0.62*((1/s)^4.0)*MeijerG([[-3.5,3.0],[]],[[-2.0,2.0],[-4.0]],4.0/s)

and I want to plot this equation for various values of s.

Do I need to specify my special functions using "with" as we do in case of with(plot)  orany other way to remove this error?

 1 2 3 Page 2 of 3
﻿