How to update this code to the Statistic Package for the graph below?for n from 1 to 15 do num:=convert(evalf(n), string): tracker[n]:=textplot([18,0.3,`Lambda is `.num],color=blue): H[n]:=ProbHist(PoissonPDF(n,x),0..25,26): N[n]:=plot(NormalPDF(n,n,x),x = 0..25):  P[n]:=display({H[n],N[n],tracker[n]}): od: display([seq(P[n], n=1..15)], insequence=true,title="Normal Approx. to the Poisson. Lambda is increasing from 1 to 15");

One of the most basic decisions a baseball manager has to make is how to arrange the batting order.  There are many heuristics and models for estimating the productivity of a given order.  My personal favourite is the use of simulation, but by far the most elegant solution from a mathematical perspective uses probability matrices and Markov chains.  An excellent treatment of this topic can be found in Dr. Joel S. Sokol's article,

I suffer from floaters, (objects that appear in ones vision), a condition that is worsened by looking at white surfaces.  Is there a way of changing the background colour of a maple worksheet from white to a darker shade, perhaps a light grey or something?

Hi,

I have this equation:

gammainc(5,x) = pf

pf = [24*10^-4    120*10^-4    24*10^-3    120*10^-3    24*10^-2    120*10^-2    24*10^-1    120*10^-1    24*10^0]how to find the values of x for the corressponding pf?

Thank you

Jamal

Hello,

I like to plot the Bode diagram of function transfer. For that purpose, i use Dynamicsystem package.

I obtain my transfer function by 2 different methods : 
- direct method
- state space method.

The issue is the fact that i don't obtain the same results with the 2 methods.

 state_space_metho.pdf

I want construst a vector with 2^62 entries. Only 2^20+1=2049 entries are non-null complex numbers and others are all zero. Maple tells me the memory is not enough.

Recently I was able to try Maple 17 beside Maple 12.  There were little wiggles in 17's plot but not so in 12's.  I thought what if I dragged 12's plot into 17's?  What would I get?  Well here is the result.

I made a plot of x^2+10 in Maple 17 and a plot of x^2 in Maple 12 then dragged the line over. And this is the visual I got

Hi all,

This is just for fun.

I am originally from Shanghai China and I went abroad to UK in 2006. Maple has always been there for me. Maple was there for my senior year, for my foundation year, for all my undergraduate year, for my master's degree and for my PhD!

Somehow, I lost Maple 10. I will see if I can find it from my backup drive again!

A reddit user posted an interesting question:

http://www.reddit.com/r/math/comments/1eyx0f/help_with_simplifying_in_maple/

The problem is how to simplify, e.g. 2/3*37^(1/2)*sin(1/6*Pi+1/3*arccos(55/1369*37^(1/2)))+2/3 to 4.  I tried a few things and got nowhere.  Any suggestions?

#### Code attached at the BOTTOM ####

If I run these seperately, they seem to be running fine as shown above. But when I run '8' after '7', it does not seem to be working as in first picture.

I am trying to solve a problem using fsolve. Where the fsolve is inside the integral. I keep getting the following message. Here is what I am trying to do 

 

int(fsolve(int(k*r*PDF(R,r),r=-infinity..x)=0,x)*PDF(Q,z),z=-1..1);

where Q is a normal distribution and R is also normally distributed with its parameter conditional on Q. Following is the error message that I get.

"Error, (in fsolve) z is in the equation, and is not solved for"

Please, explain me, why I got this? And how to get normal one?

restart;

with(ExcelTools):

with(ListTools):

with(DynamicSystems):

filename := "MSFT";

close3 := Import(cat(cat("C://US//",filename),".xls"), filename, "E2:E100");

this i usually use ln(close3[n]/close[n-1]) as the original series

as i do not know whether Sample in statistics package can use this directly,

or 

need to use counting to classify them into some group and got a distribution graph

genDE := Diff(g(t),t)^2 - 1 - x*(g(t))^2 - (y-x^2)/12*(g(t))^4 = 0;gensol := dsolve(genDE,g(t));genfun := rhs(gensol[1]);

from above, i got generating function and then continue to derive got the following error

u = ln(-sqrt((-z^2+y)*(-3*z+sqrt(12*z^2-3*y)))*z^2/(sqrt(-z^3*(-3+2*sqrt(3)*csgn(z)))*(-z^2+y)))/x;

tau := subs(u=z,solve(q,z));

it said solution may have lost

if this can not be direct to...

restart;with(Statistics):

MartinPoisson := (p/((p+2*q*T)/(1+p*T+q*T^2)-(p*T+q*T^2)*(p+2*q*T)/(1+p*T+q*T^2)^2))^r*r^2*T^n/factorial(n);

Dist := subs(T=t,MartinPoisson) assuming t > 0;

MartinDist := Distribution(CDF = unapply(piecewise(t>1,Dist,0),t));

MartinDist := Distribution(CDF = unapply(Dist, t)) assuming t > 0; #Change 2

X:=RandomVariable(MartinDist);

MartinDensity := PDF(X,t);

MartinDensity := subs(t=x, MartinDensity);