I'm trying to use Sagemath notebook environment without much prior knowledge of Sage. What I tried to do was define some variables (R, a, x, y, x, t...) and then some functions expressed in these variables. The expressions get quite large and complex this way as they are functions built of other functions. Last line of my code I'm trying to expand a large complex-valued symbolic expression using the variables of the functions (so not their values). I use expand() and collectterms(x,y,z,t) to get a result that is sorted by powers of x,y,z,t. When I used this with less complicated expressions it seemed to work OK giving me some reassurance that my syntax should be correct. But when the expressions get larger (see my code below) the outcome isn't correct anymore. I get a result without error messages, but the result is not correct. For example : the code below gives me a result without any imaginary numbers, also the real part is wrong. Does anyone know what is going on here? Is my syntax wrong? Or does Sage fail without error when expressions get large? And if so: are there ways to do this correctly? Thanks in advance!

` `

clos(R,a)=((1/(2*R))+R/2)*cos(a)+I*((1/(2*R))-R/2)*sin(a)
slin(R,a)=((1/(2*R))+R/2)

*sin(a)-I*((1/(2

*R))-R/2)*cos(a)

a_1(R,a)=(1/5)

*sqrt(5)*

a_2(R,a)=-(1/20)sqrt(5)-(1/4)+(1/2)

a_2(R,a)=-(1/20)

*clos(R,a)*

a_3(R,a)=-(1/20)sqrt(5)+(1/4)+(1/2)

a_3(R,a)=-(1/20)

*slin(R,a)*

a_4(R,a)=-(1/20)sqrt(5)-(1/4)-(1/2)

a_4(R,a)=-(1/20)

*clos(R,a)*

a_5(R,a)=-(1/20)sqrt(5)+(1/4)-(1/2)

a_5(R,a)=-(1/20)

*slin(R,a)*

b_1(R,a)=(1/5)sqrt(5)

b_1(R,a)=(1/5)

b_2(R,a)=-(1/20)

*sqrt(5)+(1/4)-(1/2)*slin(R,a)

b_3(R,a)=-(1/20)

*sqrt(5)-(1/4)+(1/2)*clos(R,a)

b_4(R,a)=-(1/20)

*sqrt(5)+(1/4)+(1/2)*slin(R,a)

b_5(R,a)=-(1/20)

*sqrt(5)-(1/4)-(1/2)*clos(R,a)

cc_1(R,a)=(1/5)

*sqrt(5)*

cc_2(R,a)=-(1/20)sqrt(5)-(1/4)-(1/2)

cc_2(R,a)=-(1/20)

*clos(R,a)*

cc_3(R,a)=-(1/20)sqrt(5)+(1/4)-(1/2)

cc_3(R,a)=-(1/20)

*slin(R,a)*

cc_4(R,a)=-(1/20)sqrt(5)-(1/4)+(1/2)

cc_4(R,a)=-(1/20)

*clos(R,a)*

cc_5(R,a)=-(1/20)sqrt(5)+(1/4)+(1/2)

cc_5(R,a)=-(1/20)

*slin(R,a)*

d_1(R,a)=(1/5)sqrt(5)

d_1(R,a)=(1/5)

d_2(R,a)=-(1/20)

*sqrt(5)+(1/4)+(1/2)*slin(R,a)

d_3(R,a)=-(1/20)

*sqrt(5)-(1/4)-(1/2)*clos(R,a)

d_4(R,a)=-(1/20)

*sqrt(5)+(1/4)-(1/2)*slin(R,a)

d_5(R,a)=-(1/20)

*sqrt(5)-(1/4)+(1/2)*clos(R,a)

r_1(R,a,x,y,z,t)=-c_1

*(1/5)+a_1(R,a)*x+ b_1(R,a)

*t+ cc_1(R,a)*z+ d_1(R,a)

*y*

r_2(R,a,x,y,z,t)=-c_1(1/5)+a_2(R,a)

r_2(R,a,x,y,z,t)=-c_1

*x+ b_2(R,a)*t+ cc_2(R,a)

*z+ d_2(R,a)*y

r_3(R,a,x,y,z,t)=-c_1

*(1/5)+a_3(R,a)*x+ b_3(R,a)

*t+ cc_3(R,a)*z+ d_3(R,a)

*y*

r_4(R,a,x,y,z,t)=-c_1(1/5)+a_4(R,a)

r_4(R,a,x,y,z,t)=-c_1

*x+ b_4(R,a)*t+ cc_4(R,a)

*z+ d_4(R,a)*y

r_5(R,a,x,y,z,t)=-c_1

*(1/5)+a_5(R,a)*x+ b_5(R,a)

*t+ cc_5(R,a)*z+ d_5(R,a)

*y*

uu(R,a,x,y,z,t)=r_1(R,a,x,y,z,t)r_2(R,a,x,y,z,t)

uu(R,a,x,y,z,t)=r_1(R,a,x,y,z,t)

*r_3(R,a,x,y,z,t)*r_4(R,a,x,y,z,t)+ r_1(R,a,x,y,z,t)

*r_2(R,a,x,y,z,t)*r_3(R,a,x,y,z,t)

*r_5(R,a,x,y,z,t)+ r_1(R,a,x,y,z,t)*r_2(R,a,x,y,z,t)

*r_4(R,a,x,y,z,t)*r_5(R,a,x,y,z,t) + r_1(R,a,x,y,z,t)

*r_3(R,a,x,y,z,t)*r_4(R,a,x,y,z,t)

*r_5(R,a,x,y,z,t) + r_2(R,a,x,y,z,t)*r_3(R,a,x,y,z,t)

*r_4(R,a,x,y,z,t)*r_5(R,a,x,y,z,t)

(uu(i,j,x,y,z,t)).expand().maxima_methods().collectterms(x,y,z,t)