function xpr = test(t, x)
xpr = cos(t) / (2 * (x - 1));

