function lineq
    grid on
    x0=0;y0=1;
    [T,y]=ode23s(@ff,[x0,x0+10],y0);
    steps=length(T)-1
    plot(T,y,'-o')
    function z=ff(x,y)
        z=-20*y+20*x+21;
    end
   
end

