# How to use a convolution integral with NeuralPDE.discretize

**URL:** https://discourse.julialang.org/t/how-to-use-a-convolution-integral-with-neuralpde-discretize/137997
**Category:** Modelling & Simulations
**Created:** [July 7, 2026, 10:53pm UTC](https://discourse.julialang.org/t/how-to-use-a-convolution-integral-with-neuralpde-discretize/137997 "2026-07-07T22:53:27Z")
**Posts on this page:** 1
**Showing post:** 3

<div class="post-metadata">

### Author: ![cmichelenstrofer](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/cmichelenstrofer/32/44392_2.png) [@cmichelenstrofer](https://discourse.julialang.org/u/cmichelenstrofer)
#### Post date: [July 8, 2026, 3:53pm UTC](https://discourse.julialang.org/t/how-to-use-a-convolution-integral-with-neuralpde-discretize/137997/3 "2026-07-08T15:53:07Z")

</div>

[Similar issue](https://discourse.julialang.org/t/way-to-write-the-integral-with-modelingtoolkit/82347/5) from a few years ago, but doesn’t seem to have a resolution either.

> [@Way to write the integral with ModelingToolkit](https://discourse.julialang.org/t/way-to-write-the-integral-with-modelingtoolkit/82347/5):
>
> Hi @ChrisRackauckas I tried what you suggested but I am unable to discretize the equation, below is the code and the error. I am not sure where to pass the parameter x. using NeuralPDE, Flux, ModelingToolkit, DiffEqFlux, DomainSets import ModelingToolkit: Interval @parameters t, x @variables u(..) lambda = -1.5 Di = Differential(t) Ii = Integral(x in DomainSets.ClosedInterval(0, t)) eq = Di(u(t)) + Ii((t-x)^lambda\*u(x)) ~ 1.0 bcs = [u(0.0) ~ 0.0] domains = [t ∈ Interval(0.0,2.0)] chain = Chain…

---

_[View the full topic](https://discourse.julialang.org/t/how-to-use-a-convolution-integral-with-neuralpde-discretize/137997)._