code: def linear_multistep_coeff(order, t, i, j):
      if order - 1 > i:
          raise ValueError(f"Order {order} too high for step {i}")
  
      def fn(tau):
          prod = 1.0
          for k in range(order):
              if j == k:
                  continue
              prod *= (tau - t[i - k]) / (t[i - j] - t[i - k])
          return prod
  
      return integrate.quad(fn, t[i], t[i + 1], epsrel=1e-4)[0]
  
imports: import math
  import torch
  import torchsde
  from scipy import integrate
  from torch import nn
  from torchdiffeq import odeint
  from tqdm.auto import tqdm, trange
label_desc: Use image node_id 15 for calling a global function with example usage: linear_multistep_coeff(order, t, i, j) and returns: integrate, prod