Both the equation defining a PDE and its boundary conditions involve functions that are the result of the numerical resolution process of an ODE system.

The PDE solution is obtained. However, when I want to calculate the value of the solution at some point the following error message is displayed: "Error, (in solnproc) unable to compute solution for t>HFloat(1.8936825887327868e15): unable to store HFloat(undefined)+HFloat(undefined)*I when datatype=float[8]".

I tried to formulate the problem differently. However, the error appeared again.

Where am I going wrong?

I've attached the PDE.mw file to the conference. 

Thank you for your help.

Is it possible to solve numerically a PDE, where the IBC is a function resulting from the numerical resolution of an ODE?

This differential equation has an analytical solution.

However, I am looking for a numerical solution for it.

A single boundary condition with respect to t would be sufficient to solve the problem, but the command does not accept this.

In the end, I try to insert boundary conditions in relation to the variables x and t, but again this does not work.

Where am I going wrong?