exactsol = Simplify[
k[t]/.DSolve[{k'[t]==s ((k[t])^α) -(δ+n) k[t],k[0]==a},k[t],t][[1]]]

微分方程式から導かれるベクトル場を描いてみる．

VectorPlot[{1,s (k^α) -(δ+n)  k/.{δ->0.1,n->0.01,s->0.15,α->0.25}},{t,0,30},{k,0,2.5},VectorPoints->Fine]

exactsol/.{α->0.25,δ->0.1,n->0.01,s->0.1}//Simplify

初期条件 a をいろいろ変えてグラフを描いてみる．

Plot[Table[exactsol/.{α->0.25,δ->0.1,n->0.01,s->0.1},{a,0.1,2,0.1}],{t,0,40},PlotRange->Full]

今度は，s を色々変えてグラフを描いてみる．

Plot[Table[exactsol/.{δ->0.1,n->0.01,α->0.25,a->0.1},{s,0.1,0.9,0.1}],{t,0,80},PlotRange->Full]