function y=rhsfu(x)
y=100*exp(-10*x);
