kill(all)$ load(ntrig)$ load(atrig1)$

van=facebook.com/LearningCalculus$

rad(theta):=(%pi/180)*theta;

f(x):=trigsimp((1+tan(rad(x))+1/cos(rad(x)))*(1+1/tan(rad(x))+1/sin(rad(x))));

'sum(''trigsimp(f(x)),x,1,50);

float(sum(f(k),k,1,50));