in main:

    nx = 1066 #100  # NUMBER OF X-STEPS
    nt = 2000 #1500  # NUMBER OF TIME STEPS
    endtime = 0.312
    width = 1
    method = ut.LaxWendroff


INITIAL STATE:
    if initialConditions =="sine3":
        for i in range(nx):
            V[0][i][0] += 0.01*np.sin(4*np.pi*i/nx)
            V[0][i][1] += 0.002*np.sin(np.pi*i/nx)