function y=rhsfu(x)
y=-exp(x-1);
