So this is the function I have for like a standing wave solution to generate the current u:
def standing_wave_solution(t,nx,mx,my,x,y): dx = (1/nx-1)