Lyapunov exponent with respect to a parameter

I want to investigate the chaotic nature of Lorenz system when the parameter R varies. I am using the code which was previously posted on Mathworks. Now, for every R, I want to take the Lyapunov exponent at the end time. Then, I would like to store that R and Lyapunov exponent and plot them together. To get this kind of plot, I made some changes, but it gives some errors. Please help.

[Haneche]

Hi my brother, i try this three codes in my PC but an error appear , can you help me ?

[DGM]

You ran some "modified" code that was known to have some undescribed errors, and you got some undescribed errors.

FWIW, the unmodified files are here:

https://www.mathworks.com/matlabcentral/fileexchange/4628-calculation-lyapunov-exponents-for-ode

If you want to conduct the same modification that OP wanted, you'll have to figure out what OP wanted, because the changes that were made don't really make sense. R is not defined before it's used, so you'll have to figure out where and how it should be defined and possibly pass it between functions if and when necessary.