function ex403
    [x,y]=ode45(@f,[0,pi],0);
    plot(x,y,'o-');grid('on');
function xy=f(x)
    xy=cos(x);