clear all; close all n=100; k=[1:n]'; c=ones(n,1); b=ones(n-1,1); a=cos(k/5)+1; tdes=1; Tdes=cos(k/13); Tmax=2.5;