ASKSAGE: Sage Q&A Forum - RSS feedhttps://ask.sagemath.org/questions/Q&A Forum for SageenCopyright Sage, 2010. Some rights reserved under creative commons license.Fri, 06 Sep 2013 14:20:20 +0200Fermi-Dirac integral of half orderhttps://ask.sagemath.org/question/10517/fermi-dirac-integral-of-half-order/I am trying to implement a 1D model for semiconductor pn-junctions. This involves evaluating the Fermi-Dirac integral of half order, which can only be done numerically. I saw that the GSL library has an algorithm to perform this integration, but I can't figure out how to call it in Sage, is this possible? If not, what are some alternatives to evaluate this integral? Thanks in advance.wlp2Fri, 06 Sep 2013 14:20:20 +0200https://ask.sagemath.org/question/10517/How can one use maxima kummer confluent functions in sagehttps://ask.sagemath.org/question/8727/how-can-one-use-maxima-kummer-confluent-functions-in-sage/Hi, here is my piece of code
var('x,m')
assume(m, 'integer')
y = function('y', x)
desolve(diff(y,x,2) + 2*x*diff(y,x) - 4*m*y, y,contrib_ode=true,ivar=x)
which yields
[y(x) == k1*kummer_m(-m, 1/2, -x^2) + k2*kummer_u(-m, 1/2, -x^2)]
Internet(Wikipedia) sais that the solutions to this differential equation are the Kummer functions or Confluent_hypergeometric_functions. I now want to know how I can use these functions in sage since they must be defined in maxima.
Thanks in advancegod.oneTue, 21 Feb 2012 05:20:27 +0100https://ask.sagemath.org/question/8727/