fig snapshots of the dimensionless solution of scaled problem rhs[ : - (beta/dx** )*( [ : + * [ :nu[ : - ] ( [ : ]tthis rewrite speeds up the code by about factor of complete code is found in the file rod_fe_vec py using odespy to solve the system of odes let us now show how to apply general ode package like odespy (see sect to solve our diffusion problem as long as we have defined right-hand side function rhs this is very straightforwardimport odespy solver odespy rkfehlberg(rhssolver set_initial_condition(u_ n_t int(round( /float(dt))time_points linspace( tn_t+ ut solver solve(time_points