请求：计算混沌系统李雅普诺夫指数的MATLAB程序

文献中的混沌系统为

参数及初始条件如下，里面的李雅普诺夫指数是怎么计算的？MATLAB程序是什么？