I have to solve the following problem :

enter image description here

I wrote the second hand :

def dN_dthyb(t,N,param):

dN = np.zeros(m)
dN[0] = r0-alpha[0]*N[0]
for i in range(1, m-1):
    dN[i] = alpha[i-1]*N[i-1]-alpha[i]*N[i]
dN[m-1] = alpha[m-2]*N[m-2]
return dN

Where param is an array with the values of alpha and the value of r(t) (here r is constant) I already used this function in order to solve the problem as a Cauchy problem with ode and now I want to have the boundary coundition such as N(0) = N(T), and r is a periodic function with amplitude equal to r0.

First, in order to use solve_bvp, I rewrore the function :

def dN_dthyb(t,N,param):
dN = np.zeros((m+n,len(t)))
for i in range(1, m-1):
    dN[i,:] = alpha[i-1]*N[i-1,:]-alpha[i]*N[i,:]
dN[m-1,:] = alpha[m-2]*N[m-2,:]
dN[m,:] = -beta[0]*N[m,:]
for i in range(1, n-1):
return dN

Is it correct to see N as a matrix where each row, for example row i represents N_i(t) ?

Secondly, I don't understand how to set the boundary condition. I don't understand the nature of the return value. What does this value correspond to?

Thanks for your help,

  • $\begingroup$ from scipy.integrate import solve_bvp, note the "from". Is there a reason to suspect that periodic solutions exist? You would need strong forcing terms $p(t)$, $r(t)$ to overcome the mostly exponential behavior from the remaining terms $\endgroup$ Aug 4 at 14:47
  • $\begingroup$ Please clarify your specific problem or provide additional details to highlight exactly what you need. As it's currently written, it's hard to tell exactly what you're asking. $\endgroup$
    – Community Bot
    Aug 4 at 14:49
  • $\begingroup$ @Lutz Lehmann yes, I forgot to specify that in this case, r is a function, I'm sorry $\endgroup$
    – Insolence
    Aug 4 at 15:23
  • $\begingroup$ As the equation is affine-linear, the propagator from 0 to T will also be affine-linear. So from a sufficient number of independent solutions you could construct a (convex, weights sum to one) linear combination that solves the boundary conditions. $\endgroup$ Aug 4 at 15:46


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.