Skip to main content
Mathematica

Return to Question

boarding -> boundary
Source Link
flinty
flinty
  • 25.9k
  • 2
  • 22
  • 92

I have the following code:

sys = {
    D[f[ϵ, x], {x, 2}] + Iee[ϵ, x] + Iph[ϵ, x] == 0,  
    f[ϵ, -L/2] == fF[ϵ - eV/2], 
    f[ϵ, L/2] == fF[ϵ + eV/2]
}
where:
Iee[ϵ_, x_] := a Integrate[
    Integrate[
        (f[ϵ, x] f[ϵ1 - ω, x] (1 - f[ϵ - ω, x]) (1 - f[ϵ1, x]) - 
         f[ϵ - ω, x] f[ϵ1, x] (1 - f[ϵ, x]) (1 - f[ϵ1 - ω, x])), 
        {ϵ1, -∞, ∞}, 
        {ω, -∞, ∞}
    ]
]

Iph[ϵ_, x_] := b Integrate[
    ω^2 ((1 - f[ϵ, x]) f[ϵ + ω, x] (1 + n[ω]) + 
         (1 - f[ϵ, x]) f[ϵ - ω] n[ω] - 
         f[ϵ, x] (1 - f[ϵ - ω]) (1 + n[ω]) - 
         f[ϵ, x] (1 - f[ϵ + ω]) n[ω]), 
    {ω, 0, ∞}
]
n[ω_] = 1/(Exp[ω/T] - 1);
fF[ϵ_] = 1/(1 + Exp[ϵ/T]);

where a,b,T,L,eV are constants.

I need to solve this equation numerically. As I understand, simple functions as well as NDSolve don't work here, and I need to use The Numerical Method of Lines (like here Integrating over $x$ in numerically solving a partial integrodifferential equation)

But anyway I still dont understand what to do with limits of integration and how to connect it with my boardingboundary conditions.

I have the following code:

sys = {
    D[f[ϵ, x], {x, 2}] + Iee[ϵ, x] + Iph[ϵ, x] == 0,  
    f[ϵ, -L/2] == fF[ϵ - eV/2], 
    f[ϵ, L/2] == fF[ϵ + eV/2]
}
where:
Iee[ϵ_, x_] := a Integrate[
    Integrate[
        (f[ϵ, x] f[ϵ1 - ω, x] (1 - f[ϵ - ω, x]) (1 - f[ϵ1, x]) - 
         f[ϵ - ω, x] f[ϵ1, x] (1 - f[ϵ, x]) (1 - f[ϵ1 - ω, x])), 
        {ϵ1, -∞, ∞}, 
        {ω, -∞, ∞}
    ]
]

Iph[ϵ_, x_] := b Integrate[
    ω^2 ((1 - f[ϵ, x]) f[ϵ + ω, x] (1 + n[ω]) + 
         (1 - f[ϵ, x]) f[ϵ - ω] n[ω] - 
         f[ϵ, x] (1 - f[ϵ - ω]) (1 + n[ω]) - 
         f[ϵ, x] (1 - f[ϵ + ω]) n[ω]), 
    {ω, 0, ∞}
]
n[ω_] = 1/(Exp[ω/T] - 1);
fF[ϵ_] = 1/(1 + Exp[ϵ/T]);

where a,b,T,L,eV are constants.

I need to solve this equation numerically. As I understand, simple functions as well as NDSolve don't work here, and I need to use The Numerical Method of Lines (like here Integrating over $x$ in numerically solving a partial integrodifferential equation)

But anyway I still dont understand what to do with limits of integration and how to connect it with my boarding conditions.

I have the following code:

sys = {
    D[f[ϵ, x], {x, 2}] + Iee[ϵ, x] + Iph[ϵ, x] == 0,  
    f[ϵ, -L/2] == fF[ϵ - eV/2], 
    f[ϵ, L/2] == fF[ϵ + eV/2]
}
where:
Iee[ϵ_, x_] := a Integrate[
    Integrate[
        (f[ϵ, x] f[ϵ1 - ω, x] (1 - f[ϵ - ω, x]) (1 - f[ϵ1, x]) - 
         f[ϵ - ω, x] f[ϵ1, x] (1 - f[ϵ, x]) (1 - f[ϵ1 - ω, x])), 
        {ϵ1, -∞, ∞}, 
        {ω, -∞, ∞}
    ]
]

Iph[ϵ_, x_] := b Integrate[
    ω^2 ((1 - f[ϵ, x]) f[ϵ + ω, x] (1 + n[ω]) + 
         (1 - f[ϵ, x]) f[ϵ - ω] n[ω] - 
         f[ϵ, x] (1 - f[ϵ - ω]) (1 + n[ω]) - 
         f[ϵ, x] (1 - f[ϵ + ω]) n[ω]), 
    {ω, 0, ∞}
]
n[ω_] = 1/(Exp[ω/T] - 1);
fF[ϵ_] = 1/(1 + Exp[ϵ/T]);

where a,b,T,L,eV are constants.

I need to solve this equation numerically. As I understand, simple functions as well as NDSolve don't work here, and I need to use The Numerical Method of Lines (like here Integrating over $x$ in numerically solving a partial integrodifferential equation)

But anyway I still dont understand what to do with limits of integration and how to connect it with my boundary conditions.

deleted 9 characters in body; edited title
Source Link
Domen
Domen
  • 33.4k
  • 3
  • 47
  • 66

Numerically Numerical solution jfof integro-differential equation where the integral term contains a double integral over one of the variables in infinite limits

So, i have:

I have the following code:

sys = {
    D[f[ϵ, x], {x, 2}] + Iee[ϵ, x] + Iph[ϵ, x] == 0,  
    f[ϵ, -L/2] == fF[ϵ - eV/2], 
    f[ϵ, L/2] == fF[ϵ + eV/2]
}
where:
Iee[ϵ_, x_] := a Integrate[
    Integrate[
        (f[ϵ, x] f[ϵ1 - ω, x] (1 - f[ϵ - ω, x]) (1 - f[ϵ1, x]) - 
         f[ϵ - ω, x] f[ϵ1, x] (1 - f[ϵ, x]) (1 - f[ϵ1 - ω, x])), 
        {ϵ1, -∞, ∞}, 
        {ω, -∞, ∞}
    ]
]

Iph[ϵ_, x_] := b Integrate[
    ω^2 ((1 - f[ϵ, x]) f[ϵ + ω, x] (1 + n[ω]) + 
         (1 - f[ϵ, x]) f[ϵ - ω] n[ω] - 
         f[ϵ, x] (1 - f[ϵ - ω]) (1 + n[ω]) - 
         f[ϵ, x] (1 - f[ϵ + ω]) n[ω]), 
    {ω, 0, ∞}
]
n[ω_] = 1/(Exp[ω/T] - 1);
fF[ϵ_] = 1/(1 + Exp[ϵ/T]);

a,b,T,L,eV - constants.

I need to solve this equation numerically. As i understand, simple functions as well as NDSolve don't work here, and i need to use The Numerical Method of Lines (like here Integrating over $x$ in numerically solving a partial integrodifferential equation) But anyway i still dont understand what to do with limits of integration and how to connect it with my boarding conditions. Help please.

where a,b,T,L,eV are constants.

I need to solve this equation numerically. As I understand, simple functions as well as NDSolve don't work here, and I need to use The Numerical Method of Lines (like here Integrating over $x$ in numerically solving a partial integrodifferential equation)

But anyway I still dont understand what to do with limits of integration and how to connect it with my boarding conditions.

Numerically solution jf integro-differential equation where the integral term contains a double integral over one of the variables in infinite limits

So, i have:

sys = {
    D[f[ϵ, x], {x, 2}] + Iee[ϵ, x] + Iph[ϵ, x] == 0,  
    f[ϵ, -L/2] == fF[ϵ - eV/2], 
    f[ϵ, L/2] == fF[ϵ + eV/2]
}
where:
Iee[ϵ_, x_] := a Integrate[
    Integrate[
        (f[ϵ, x] f[ϵ1 - ω, x] (1 - f[ϵ - ω, x]) (1 - f[ϵ1, x]) - 
         f[ϵ - ω, x] f[ϵ1, x] (1 - f[ϵ, x]) (1 - f[ϵ1 - ω, x])), 
        {ϵ1, -∞, ∞}, 
        {ω, -∞, ∞}
    ]
]

Iph[ϵ_, x_] := b Integrate[
    ω^2 ((1 - f[ϵ, x]) f[ϵ + ω, x] (1 + n[ω]) + 
         (1 - f[ϵ, x]) f[ϵ - ω] n[ω] - 
         f[ϵ, x] (1 - f[ϵ - ω]) (1 + n[ω]) - 
         f[ϵ, x] (1 - f[ϵ + ω]) n[ω]), 
    {ω, 0, ∞}
]
n[ω_] = 1/(Exp[ω/T] - 1);
fF[ϵ_] = 1/(1 + Exp[ϵ/T]);

a,b,T,L,eV - constants.

I need to solve this equation numerically. As i understand, simple functions as well as NDSolve don't work here, and i need to use The Numerical Method of Lines (like here Integrating over $x$ in numerically solving a partial integrodifferential equation) But anyway i still dont understand what to do with limits of integration and how to connect it with my boarding conditions. Help please.

Numerical solution of integro-differential equation where the integral term contains a double integral

I have the following code:

sys = {
    D[f[ϵ, x], {x, 2}] + Iee[ϵ, x] + Iph[ϵ, x] == 0,  
    f[ϵ, -L/2] == fF[ϵ - eV/2], 
    f[ϵ, L/2] == fF[ϵ + eV/2]
}
where:
Iee[ϵ_, x_] := a Integrate[
    Integrate[
        (f[ϵ, x] f[ϵ1 - ω, x] (1 - f[ϵ - ω, x]) (1 - f[ϵ1, x]) - 
         f[ϵ - ω, x] f[ϵ1, x] (1 - f[ϵ, x]) (1 - f[ϵ1 - ω, x])), 
        {ϵ1, -∞, ∞}, 
        {ω, -∞, ∞}
    ]
]

Iph[ϵ_, x_] := b Integrate[
    ω^2 ((1 - f[ϵ, x]) f[ϵ + ω, x] (1 + n[ω]) + 
         (1 - f[ϵ, x]) f[ϵ - ω] n[ω] - 
         f[ϵ, x] (1 - f[ϵ - ω]) (1 + n[ω]) - 
         f[ϵ, x] (1 - f[ϵ + ω]) n[ω]), 
    {ω, 0, ∞}
]
n[ω_] = 1/(Exp[ω/T] - 1);
fF[ϵ_] = 1/(1 + Exp[ϵ/T]);

where a,b,T,L,eV are constants.

I need to solve this equation numerically. As I understand, simple functions as well as NDSolve don't work here, and I need to use The Numerical Method of Lines (like here Integrating over $x$ in numerically solving a partial integrodifferential equation)

But anyway I still dont understand what to do with limits of integration and how to connect it with my boarding conditions.

Source Link
Pablo Escobar
Pablo Escobar
  • 1
  • 1

Numerically solution jf integro-differential equation where the integral term contains a double integral over one of the variables in infinite limits

So, i have:

sys = {
    D[f[ϵ, x], {x, 2}] + Iee[ϵ, x] + Iph[ϵ, x] == 0,  
    f[ϵ, -L/2] == fF[ϵ - eV/2], 
    f[ϵ, L/2] == fF[ϵ + eV/2]
}
where:
Iee[ϵ_, x_] := a Integrate[
    Integrate[
        (f[ϵ, x] f[ϵ1 - ω, x] (1 - f[ϵ - ω, x]) (1 - f[ϵ1, x]) - 
         f[ϵ - ω, x] f[ϵ1, x] (1 - f[ϵ, x]) (1 - f[ϵ1 - ω, x])), 
        {ϵ1, -∞, ∞}, 
        {ω, -∞, ∞}
    ]
]

Iph[ϵ_, x_] := b Integrate[
    ω^2 ((1 - f[ϵ, x]) f[ϵ + ω, x] (1 + n[ω]) + 
         (1 - f[ϵ, x]) f[ϵ - ω] n[ω] - 
         f[ϵ, x] (1 - f[ϵ - ω]) (1 + n[ω]) - 
         f[ϵ, x] (1 - f[ϵ + ω]) n[ω]), 
    {ω, 0, ∞}
]
n[ω_] = 1/(Exp[ω/T] - 1);
fF[ϵ_] = 1/(1 + Exp[ϵ/T]);

a,b,T,L,eV - constants.

I need to solve this equation numerically. As i understand, simple functions as well as NDSolve don't work here, and i need to use The Numerical Method of Lines (like here Integrating over $x$ in numerically solving a partial integrodifferential equation) But anyway i still dont understand what to do with limits of integration and how to connect it with my boarding conditions. Help please.