MATLAB中u1([2:M])= sin(pi*x)+x.*(1 -