% test myT5
clf;

x = (-1:.001:1)';
y = myT5(x,4);
subplot(121); plot(x,y);
title('Formula 2');

x = (-1.05:.001:1.05)';
subplot(122); 
for flag = [1 3:6]
    y = myT5(x,flag);
    plot(x,y); 
    hold on
end
title('Formulas 1,3,4,5,6');
hold off