I need to plot the function f(k) which is defined as a square root of the time integral of the square of the Euclidean norm of the system state vector x(t) during the first 10 seconds after the start of movement. Input is matrix A, vector b and inital state x0 for a system of differential equations dx/dt = Ax + bu. u is called control signal and is of the form u = kx. Dimension of vector x is 2. If system is unstable f(k) should return -1.