 | Find problems with this code (comp is a Fortran function in the extension module pde):
x = arange(0, 1, 0.01)
b = myfunc1(x) # compute b array of size (n,n)
u = myfunc2(x) # compute u array of size (n,n)
c = myfunc3(x) # compute c array of size (n,n)
dt = 0.05
for i in range(n)
u = pde.comp(u, b, c, i*dt)
|
