用matlab解带参数的如下的常微分方程：
function f=fun(t,x)
dx1=x(2);
dx2=-sin(kk*t)-0.2*x(2)-x(1);
f=[dx1;dx2];
end
参数kk=0:0.01:2变化，请问如何在程序中实现，求出对应于每一个kk下的x的最大值，我只会在操作中输入：
[t,x]=ode45(@fun,[0:0.1:300],[0 0]);
a=max(abs(x))
能不能在程序中自动求出所有kk下的x的最大值，然后绘出max(X)与kk的关系图，谢谢