id='angtan(x,y)' algoritmo='sí' expresión='at' local='at' inicio='at=(abs(x)>0.00001)?atan(abs(y/x)):pi/2' hacer='at=((x<0)&(y>=0))?pi-at:at;at=((x<0)&(y<0))?pi+at:at;at=((x>=0)&(y<0))?2*pi-at:at;' tipo='función' 
id='Angle2D(ax,ay,bx,by)' algoritmo='sí' expresión='acos((ax*bx+ay*by)/(a*b))%pi' local='a,b' hacer='a=sqrt(ax^2+ay^2);b=sqrt(bx^2+by^2)' tipo='función' 
id='Angle3D(ax,ay,az,bx,by,bz)' algoritmo='sí' expresión='ang' local='a,b,ang' hacer='a=sqrt(ax^2+ay^2+az^2);b=sqrt(bx^2+by^2+bz^2);ang=acos((ax*bx+ay*by+az*bz)/(a*b))%pi;' tipo='función' 
id='Det2x2(a11,a12,a21,a22)' algoritmo='no' expresión='a11*a22-a12*a21' tipo='función' 
id='solve2x2(a1,a2,b1,b2,c1,c2)' algoritmo='sí' expresión='ok' local='det,ok' hacer='det=Det2x2(a1,a2,b1,b2);ok=(abs(det)>0.000001);tt=(ok)?Det2x2(c1,c2,b1,b2)/det:&squot;NaN&squot;;ss=(ok)?Det2x2(a1,a2,c1,c2)/det:&squot;NaN&squot;;' tipo='función' 
id='Det3x3(a11,a12,a13,a21,a22,a23,a31,a32,a33)' algoritmo='no' expresión='a11*(a22*a33-a32*a23)+a12*(a23*a31-a21*a33)+a13*(a21*a32-a31*a22)' tipo='función' 
id='Det4x4(x1,y1,z1,w1,x2,y2,z2,w2,x3,y3,z3,w3,x4,y4,z4,w4)' algoritmo='no' expresión='x4*Det3x3(y1,z1,w1,y2,z2,w2,y3,z3,w3)-y4*Det3x3(z1,w1,x1,z2,w2,x2,z3,w3,x3)+z4*Det3x3(w1,x1,y1,w2,x2,y2,w3,x3,y3)-w4*Det3x3(x1,y1,z1,x2,y2,z2,x3,y3,z3)' hacer=';' tipo='función' 
id='NewtonRaphson()' algoritmo='sí' expresión='(diff2<=eps)' local='eps,x,y,diff2,GF2,F0,Fx,Fy,n' inicio='eps=1/2^25;n=0' hacer='x=xi;y=yi;F0=_Eval_(F);x=xi+epsilon;Fx=(_Eval_(F)-F0)/epsilon;x=xi;y=yi+epsilon;Fy=(_Eval_(F)-F0)/epsilon;GF2=(Fx^2+Fy^2);x=xi-Fx*F0/GF2;y=yi-Fy*F0/GF2;diff2=(x-xi)^2+(y-yi)^2;xi=x;yi=y;n=n+1;' mientras='(diff2>eps)&(n<256)' tipo='función'