function y = fum_temp(x)

y = x^2 + sin(x);