function u = rho(x)
    u = sin(pi*x).^4;
end
