help links subscribe changes refresh edit

	SandBoxFunctionalSpecialFunction
	↑

SandBoxDiracDelta

last edited 9 years ago by Bill page

Edit detail for SandBoxDiracDelta revision 38 of 40

	1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40
Editor: Bill page Time: 2014/08/19 18:03:06 GMT+0
Note: problems

removed:
---integrate(log(x)*unitDoublet(x-t),x=%minusInfinity..%plusInfinity,"noPole")

added:
-- loops forever on %plusInfinity limit of indefinite integral
--integrate(log(x)*unitDoublet(x-t),x=%minusInfinity..%plusInfinity,"noPole")

Distributions

We represent distributions as a complex function in the limit where the imaginary part goes to 0. Signum (sign), delta and doublet functions arise from the derivative of the absolute value function with respect to its real part.

Ref.

Derivatives of Distributions (Wikipedia)

Use definition of conjugate and derivative of abs(x) from SandBoxFunctionalSpecialFunction?

fricas

)lib FSPECX

   FunctionalSpecialFunction is now explicitly exposed in frame initial

   FunctionalSpecialFunction will be automatically loaded when needed 
      from /var/aw/var/LatexWiki/FSPECX.NRLIB/FSPECX
wirtingerD(ex,z) == eval(D(eval(ex,z=%conjugate),%conjugate),%conjugate=z)

Type: Void

Absolute value as a function of real or complex variables

fricas

abs(z)==sqrt(z*conjugate(z)); abs(z)

fricas

Compiling function abs with type Variable(z) -> Expression(Integer)

$\label{eq1}\sqrt{z \ {\overline z}}$

(1)

Type: Expression(Integer)

fricas

abz(a)==sqrt((a+%i*b)*(a-%i*b)); abz(x)

fricas

Compiling function abz with type Variable(x) -> Expression(Complex(
      Integer))

$\label{eq2}\sqrt{{{x}^{2}}+{{b}^{2}}}$

(2)

Type: Expression(Complex(Integer))

Total Wirtinger derivatives of this complex function (CR-calculus)

fricas

abs':=wirtingerD(abs(z),z)+wirtingerD(abs(z),conjugate z)

fricas

Compiling function wirtingerD with type (Expression(Integer),
      Variable(z)) -> Expression(Integer)

fricas

Compiling function wirtingerD with type (Expression(Integer),
      Expression(Integer)) -> Expression(Integer)

$\label{eq3}{{\overline z}+ z}\over{2 \ {\sqrt{z \ {\overline z}}}}$

(3)

Type: Expression(Integer)

fricas

abs'':=wirtingerD(abs',z)+wirtingerD(abs',conjugate z)

$\label{eq4}{-{{\overline z}^{2}}+{2 \ z \ {\overline z}}-{{z}^{2}}}\over{4 \ z \ {\overline z}\ {\sqrt{z \ {\overline z}}}}$

(4)

Type: Expression(Integer)

fricas

abs''':=wirtingerD(abs'',z)+wirtingerD(abs'',conjugate z)

$\label{eq5}{{3 \ {{\overline z}^{3}}}-{3 \ z \ {{\overline z}^{2}}}-{3 \ {{z}^{2}}\ {\overline z}}+{3 \ {{z}^{3}}}}\over{8 \ {{z}^{2}}\ {{\overline z}^{2}}\ {\sqrt{z \ {\overline z}}}}$

(5)

Type: Expression(Integer)

Consider these as complex functions of for real and .

fricas

realSignum:=eval(eval(abs',z=a+%i*b),[conjugate(a)=a,conjugate(b)=b])

$\label{eq6}a \over{\sqrt{{{b}^{2}}+{{a}^{2}}}}$

(6)

Type: Expression(Complex(Integer))

fricas

test(realSignum=D(abz(a),a))

fricas

Compiling function abz with type Variable(a) -> Expression(Complex(
      Integer))

$\label{eq7} \mbox{\rm true}$

(7)

Type: Boolean

fricas

signum(x)==eval(realSignum,a=x); signum(x)

fricas

Compiling function signum with type Variable(x) -> Expression(
      Complex(Integer))

$\label{eq8}x \over{\sqrt{{{x}^{2}}+{{b}^{2}}}}$

(8)

Type: Expression(Complex(Integer))

fricas

--
realDirac:=eval(eval(abs'',z=a+%i*b),[conjugate(a)=a,conjugate(b)=b])

$\label{eq9}{{b}^{2}}\over{{\left({{b}^{2}}+{{a}^{2}}\right)}\ {\sqrt{{{b}^{2}}+{{a}^{2}}}}}$

(9)

Type: Expression(Complex(Integer))

fricas

test(realDirac=D(abz(a),[a,a]))

$\label{eq10} \mbox{\rm true}$

(10)

Type: Boolean

fricas

diracDelta(x)==eval(realDirac/2,a=x); diracDelta(x)

fricas

Compiling function diracDelta with type Variable(x) -> Expression(
      Complex(Integer))

$\label{eq11}{{b}^{2}}\over{{\left({2 \ {{x}^{2}}}+{2 \ {{b}^{2}}}\right)}\ {\sqrt{{{x}^{2}}+{{b}^{2}}}}}$

(11)

Type: Expression(Complex(Integer))

fricas

--
realUnitDoublet:=eval(eval(abs''',z=a+%i*b),[conjugate(a)=a,conjugate(b)=b])

$\label{eq12}-{{3 \ a \ {{b}^{2}}}\over{{\left({{b}^{4}}+{2 \ {{a}^{2}}\ {{b}^{2}}}+{{a}^{4}}\right)}\ {\sqrt{{{b}^{2}}+{{a}^{2}}}}}}$

(12)

Type: Expression(Complex(Integer))

fricas

test(realUnitDoublet=D(abz(a),[a,a,a]))

$\label{eq13} \mbox{\rm true}$

(13)

Type: Boolean

fricas

unitDoublet(x)==eval(-realUnitDoublet/2,a=x); unitDoublet(x)

fricas

Compiling function unitDoublet with type Variable(x) -> Expression(
      Complex(Integer))

$\label{eq14}{3 \ {{b}^{2}}\ x}\over{{\left({2 \ {{x}^{4}}}+{4 \ {{b}^{2}}\ {{x}^{2}}}+{2 \ {{b}^{4}}}\right)}\ {\sqrt{{{x}^{2}}+{{b}^{2}}}}}$

(14)

Type: Expression(Complex(Integer))

We are interested in the behaviour of these functions in the limit as $b \to 0$ .

signum Properties

fricas

limit( signum(1) ,b=0)

fricas

Compiling function signum with type PositiveInteger -> Expression(
      Complex(Integer))

$\label{eq15}1$

(15)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( signum(-1) ,b=0)

fricas

Compiling function signum with type Integer -> Expression(Complex(
      Integer))

$\label{eq16}- 1$

(16)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( signum(0) ,b=0)

fricas

Compiling function signum with type NonNegativeInteger -> Expression
      (Complex(Integer))

$\label{eq17}0$

(17)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( signum(x)^2 ,b=0)

$\label{eq18}1$

(18)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

-- odd function
test( signum(x)=-signum(-x))

fricas

Compiling function signum with type Polynomial(Integer) -> 
      Expression(Complex(Integer))

$\label{eq19} \mbox{\rm true}$

(19)

Type: Boolean

diracDelta Properties

http://en.wikipedia.org/wiki/Dirac_delta_function

fricas

limit( diracDelta(1) ,b=0)

fricas

Compiling function diracDelta with type PositiveInteger -> 
      Expression(Complex(Integer))

$\label{eq20}0$

(20)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( diracDelta(-1) ,b=0)

fricas

Compiling function diracDelta with type Integer -> Expression(
      Complex(Integer))

$\label{eq21}0$

(21)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( diracDelta(0) ,b=0)

fricas

Compiling function diracDelta with type NonNegativeInteger -> 
      Expression(Complex(Integer))

$\label{eq22}+ \infty$

(22)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

-- even function
test( diracDelta(x)=diracDelta(-x) )

fricas

Compiling function diracDelta with type Polynomial(Integer) -> 
      Expression(Complex(Integer))

$\label{eq23} \mbox{\rm true}$

(23)

Type: Boolean

fricas

limit( diracDelta(x)^2 ,b=0)

$\label{eq24}0$

(24)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

integrate(diracDelta(x),x=%minusInfinity..%plusInfinity,"noPole")

$\label{eq25}1$

(25)

Type: Union(f1: OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

f(x)==c1*x+c2

Type: Void

fricas

test(integrate(f(x)*diracDelta(x-t),x=%minusInfinity..%plusInfinity,"noPole")=f(t))

fricas

Compiling function f with type Variable(x) -> Polynomial(Integer)

fricas

Compiling function f with type Variable(t) -> Polynomial(Integer)

$\label{eq26} \mbox{\rm true}$

(26)

Type: Boolean

fricas

limit(integrate(log(x)*diracDelta(x-t),x=%minusInfinity..%plusInfinity,"noPole"), b=0)

$\label{eq27}{{\log \left({4 \ {{t}^{2}}}\right)}-{\log \left({1 \over{4 \ {{t}^{2}}}}\right)}-{2 \ {\log \left({4}\right)}}}\over 4$

(27)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

complexNormalize complexExpand %

$\label{eq28}{-{\log \left({1 \over{4 \ {{t}^{2}}}}\right)}-{\log \left({4}\right)}}\over 2$

(28)

Type: Expression(Complex(Integer))

fricas

eval(%,t=exp(x))

$\label{eq29}{-{\log \left({1 \over{4 \ {{{e}^{x}}^{2}}}}\right)}-{\log \left({4}\right)}}\over 2$

(29)

Type: Expression(Complex(Integer))

fricas

integrate(diracDelta(2*x),x=%minusInfinity..%plusInfinity,"noPole")

$\label{eq30}1 \over 2$

(30)

Type: Union(f1: OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

--
limit( signum(x)*diracDelta(x) ,b=0)

$\label{eq31}0$

(31)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

diracDelta(x^2)

$\label{eq32}{{b}^{2}}\over{{\left({2 \ {{x}^{4}}}+{2 \ {{b}^{2}}}\right)}\ {\sqrt{{{x}^{4}}+{{b}^{2}}}}}$

(32)

Type: Expression(Complex(Integer))

fricas

limit(%,b=0)

$\label{eq33}0$

(33)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

The Unit Doublet function comes after diracDelta

http://en.wikipedia.org/wiki/Unit_doublet

fricas

test(unitDoublet(x)=3*diracDelta(x)*signum(x)/abz(x))

$\label{eq34} \mbox{\rm true}$

(34)

Type: Boolean

fricas

limit( unitDoublet(1) ,b=0)

fricas

Compiling function unitDoublet with type PositiveInteger -> 
      Expression(Complex(Integer))

$\label{eq35}0$

(35)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( unitDoublet(-1) ,b=0)

fricas

Compiling function unitDoublet with type Integer -> Expression(
      Complex(Integer))

$\label{eq36}0$

(36)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( unitDoublet(0) ,b=0)

fricas

Compiling function unitDoublet with type NonNegativeInteger -> 
      Expression(Complex(Integer))

$\label{eq37}0$

(37)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

-- odd function
test( unitDoublet(x)=-unitDoublet(-x) )

fricas

Compiling function unitDoublet with type Polynomial(Integer) -> 
      Expression(Complex(Integer))

$\label{eq38} \mbox{\rm true}$

(38)

Type: Boolean

fricas

limit( unitDoublet(x)^2 ,b=0)

$\label{eq39}0$

(39)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

integrate(unitDoublet(x),x=%minusInfinity..%plusInfinity,"noPole")

$\label{eq40}0$

(40)

Type: Union(f1: OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

f(x)==c1*x^2+c2*x+c3

   Compiled code for f has been cleared.
   1 old definition(s) deleted for function or rule f

Type: Void

fricas

test(integrate(f(x)*unitDoublet(x-t),x=%minusInfinity..%plusInfinity,"noPole")=eval(D(f(x),x),x=t))

fricas

Compiling function f with type Variable(x) -> Polynomial(Integer)

$\label{eq41} \mbox{\rm true}$

(41)

Type: Boolean

fricas

integrate(unitDoublet(2*x),x=%minusInfinity..%plusInfinity,"noPole")

$\label{eq42}0$

(42)

Type: Union(f1: OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

--
limit( signum(x)*unitDoublet(x) ,b=0)

$\label{eq43}0$

(43)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( diracDelta(x)*unitDoublet(x) ,b=0)

$\label{eq44}0$

(44)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

unitDoublet(x^2)

$\label{eq45}{3 \ {{b}^{2}}\ {{x}^{2}}}\over{{\left({2 \ {{x}^{8}}}+{4 \ {{b}^{2}}\ {{x}^{4}}}+{2 \ {{b}^{4}}}\right)}\ {\sqrt{{{x}^{4}}+{{b}^{2}}}}}$

(45)

Type: Expression(Complex(Integer))

fricas

limit(%,b=0)

$\label{eq46}0$

(46)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

Heaviside

http://en.wikipedia.org/wiki/Unit_step_function

fricas

realHeaviside:=integrate(diracDelta(a),a)::Expression Complex Integer

fricas

Compiling function diracDelta with type Variable(a) -> Expression(
      Complex(Integer))

$\label{eq47}-{{{b}^{2}}\over{{2 \ a \ {\sqrt{{{b}^{2}}+{{a}^{2}}}}}-{2 \ {{b}^{2}}}-{2 \ {{a}^{2}}}}}$

(47)

Type: Expression(Complex(Integer))

fricas

heavisideStep(x)==eval(realHeaviside,a=x)

Type: Void

fricas

limit( heavisideStep(1), b=0)

fricas

Compiling function heavisideStep with type PositiveInteger -> 
      Expression(Complex(Integer))

$\label{eq48}1$

(48)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( heavisideStep(-1) ,b=0)

fricas

Compiling function heavisideStep with type Integer -> Expression(
      Complex(Integer))

$\label{eq49}0$

(49)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( heavisideStep(0) ,b=0)

fricas

Compiling function heavisideStep with type NonNegativeInteger -> 
      Expression(Complex(Integer))

$\label{eq50}1 \over 2$

(50)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( heavisideStep(x), x=0)

fricas

Compiling function heavisideStep with type Variable(x) -> Expression
      (Complex(Integer))

$\label{eq51}1 \over 2$

(51)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( heavisideStep(x) ,x=%plusInfinity)

$\label{eq52}1$

(52)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit( heavisideStep(x) ,x=%minusInfinity)

$\label{eq53}0$

(53)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

test(simplify( (1+signum(x))/2 - heavisideStep(x) ) = 0)

$\label{eq54} \mbox{\rm true}$

(54)

Type: Boolean

fricas

test(limit( integrate(heavisideStep(t),t=%minusInfinity..x, "noPole") - x*heavisideStep(x), b=0) =0)

fricas

Compiling function heavisideStep with type Variable(t) -> Expression
      (Complex(Integer))

$\label{eq55} \mbox{\rm true}$

(55)

Type: Boolean

fricas

-- diracDelta?
d1:=integrate(heavisideStep(t)*unitDoublet(t),t=%minusInfinity..x, "noPole")

fricas

Compiling function unitDoublet with type Variable(t) -> Expression(
      Complex(Integer))

$\label{eq56}{\left( \begin{array}{@{}l} \displaystyle {{\left({ \begin{array}{@{}l} \displaystyle { \begin{array}{@{}l} \displaystyle {\left({{{96}\ i \ {{x}^{7}}}+{{264}\ i \ {{b}^{2}}\ {{x}^{5}}}+{{240}\ i \ {{b}^{4}}\ {{x}^{3}}}+{{72}\ i \ {{b}^{6}}\ x}}\right)}\ \cdot \ \ \displaystyle {\sqrt{1 \over{{1024}\ {{b}^{2}}}}}\ \cdot \ \ \displaystyle {\sqrt{{{x}^{2}}+{{b}^{2}}}}$

(56)

Type: Union(f1: OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

limit(d1,b=0)

$\label{eq57}0$

(57)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

normalize(d1-diracDelta(x))

$\label{eq58}{\left( \begin{array}{@{}l} \displaystyle {{\left({ \begin{array}{@{}l} \displaystyle {{\left({ \begin{array}{@{}l} \displaystyle {{12}\ i \ {{x}^{6}}}+ \ \ \displaystyle {{27}\ i \ {{b}^{2}}\ {{x}^{4}}}+ \ \ \displaystyle {{18}\ i \ {{b}^{4}}\ {{x}^{2}}}+ \ \ \displaystyle {3 \ i \ {{b}^{6}}}$

(58)

Type: Expression(Complex(Integer))

Problems?

fricas

-- expected 1/abs(c)
integrate(diracDelta(c*x),x=%minusInfinity..%plusInfinity,"noPole")

$\label{eq59}0$

(59)

Type: Union(f1: OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

f(x)==c1*x^2+c2*x+c3

   Compiled code for f has been cleared.
   1 old definition(s) deleted for function or rule f

Type: Void

fricas

integrate(f(x)*diracDelta(x-t),x=%minusInfinity..%plusInfinity,"noPole")

fricas

Compiling function f with type Variable(x) -> Polynomial(Integer)

$\label{eq60}\mbox{\tt "failed"}$

(60)

Type: Union(fail: failed,...)

fricas

-- expected
f(t)

fricas

Compiling function f with type Variable(t) -> Polynomial(Integer)

$\label{eq61}{c 1 \ {{t}^{2}}}+{c 2 \ t}+ c 3$

(61)

Type: Polynomial(Integer)

fricas

-- expected 1
integrate(exp(x)*diracDelta(x),x=%minusInfinity..%plusInfinity,"noPole")

$\label{eq62}\mbox{\tt "failed"}$

(62)

Type: Union(fail: failed,...)

fricas

-- expected 1
integrate(exp(x)*unitDoublet(x),x=%minusInfinity..%plusInfinity,"noPole")

$\label{eq63}\mbox{\tt "failed"}$

(63)

Type: Union(fail: failed,...)

fricas

-- loops forever on %plusInfinity limit of indefinite integral
--integrate(log(x)*unitDoublet(x-t),x=%minusInfinity..%plusInfinity,"noPole")
-- expected x*abs(x)/2
integrate(abs(x),x)

fricas

Compiling function abs with type Variable(x) -> Expression(Integer)

$\label{eq64}\int^{ \displaystyle x}{{\sqrt{\%A \ {\overline \%A}}}\ {d \%A}}$

(64)

Type: Union(Expression(Integer),...)

fricas

integrate(abz(x),x)

$\label{eq65}{\left( \begin{array}{@{}l} \displaystyle {{\left({ \begin{array}{@{}l} \displaystyle -{2 \ {{b}^{2}}\ x \ {\sqrt{{{x}^{2}}+{{b}^{2}}}}}+ \ \ \displaystyle {2 \ {{b}^{2}}\ {{x}^{2}}}+{{b}^{4}}$

(65)

Type: Union(Expression(Complex(Integer)),...)

fricas

limit(%,b=0)

$\label{eq66}-{{{x}^{2}}\over 2}$

(66)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

-- expected abs(x)
integrate(abs',z)

$\label{eq67}\int^{ \displaystyle z}{{{{\overline \%A}+ \%A}\over{2 \ {\sqrt{\%A \ {\overline \%A}}}}}\ {d \%A}}$

(67)

Type: Union(Expression(Integer),...)

fricas

integrate(signum(x),x)

$\label{eq68}{-{x \ {\sqrt{{{x}^{2}}+{{b}^{2}}}}}+{{x}^{2}}+{{b}^{2}}}\over{{\sqrt{{{x}^{2}}+{{b}^{2}}}}- x}$

(68)

Type: Union(Expression(Complex(Integer)),...)

fricas

limit(%,b=0)

$\label{eq69}- x$

(69)

Type: Union(OrderedCompletion?(Expression(Complex(Integer))),...)

fricas

-- expected signum(x)/2
integrate(abs'',z)

$\label{eq70}\int^{ \displaystyle z}{{{-{{\overline \%A}^{2}}+{2 \ \%A \ {\overline \%A}}-{{\%A}^{2}}}\over{4 \ \%A \ {\overline \%A}\ {\sqrt{\%A \ {\overline \%A}}}}}\ {d \%A}}$

(70)

Type: Union(Expression(Integer),...)

fricas

integrate(diracDelta(x^2+y^2), x=%minusInfinity..%plusInfinity, "noPole")

$\label{eq71}\mbox{\tt "failed"}$

(71)

Type: Union(fail: failed,...)

Convolution

fricas

g(x)==abz(x)/2; g(y)

fricas

Compiling function abz with type Variable(y) -> Expression(Complex(
      Integer))

fricas

Compiling function g with type Variable(y) -> Expression(Complex(
      Integer))

$\label{eq72}{\sqrt{{{y}^{2}}+{{b}^{2}}}}\over 2$

(72)

Type: Expression(Complex(Integer))

fricas

h(x)==heavisideStep(x+1/2)-heavisideStep(x-1/2); h(x)

fricas

Compiling function heavisideStep with type Polynomial(Fraction(
      Integer)) -> Expression(Complex(Integer))

fricas

Compiling function h with type Variable(x) -> Expression(Complex(
      Integer))

$\label{eq73}{\left( \begin{array}{@{}l} \displaystyle {{\left({4 \ {{b}^{2}}\ x}+{2 \ {{b}^{2}}}\right)}\ {\sqrt{{4 \ {{x}^{2}}}+{4 \ x}+{4 \ {{b}^{2}}}+ 1}}}+ \ \ \displaystyle {{\left(-{4 \ {{b}^{2}}\ x}+{2 \ {{b}^{2}}}\right)}\ {\sqrt{{4 \ {{x}^{2}}}-{4 \ x}+{4 \ {{b}^{2}}}+ 1}}}-{{16}\ {{b}^{2}}\ x}$

(73)

Type: Expression(Complex(Integer))

fricas

conv:=integrate(h(x-y)*g(y),y=%minusInfinity..%plusInfinity,"noPole")

fricas

Compiling function h with type Polynomial(Integer) -> Expression(
      Complex(Integer))

$\label{eq74}\mbox{\tt "failed"}$

(74)

Type: Union(fail: failed,...)

Representing a distribution as a series of bump functions

fricas

)lib GDRAW

   GnuDraw is now explicitly exposed in frame initial 
   GnuDraw will be automatically loaded when needed from 
      /var/aw/var/LatexWiki/GDRAW.NRLIB/GDRAW
X:=[(x/10)::DFLOAT for x in -100..100 by 1];

Type: List(DoubleFloat?)

fricas

Y:=[eval(unitDoublet(x),b=1.0)::DFLOAT for x in X];

fricas

Compiling function unitDoublet with type DoubleFloat -> Expression(
      Complex(DoubleFloat))

Type: List(DoubleFloat?)

fricas

gnuDraw(X,Y,"SandBoxDiracDeltaX2.dat")

   Graph data being transmitted to the viewport manager...
   FriCAS2D data being transmitted to the viewport manager...

Type: Void

[terminal=pslatex,terminaloptions=color,scale=1.3]
load "SandBoxDiracDeltaX2.dat"