TransWikia.com

Can the error function be expressed in terms of other special functions?

Mathematica Asked on January 30, 2021

I obtained with Mathematica some results written in terms of the error function Erfi[x].

Is there is a way to transform the error function into other special functions e.g. Bessel functions or others?

I wonder if it can be done with Mathematica.
Any suggestions are welcome.

U= Erfi[((1/2 + I/2) (R – z))/Sqrt[k R]] + Erfi[((1/2 + I/2) (R – Sqrt[D^2 + z^2]))/Sqrt[k R]]

2 Answers

There are many identities reminding one in the question (unclear what kind of relation has been intended), e.g. expressing FresnelS in terms of Erfi

FullSimplify[-1/4 (1 + I)( Erfi[(1 + I)Sqrt[Pi]z/2]- I Erfi[(1 - I)Sqrt[Pi]z/2])]
FresnelS[z]

Let's demonstrate other relations:

FullSimplify[{-1/4 (1+I)(I  Erfi[(1+I)Sqrt[Pi]z/2]- Erfi[(1-I)Sqrt[Pi]z/2]),
              -(Sqrt[-z^2]/z) - (z/Sqrt[Pi]) ExpIntegralE[1/2, -z^2],
              (z/Sqrt[-z^2])(1 - (E^z^2/Sqrt[Pi])HypergeometricU[1/2, 1/2, -z^2])
              }]//Column
  FresnelC[z]
  Erfi[z]
  Erfi[z]}

Such identities can be found exploiting MathematicalFunctionData and MathematicalFunction (the latter new in version 12), nontheless one can start with

MeijerGReduce[Erfi[x],x]

enter image description here

 Entity["MathematicalFunction","Erfi"]["AlternativeRepresentations"]

enter image description here enter image description here

as well as

Entity["MathematicalFunction","Erfi"]["HypergeometricRepresentations"]

enter image description here

enter image description here

A convenient way of exploring mathematical data involves Manipulate, e.g.

Manipulate[ Entity["MathematicalFunction","Erfi"][z], 
            {z, Entity["MathematicalFunction","Erfi"]["Properties"]}]

various items of the following can yield other identities:

MathematicalFunctionData["Properties"]

enter image description here

Answered by Artes on January 30, 2021

U = Erfi[((1/2 + I/2) (R - z))/Sqrt[k R]] + 
   Erfi[((1/2 + I/2) (R - Sqrt[D^2 + z^2]))/Sqrt[k R]];

Use ComplexityFunction to penalize the use of Erfi

U2 = FullSimplify[U, 
  ComplexityFunction -> (LeafCount[#] + 
      1000 Count[#, _Erfi, {0, Infinity}] &)]

(* (1 + I) (FresnelC[(R - z)/(Sqrt[π] Sqrt[k R])] + 
   FresnelC[(R - Sqrt[D^2 + z^2])/(Sqrt[π] Sqrt[k R])] + 
   I (FresnelS[(R - z)/(Sqrt[π] Sqrt[k R])] + 
      FresnelS[(R - Sqrt[D^2 + z^2])/(Sqrt[π] Sqrt[k R])])) *)

The expressions are equivalent

U == U2 // FullSimplify

(* True *)

However, you now have four special functions rather than two.

Answered by Bob Hanlon on January 30, 2021

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP