02-152824.py

from nutils import mesh, function as fn, plot
import numpy as np


xpts = np.linspace(0, 4, 5)
domain, geom = mesh.rectilinear([xpts])
basis = domain.basis('spline', 3, knotmultiplicities=[[4, 1, 3, 1, 4]])
dofs = [0, 0, 0, .5, 1, .5, 0, 0, 0]
func = basis.dot(dofs)

x, y = domain.elem_eval([geom, func], separate=True, ischeme='bezier40')
with plot.PyPlot('test', index=0) as plt:
    plt.mesh(x, y)
    plt.xlim(0, 4)