function y=pr03b_fce(x)
y=x.^5-9.*x.^3+5-2.*x;