Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 2,046 Was thanked: 1176 time(s) in 752 post(s)
|
Originally Posted by: omorr I can not figured out what is going on here by looking at the log file. Please look at the picture and the attached file. If I use numbers inside distrib package functions I have the result, but if I use variables - problems arise. Am I missing something here? This is a distrib package issue. By whatever reason, some functions in distrib seem to insist on floating point values for some of their arguments in order to issue a numeric result. Otherwise they return the unevaluated function form, which you do not see due to numeric optimization of your result display regions. Maxima(quantile_student_t(α/2,30))=quantile.student.t(1/40,30)Maxima(quantile_student_t(α/2,30.0))=-2.0423In the current experimental SVN-version of Maxima plugin, there is no conversion of underscores to periods any more. Here is the code of quantile_student_t:
/* R: qt(q,n) */ quantile_student_t(q,n):=block([cp:control1(n), cq:controlq(q), fq:float(q), aux, sgn], if cp=-1 or cq=-1 then error("Illegal parameter" , if cp=0 or cq=0 then return(funmake('quantile_student_t,[q,n])), if fq=0.0 then return('minf), if fq=1.0 then return('inf), if fq=0.5 then return(0), if numberp(q) and numberp(n) and (numer or ?floatp(q) or ?floatp(n)) then(if fq<0.5 then (aux : 2*fq, sgn : -1) else (aux : 2*(1-fq), sgn : 1), sgn*sqrt(n*(1 / ?iibeta(aux,float(n/2),0.5)-1)) ) else funmake('quantile_student_t,[q,n]) )$
|
|