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;