BARZ-oldies
Digits:=30;
Speciaux
special:= proc(nn, c) option remember; special||nn(c) end;
X368, solution d'une \303\251quation du 5\303\250me degr\303\251
K018
# eqcub:= subs(zipq(pp, pX(6)), zipq(pu, pX(523)), k=0, cub1);
anticomplem de Kieppert hyperbola
#conicir(pX(523)): subs(zipq(px, complem(px)), %): eqcon:= collect(numer(%), px);
special368:= proc(_c_) local eqcub,eqcon,eqs0,eqs1,eli,solz;
eqcub := x*(b^2*z^2+y^2*c^2)*(b^2-c^2)+y*(z^2*a^2+c^2*x^2)*(-a^2+c^2)+z*(a^2*y^2+b^2*x^2)*(a^2-b^2);
eqcon := (b^2-c^2)*x^2+(-a^2+c^2)*y^2+(a^2-b^2)*z^2;
eqs0:= {eqcon, eqcub, x+y+z=1}: (cat)('en',_c_,"y_");
eqs1:= subs(%, eqs0);
eli:= eliminate(eqs1, {x,y}): solz:= fsolve(op(eli[2])/(3*z-1),{z});
subs(eli[1], solz, [x,y,z]):;
end:
# special(368,"c"); pX(368);
X369 et X3232, solution d'une equation du 3\303\250me degr\303\251
special369:= proc(_c_) local yffeq, yffsol;
yffeq:= subs(t=K, 2*t^3 - 3*(a + b + c)*t^2 + (a^2 + b^2 + c^2 + 8*b*c + 8*c*a + 8*a*b)*t - (c*b^2 + a*c^2 + b*a^2 + 5*b*c**2 + 5*c*a**2 + 5*a*b**2 + 9*a*b*c));
yffsol:= (K^2-(2*c + a)*K + (- a**2 + b**2 + 2*c**2 + 2*b*c + 3*c*a + 2*a*b));
(cat)('en',_c_,"y_");
subs(fsolve(subs(%, yffeq), {K}), %, [rot3](yffsol)):
%/add(k, k=%); end:
#special(369,"c");%-pX(369);
special3232:= proc(_c_); map(1/id, special(369,_c_));%/add(k, k=%); end:
#special(3232,"c"); %-pX(3232);
# latexx(yffeq); latexy(yffsol); latexz(pX(369)), ency(%);
X370
special370:= proc(_c_) local con1, eqs, eli; con1 :=
(x*y*(a^2-b^2+c^2)+x*z*(a^2+b^2-c^2)+2*z*y*a^2)*sqrt(3)-a*b*c*x*(y+z+2*x)/R;
(cat)('en',_c_,"y_");
eqs:= subs(%, {con1, rot(con1)}); eli:= [eliminate](%, {y,z}):
evalf(eli[1]); (subs)(%[1], [1,Re(y/x), Re(z/x)]); %/add(k,k=%); end:
# 370: special(%,"c"); %-pX(%%);
X1144
special1144:= proc(_c_) local eqs, fso,tmpbar;
eqs:= a**2/(a - L) + b**2/(b - L) + c**2/(c - L) - 2*S/L;
-(numer@subs)(ency_,%): collect(%, L); fsolve(%, L); fso:= min(%);
tmpbar:= [a**2/(a - L), b**2/(b - L), c**2/(c - L)];
(cat)('en',_c_,"y_");
subs(L=fso, %, tmpbar): %/add(k,k=%); end:
# 1144: special(%,"c"); %-pX(%%);
# con1:= conicir(pX(649));
con1 := (b-c)*a^2/x+(c-a)*b^2/y+c^2*(-b+a)/z;
(simplify@subs)(zipq([x,y,z],special(1144,"c")), ency_, con1),
(simplify@subs)(zipq([x,y,z],special(1144,"z")), enzy_, con1);
X2061 trop long
special2061:= proc(_c_) local orion, pX_40;
orion:= `@`(proc (p, q, r) options operator, arrow; [-(-b^2*p^2*r^2-p^2*r*b^2*q+p^2*r*a^2*q-q^2*c^2*p^2-p^2*r*c^2*q+q^2*a^2*r^2)*p, (-b^2*p^2*r^2+q^2*c^2*p^2-p*r*b^2*q^2+p*r*c^2*q^2+p*r*a^2*q^2+q^2*a^2*r^2)*q, (b^2*p^2*r^2-q^2*c^2*p^2+p*r^2*b^2*q+p*r^2*a^2*q-p*r^2*c^2*q+q^2*a^2*r^2)*r] end proc, OP):
pX_40:= [a*(a^3+b*a^2-a*b^2-b^3+c*a^2-2*a*b*c+b^2*c-a*c^2+b*c^2-c^3), b*(b^3+b^2*c-b*c^2-c^3+a*b^2-2*a*b*c+a*c^2-b*a^2+c*a^2-a^3), c*(c^3+a*c^2-c*a^2-a^3+b*c^2-2*a*b*c+b*a^2-b^2*c+a*b^2-b^3)];
(cat)('en',_c_,"y_");
subs(%, orion(pX_40)); evalf(%/add(k,k=%));
end:
# 2061: special(%,"c"); %-pX(%%);
Calcul de ff (point sur infty line)
Infinity--- erreur, due \303\240 des tailles absolues tr\303\250s petites
ckoi:= [ 30, 511, 512, 513, 514, 515, 516, 517, 518, 519, 520, 521, 522, 523, 524, 525, 526, 527, 528, 529, 530, 531, 532, 533, 534, 535, 536, 537, 538, 539, 540, 541, 542, 543, 544, 545, 674, 680, 688, 690, 696, 698, 700, 702, 704, 706, 708, 710, 712, 714, 716, 718, 720, 722, 724, 726, 730, 732, 734, 736, 740, 742, 744, 746, 752, 754, 758, 760, 766, 768, 772, 776, 778, 780, 782, 784, 786, 788, 790, 792, 794, 796, 802, 804, 806, 808, 812, 814, 816, 818, 824, 826, 829, 830, 832, 834, 838, 888, 891, 900, 912, 916, 918, 924, 926, 928, 952, 971, 1105, 1154, 1288, 1499, 1503, 1510, 1912, 1938, 1946, 2121, 2131, 2133, 2385, 2386, 2387, 2388, 2389, 2390, 2391, 2392, 2393, 2574, 2575, 2697, 2706, 2771, 2772, 2773, 2774, 2775, 2776, 2777, 2778, 2779, 2780, 2781, 2782, 2783, 2784, 2785, 2786, 2787, 2788, 2789, 2790, 2791, 2792, 2793, 2794, 2795, 2796, 2797, 2798, 2799, 2800, 2801, 2802, 2803, 2804, 2805, 2806, 2807, 2808, 2809, 2810, 2811, 2812, 2813, 2814, 2815, 2816, 2817, 2818, 2819, 2820, 2821, 2822, 2823, 2824, 2825, 2826, 2827, 2828, 2829, 2830, 2831, 2832, 2833, 2834, 2835, 2836, 2837, 2838, 2839, 2840, 2841, 2842, 2843, 2844, 2845, 2846, 2847, 2848, 2849, 2850, 2851, 2852, 2853, 2854, 2867, 2869, 2870, 2871, 2872, 2873, 2874, 2875, 2876, 2877, 2878, 2879, 2880, 2881, 2882, 3221, 3307, 3308, 3309, 3413, 3414, 3479]: nops(%);
Le codage \303\251tait mauvais....
convert(ckoi,set) minus {seq}(`if`(ff[j]=0, NULL,j), j=1..jmax);
Calcul de qff --- pr\303\251c\303\250de le calcul des xyz -- tient compte de la taille absolue des nombres
qff:= table([seq](j=0, j=1..jmax)):
for j from 1 to jmax do
if fdat[j]="special" then next fi; [rot3](parse(fdat[j])):
tmpc, tmpz:= evalf(subs(ency_, %)), evalf(subs(enzy_,%));
abs(add(k,k=tmpc))/add(abs(k),k=tmpc)+abs(add(k,k=tmpz))/add(abs(k),k=tmpz);
if %< Float(1,-7) then qff[j]:=1 fi;
od:
# seq(`if`(qff[j]=ff[j], NULL,j), j=1..jmax);
[seq](`if`(qff[j]=1,j,NULL),j=1..jmax): nops(%);
ff:= eval(qff):
Calcul des xyz, zxyz
_a_,_b_,_c_:= (op@subs)(ency_, [a,b,c]);
sqrt((a+b+c)*(b+c-a)*(a+c-b)*(a+b-c)): subs(a=_a_,b=_b_,c=_c_,%/_a_/2):
fac:= evalf(%);
xyz, zxyz, sk:= table(), table(), table():
for j to jmax do
if fdat[j]="special" then xyz[j]:= special(j,"c"); zxyz[j]:= special(j,"z");
sk[j]:= evalf(xyz[j][1]*fac);
next fi;
tmq:= [rot3](parse(fdat[j])):
if ff[j]=1 then
tmp:= evalf(subs(ency_, tmq)); xyz[j]:= evalf(tmp*add(1/k, k=tmp));
sk[j]:= evalf(tmp[1]/_a_*(_a_/tmp[1]+_b_/tmp[2]+_c_/tmp[3]));
tmp:= evalf(subs(enzy_, tmq)); zxyz[j]:= evalf(tmp*add(1/k, k=tmp));
else
tmp:= evalf(subs(ency_, tmq)); xyz[j]:= evalf(tmp/add(k, k=tmp));
sk[j]:= evalf(xyz[j][1]*fac);
tmp:= evalf(subs(enzy_, tmq)); zxyz[j]:= evalf(tmp/add(k, k=tmp));
fi;
# if j mod 200 = 0 then lprint("--",j) fi
if j=3519 then print(searchkey,j,sk[j]) fi;
od:
[seq](`if`( abs(sk[j]-enc_dat[j])<Float(1,-7), NULL,j), j=1..3514);:
[seq]([j, abs(sk[j]-enc_dat[j])], j=%);
enc_dat:= [seq(sk[j], j=1..jmax)]: siz_enc:= nops(%): enc_sort:= sort(enc_dat):
On teste l'int\303\251gration des nouveaux points (et on ne passe pas par pX, qui plante ASSERT .
# pXX:= proc (n) ([rot3])(parse(fdat[n])); if ency(%) <> n then print(n) fi end proc:
# for j from 3515 to jmax do pXX(j); od:
Calcul de gg,tt,cc,aa
forget(pX);
qgg,qtt,qcc,qac:=
table([seq](j=0, j=1..jmax)), table([seq](j=0, j=1..jmax)),
table([seq](j=0, j=1..jmax)), table([seq](j=0, j=1..jmax)):
for j to jmax do tmq:= xyz[j];
tmp:= ency( isogon(tmq)); if tmp<> `?` then qgg[j]:= tmp; fi;
tmp:= ency( isotom(tmq)); if tmp<> `?` then qtt[j]:= tmp; fi;
tmp:= ency( complem(tmq)); if tmp<> `?` then qcc[j]:= tmp; fi;
tmp:= ency(anticomplem(tmq)); if tmp<> `?` then qac[j]:= tmp; fi;
if j mod 1000 = 0 then lprint("--",j) fi
od: # j:= 'j':
jjmax:=3514:
# seq(`if`(qff[j]=ff[j], NULL,j), j=1..jjmax);
seq(`if`(qgg[j]=gg[j], NULL,j), j=1..jjmax);
seq(`if`(qtt[j]=tt[j], NULL,j), j=1..jjmax);
seq(`if`(qcc[j]=cc[j], NULL,j), j=1..jjmax);
seq(`if`(qac[j]=ac[j], NULL,j), j=1..jjmax);
Mauvaise identification, s'il y en a
seq(`if`(qgg[j]=0 or qgg[qgg[j]]=j, NULL,j), j=1..jmax);
seq(`if`(qtt[j]=0 or qtt[qtt[j]]=j, NULL,j), j=1..jmax);
seq(`if`(qcc[j]=0 or qac[qcc[j]]=j, NULL,j), j=1..jmax);
seq(`if`(qac[j]=0 or qcc[qac[j]]=j, NULL,j), j=1..jmax);
xres_ini(30, [` j`, ` ETC`, ` ok`]); for j to jjmax do if ac[j]<>0 and ac[j] <> qac[j] then xres_run([j, ac[j],qac[j]]) fi; od: xres_cut(); Dimension(Matrix(%));
# latexx(0), latexz(%%);
Il faut envoyer cela dans un r\303\251pertoire lisible et navigable
if false then
try close(fd); catch : end try;
fd := open("/home/douillet/public_html/etc/alt_igtca.csv", WRITE):
for j to jmax do
fprintf(fd, """%d"";""%d"";""%d"";""%d"";""%d"";""%d""\134n",
j, ff[j], qgg[j], qtt[j], qcc[j], qac[j]);
od: close(fd); j:='j':
fi:
# j, ff[j], gg[j], tt[j], cc[j], ac[j];
save
if false then
try close(fd); catch : end try;
fd := open("/home/douillet/public_html/etc/alt_igtca.csv", WRITE):
for j to jmax do
map(U-> RegSubs("0*0$"="", sprintf("%6.25f",U)),
evalf([op(xyz[j]),op(zxyz[j]), sk[j]],20));
fprintf(fd,
"""%d"";""%d"";""%d"";""%d"";""%d"";""%d"";""%s"";""%s"";""%s"";""%s"";""%s"";""%s"";""%s"";""%s"";""%s""\134n",
j, ff[j], qgg[j], qtt[j], qcc[j], qac[j], fdat[j], op(%), "nono");
od: close(fd); j:='j':
fi:
# nono est gera
1;
R\303\251sistants
3513
Vector(zalo);
"(a-b+c)*(a+b-c)+4*r*(r^2+4*r*R)^1/2"; substax(%);
a*parse(%); subs(kitcircle,%); map(factor, %);
simplify(%) assuming a>0,b>0,c>0,R>0:
collect(%, R, factor); (factor@subs)(valR22, %):
elimifac(%): tmp:= simplify(%) assuming a>0,b>0,c>0,R>0;
(factor@isolate)(defR2, b-c+a); subs(%, tmp);
simplify(%) assuming a>0,b>0,c>0,R>0:
tmqx:= xcollect(%*rotp(sqrt(b))/(a+b+c), sqrt(-b^2+2*a*b+2*b*c-a^2-c^2+2*a*c), factor);
ency(tmqx);
3472
fac:= u/a: (expand@cevadiv)(pp,barydiv(pX(1), pu));
tmp1:= expand(%[1]*fac/rotp(p))/fac;
sub1:= zipq(pp, map(U->expand(U), evalmm(simplify(pX(20), a2toAA)/AA/BB/CC)));
barydiv(pX(1), pX(3345)): expand(evalmm(simplify(%, a2toAA)/AA/BB/CC/2));
sub2:= zipq(pu, %);
subs(sub1, sub2, tmp1); tmqx:= subs(AAtoa2, %); length(convert(%, string));
tmqx:= a^2/tmqx; ency(%);
simplify(pX(3345), a2toAA);
simplify(pX(20), a2toAA);
simplify(pX(3345), a2toAA) :
cevadiv(%%,%);
3418
www; Vector(zalo);
14,16; vertexthird(pX(%[1]),pX(%[2])): (elimifac@factor)(%[1]);
subs(valR22, collect(%, R));;
tmqx:= xcollect(%, [R, sqrt(3)], bctax@factor);
(elimifac@factor@subs)(apbpc,%);
tmqx:= xcollect(%, [R, sqrt(3)], bctax@factor);
tmqx:= a^2/collect(denom(tmqx)/sqrt(3),R, factor); ency(%);
3353
cutt();
tmqx;
# forget(pX):
www; tmp:= wedge(wedge(pX(3), pX(2131)), wedge(pX(64), pX(3346))):
tmp1:= elimifac(factor(tmp[1])); %/tmqx;
a2toAA;
(tmp1): simplify(%, a2toAA); tmp2:= collect(%, AA, factor);
subs(AAtoa2, tmp2);; length(convert(%,string));
tmp10:= (-16*4*(b^2-c^2)^2*a^8*(-a^2+b^2+c^2)^6+
16*(a^4+2*b^2*c^2-b^4-c^4)*a^2*((a^2-b^2+c^2)^2+(a^2+b^2-c^2)^2)*(b^2-c^2)^2*(-a^2+b^2+c^2)^5
+16*(a^4+2*b^2*c^2-b^4-c^4)^2*(b^2-c^2)^2*a^4*(-a^2+b^2+c^2)^4
-8*(a^4+2*b^2*c^2-b^4-c^4)^3*a^2*((a^2-b^2+c^2)^2-(a^4+2*b^2*c^2-b^4-c^4)+(a^2+b^2-c^2)^2)*(-a^2+b^2+c^2)^3
+4*(a^4+2*b^2*c^2-b^4-c^4)^4*a^4*(-a^2+b^2+c^2)^2
+4*(a^4+2*b^2*c^2-b^4-c^4)^5*a^2*(-a^2+b^2+c^2)
-(a^4+2*b^2*c^2-b^4-c^4)^6)/
(((a^2-b^2+c^2)^2+(a^2+b^2-c^2)^2)*(-a^2+b^2+c^2)^2
-(a^4+2*b^2*c^2-b^4-c^4)^2);
ency(%), length(convert(%,string));
"(16*(b^2-c^2)^2*a^8*(-a^2+b^2+c^2)^6
-4*(a^4+2*b^2*c^2-b^4-c^4)*a^2*((a^2-b^2+c^2)^2+(a^2+b^2-c^2)^2)*(b^2-c^2)^2*(-a^2+b^2+c^2)^5
-4*(a^4+2*b^2*c^2-b^4-c^4)^2*(b^2-c^2)^2*a^4*(-a^2+b^2+c^2)^4
+2*(a^4+2*b^2*c^2-b^4-c^4)^3*a^2*((a^2-b^2+c^2)^2
-(a^4+2*b^2*c^2-b^4-c^4)+(a^2+b^2-c^2)^2)*(-a^2+b^2+c^2)^3-(a^4+2*b^2*c^2-b^4-c^4)^4*a^4*(-a^2+b^2+c^2)^2-(a^4+2*b^2*c^2-b^4-c^4)^5*a^2*(-a^2+b^2+c^2)+(a^4+2*b^2*c^2-b^4-c^4)^6/4)/(((a^2-b^2+c^2)^2+(a^2+b^2-c^2)^2)*(-a^2+b^2+c^2)^2-(a^4+2*b^2*c^2-b^4-c^4)^2)":
SubstituteAll(%, "-a^2+b^2+c^2","b^2+c^2-a^2"):
length(%); tmp11:= parse(%%): ency(%);
seq(factor(op([1,j],tmp10)/op([1,j],tmp11)), j=1..7);
tmqx:= tmp11;
cutt23(): tmqx:= elimifac(%); ency(%);
fintax(tmqx); length(%);
bary2norm(pX(2096))+bary2norm(pX(3421)): tmqx:= elimifac(factor(%[1]));; ency(%);
"(5*a^12-10*(b^2+c^2)*a^10+(-9*b^4+34*b^2*c^2-9*c^4)*a^8+36*(b^2+c^2)*(b^2-c^2)^2*a^6-(29*b^4+54*b^2*c^2+29*c^4)*(b^2-c^2)^2*a^4+2*(3*b^2+c^2)*(b^2+3*c^2)*(b^2+c^2)*(b^2-c^2)^2*a^2+(b^2-c^2)^6)/(-a^2+b^2+c^2)/(a^16-8*(b^2+c^2)*a^14+(28*b^4-40*b^2*c^2+28*c^4)*a^12-56*(b^2+c^2)*(b^2-c^2)^2*a^10+2*(35*b^4+114*b^2*c^2+35*c^4)*(b^2-c^2)^2*a^8-8*(b^2+c^2)*(7*b^4+18*b^2*c^2+7*c^4)*(b^2-c^2)^2*a^6+4*(b^4+7*c^4)*(7*b^4+c^4)*(b^2-c^2)^2*a^4-8*(b^2+c^2)*(b^2-c^2)^6*a^2+(b^4+14*b^2*c^2+c^4)*(b^2-c^2)^6)": ency(parse(%)), length(%); tmqx:= parse(%%);
(factor@crossdiv)(pX(6), pX(3350)): tmqx:= elimifac(%[1]);;
reflection(pX(3),pX(2131)): tmqx:= (elimifac@factor)(%[1]);;
tmp:= (cevadiv)(pX(20),pX(84))[1];
op(1, tmp): tmqx:= subs(%= map(U-> U/rotp(-3*c^4+2*c^2*(a^2+b^2)+(a^2-b^2)^2), %), tmp);;
# tmqy:= tmqx:
tmqx:= a^2/tmqy;
3129
tmqy:=tmqx;
xcollect(tmqy/rotp(1)/sqrt(3),[R,sqrt(3)],factor):
tmqx:= collect(-subs(apbpc,%),R, factor);
tmqx:= collect(subs(valR22,tmqy)/sqrt(3), R, factor);
(factor@cevadiv)(pX(12),pX(10))[1]; tmqx:= elimifac(%); ency(%);
2902
pX(62); isol:= isolate(valR, sqrt(a+b+c)); subs(map(1/id, %), pX(61));
tmqx:= %[1];
www; (elimifac@factor)(excentral_isogon(pX(62))[1]):
tmp:= subs(isol, collect(expand(%), R, factor));
collect(subs(valR22, tmp), R, factor):;
elimifac(factor(%)):
tmqx1:= xcollect(%, [R, sqrt], factor);
#tmqx:= xcollect(%/a/b/c, [R, sqrt(3)], (U-> map(abctax,U))@factor);
(elimifac@factor@subs)(apbpc, tmqx1): tmqx:= xcollect(%, [R,sqrt(3)],(U-> map(abctax,U))@factor);
2598
Vector(zalo);
zalo[3][4..-1]; : tmp3:= substax(%);
zalo[2]; : tmp2:= substax(%);
zalo[1];"k(a,b,c)*c*a/(c+a)+h(a,b,c)*b*a/(b+a)" : tmp1:= substax(%);
(subs)(parse(tmp2), parse(tmp3), parse(tmp1))*a;
tmpq:= xcollect(expand(%), [rot3](sec(A)), distributed);
sec(A)+sec(B); %= (factor@subs)(kashi, convert(%, sincos));
subs(rot3(%), eval({rot3}(-%)), tmpq);
tmqx:= map(U->factor(U*R^0*rotp((a+b)*(a+b-c)/(a^2+b^2-c^2))*4), %);ency(%);
tmqx:= (elimifac@factor)(tmpq); ency(%);
(factor@subs)(kashi, expand(convert(tmpq, sincos))); tmqx:= elimifac(%);
sec(A-B): %=subs(kashi, map(1/id, defR2), expand(%));
rot3(factor(%)): subs(%, tmpq): tmqx:= map(U->U/4/rotp(a^3), %); ency(%);
tmpq:= (elimifac@factor@subs)(kashi, expand(%));
collect(tmpq, R, factor);
(elimifac@factor@subs)(defR2, map(id^2, defR2), map(id^3, defR2), tmpq); tmqx:=%:
2574
www;tmqx:= isoconj(pX(6), pX(1823) )[1];ency(%);
pX(1823)[1]; tmqx:= 1/xcollect(%/a, WW3, abctax); ency(%);
evalmm(pX(1382)*rotp(sqrt(a)))[1]; a^3/xcollect(%, WW4, factor); convert(%, string);
www; # tmqy:= tmqx; length(convert,%, string);
tmqx:= a^2/cos(A)/pX(1113)[1]; ency(%), length(convert(%, string));
tmqx:= a/cos(A)/pX(1113)[1]; ency(%), length(convert(%, string));
tmqx:= 1/a^2/cos(A)/pX(1113)[1]; ency(%), length(convert(%, string));
pX(4);
subs(kashi, 2*tmqx); rationalize(%); elimifac(%);
tmqx:= xcollect(%, WW3, factor);ency(%), length(convert(%, string));
tmqx:= 1/a*pX(1114)[1]; ency(%); [rot3](%%): tmqx:= isotom(%)[1]; ency(%);
# zipd([rot3](tmqy),[rot3](tmqx)): subs(ency_,%);evalmm(expand(%)/%[1]); factor(%);
cutt();
2534
kitcircle;
WW1;
WW2;
WW3:= sqrt(a^6-a^4*b^2-a^4*c^2-a^2*b^4+3*a^2*b^2*c^2-a^2*c^4+b^6-b^4*c^2-b^2*c^4+c^6);
www;
SubstituteAll(zalo[1],",JasatX(1113).",""); "(1-J/e)*cos(A)+4*cos(B)*cos(C)";
substax(%); ; tmp1:= parse(%); #tmp2:= x=x:
#substax(zalo[2]); "d=(3-2*cos(A)-2*cos(B)-2*cos(C))^(1/2)": tmp2:= parse(%);
tmp3:= subs(tmp2, valJ, kitcircle, kashi, expand(tmp1)*a); ency(%);
tmqx:= xcollect(2*tmp3*WW2, WW2, factor);
tmqx:= simplify(%) assuming a>0,b>0,c>0:
tmqx:= xcollect(%*rotp(sqrt(1)*a^2), WW2, factor); ency(%);
(factor@subs)(defR2, map(1/id, defR2), tmp3):
simplify(%) assuming a+b+c>0;
tmqx:= xcollect(elimifac(%), [WW2,R], factor); ency(%);
tmp1:= subs(kitcircle, s=(a+b+c)/2, kashi, expand(%) ); ency(%);
WW1:= sqrt(a^2*b^2+a^2*c^2+b^2*c^2);
WW2:= sqrt(a^4-a^2*b^2-a^2*c^2+b^4-b^2*c^2+c^4);
subs(map(1/id, defR2), map(1/id^2, defR2), defR2, tmp1):
tmqx:= xcollect(%*WW1*2*a*b*c, [WW2], normal); fintax(tmqx); parse(%); ency(%);
(tmp1): collect(subs(valom, kashi, expand(%)), R);
subs(map(1/id, defR2), defR2, %): factor(%): elimifac(%); subs(defR2,%);
tmqx:= xcollect(%, sqrt(b^4-b^2*c^2+c^4-a^2*b^2-a^2*c^2+a^4), factor);
ency(%);
fintax(tmqx); parse(%);
defR2; isolate(defR, S);
valom; # latexx(%);
(subs)(valom, cos(omega)^2+sin(omega)^2): collect(%,R);
(factor@simplify)(%, {defR2});
r=S/((a+b+c)/2); (factor@subs)(isolate(defR, S), %); valr;
2454
zalo;
tmp1:= parse(substax(zalo[3]));
substax(zalo[2]);
"K=(1/2)*(a^(8)+b^(8)+c^(8)-S26+a^(2)*b^(2)*c^(2)*(a^(2)+b^(2)+c^(2)))^(1/2)";
tmp2:= parse(%);
tmp2;
solve( defR, {S});;
valK:= subs(S=areaofABC, op(%)),subs(S=areaABC, op(%)), subs(tmp1, tmp2);
substax(zalo[1]); "b*c*(3*a^(2)*b*c*cos(B)*cos(C)-4*(S)^(2)-2*K)";
subs(solve( defR, {S}), K=-K, tmp2, tmp1, kashi, a*parse(%)); ency(%);
tmqx:= ((3*(c^2+a^2-b^2))*(a^2+b^2-c^2)-(1)*a^2*b^2*c^2/R^2+4*sqrt(a^8+b^8+c^8-a^2*(b^6+c^6)-b^2*(c^6+a^6)-c^2*(a^6+b^6)+a^2*b^2*c^2*(a^2+b^2+c^2))); ency(%);
fintax(tmqx); parse(%);
2146
saver[www-1]; [rot3](parse(%)):
tmp1:= isoconj(pX(2), %); indets(%);
subs(a=A^2,b=B^2,c=C^2, tmp1[1]): simplify(%, symbolic);
factor(%); elimifac(%);
subs(seq(C^(2*k)=c^k, k=[5,4,3,2,1]),
seq(B^(2*k)=b^k, k=[5,4,3,2,1]),
seq(A^(2*k)=a^k, k=[5,4,3,2,1]),
seq(C^(2*k+1)=sqrt(CC)*c^k, k=[5,4,3,2,1,0]),
seq(A^(2*k+1)=sqrt(AA)*a^k, k=[5,4,3,2,1,0]),
seq(B^(2*k+1)=sqrt(BB)*b^k, k=[5,4,3,2,1,0]), %);
xcollect(%, [sqrt(BB), sqrt(CC)], factor);
tmqx:= subs(AA=a, BB=b,CC=c,%); ency(%);
subs(A=sqrt(a),B=sqrt(b),C=sqrt(c), %): simplify(%, symbolic);
xcollect(%, [b^(3/2),c^(3/2), a]);
X(2146) = X(238)-CEVA CONJUGATE OF X(365)
tmp1:= cevadiv(pX(238),pX(365)); ency(%);
tmp1[1]: expand(%/rotp(a)); subs(seq(a^((2*k+1)/2)=sqrt(A)*a^k, k=0..4),
seq(b^((2*k+1)/2)=sqrt(B)*b^k, k=0..4), seq(c^((2*k+1)/2)=sqrt(C)*c^k, k=0..4), %):
xcollect(%, [sqrt(A), sqrt(B), sqrt(C)], factor, distributed);
-tmp1[1]/rotp((a^2-b*c)*a); xcollect(%,[a^(3/2), b^(3/2)], factor, distributed);
tmqx:= xcollect(%,[a^(3/2), b^(1/2)], factor);
2058
valR22:= defR2, map(1/id, defR2), map(id^2, defR2);
collect(orion(pX(69))[1], R, factor): elimifac(%):
tmqx:= map(abctax,%); length(convert(%, string)), ency(%);
:
collect(subs(valR22, %)/R, R, factor):
tmq1:= collect(subs(valR22, %)*R, R, factor);
(elimifac@factor)(tmq1): tmq2:= xcollect(%, [R, sqrt(3)], abctax);
coeff(tmq2,R,0): tmqx:= collect(%,a,factor)+select(has,tmq2,R);
xcollect(tmqx, R, factor);length(convert(%, string));
tmqx:= tmq2;;
length(convert(%, string));
ency(tmqx);
2038
WW:= sqrt(-a^2*c^2+c^4-b^2*c^2-a^2*b^2+a^4+b^4);
VV:= sqrt(a^2*b^2+a^2*c^2+b^2*c^2);
www; zalocom*a; (factor@subs)(kitcircle, %):
xcollect(subs(valom, kashi,expand(%)), [WW, VV, R], factor);
subs(defR2, map(1/id, defR2), map(1/id^2, defR2), %);
-elimifac(factor(%));
:
subs(defR2, map(1/id, defR2), %): subs(defR2, map(1/id, defR2),collect(%*R,R));
elimifac(factor(%)); subs(defR2, map(1/id, defR2),collect(%*R, R, factor));
elimifac(factor(%)): tmqx:= collect(%,R, factor); ency(%);
:
xcollect(-%*2*rotp(a), [WW, VV, R], factor);
tmqx:= xcollect(numer(%), [WW], factor); ency(%);
2037
zalo[2];"rr=(r^2+s^2)/(4*r)"; val_rr:= parse(%);
WW:= sqrt(-a^2*c^2+c^4-b^2*c^2-a^2*b^2+a^4+b^4);
zalo[1]; "e*a*cot(A)-(2*rr-e*s*csc(omega))*cos(A+omega)";
parse(%); tmp1:= subs(val_rr, kitcircle, valom, kashi, convert(a*expand(%),sincos));
tmp2:= xcollect(tmp1*4*R*rotp(a^2)*sqrt(a^2*b^2+a^2*c^2+b^2*c^2),[WW, R], factor);
subs(defR2, map(1/id, defR2), tmp2):
tmqx:= xcollect(%*rotp((a+b-c)/a^2)*(a+b+c)/2, WW, factor);
ency(tmqx);
2019
zalocom*a; subs(kitcircle, %): tmp1:= collect(subs(valom, kashi, expand(%)),R, factor);
collect(tmp1*R^2, R, factor): collect(subs(defR2,%)*R, R, factor):
collect(subs(defR2,%), R, factor); tmp2:= elimifac(factor(%));
tmqx:= collect(tmp2, R, factor);
ency(%);
2015
cutt();
tmqx: tmqx:= xcollect(%*R*rotp(a^5)/(a^2+b^2+c^2), [sqrt(a^6-a^4*b^2-a^2*b^4+b^6-a^4*c^2+3*b^2*a^2*c^2-b^4*c^2-c^4*a^2-b^2*c^4+c^6), sqrt(a^2*b^2+a^2*c^2+b^2*c^2)], factor);
tmqx:= xcollect(tmqx*rotp(a^2), [sqrt(a^6-a^4*b^2-a^2*b^4+b^6-a^4*c^2+3*b^2*a^2*c^2-b^4*c^2-c^4*a^2-b^2*c^4+c^6)], factor); ency(%);
subs(apbpc, %)*(a^2+b^2+c^2)/rotp(a^2):
tmqx:= xcollect(%, [sqrt(a^6-a^4*b^2-a^2*b^4+b^6-a^4*c^2+3*b^2*a^2*c^2-b^4*c^2-c^4*a^2-b^2*c^4+c^6)], factor); ency(%);
tmqx; subs(isolate(valR,sqrt(b+a+c)), %)*R:
collect(%, R, factor); subs(apbpc,%); tmqx:= collect(%/a/b/c, R, factor); ency(%);
1683
wedge(wedge(pX(3),pX(6)), wedge(pX(10), pX(1676))):
tmp1:= subs(valom, kashi, expand(%))[1]:
# ency(%);
WW:= op(indets(tmp1, sqrt));
xcollect(tmp1*R^3, [R,WW], factor):;
xcollect(subs(defR2, %), WW, factor);
elimifac(factor(%)): tmqx:= xcollect(%, 1/WW, factor); fintax(%); parse(%);
1662
zalo[1]; kitcircle;
substax(zalo[1]); parse(%)*a;
tmp1:= subs(kitcircle, s=(a+b+c)/2, kashi, expand(%) ); ency(%);
WW1:= sqrt(a^2*b^2+a^2*c^2+b^2*c^2);
WW2:= sqrt(a^4-a^2*b^2-a^2*c^2+b^4-b^2*c^2+c^4);
subs(map(1/id, defR2), map(1/id^2, defR2), defR2, tmp1):
tmqx:= xcollect(%*WW1*2*a*b*c, [WW2], normal); fintax(tmqx); parse(%); ency(%);
(tmp1): collect(subs(valom, kashi, expand(%)), R);
subs(map(1/id, defR2), defR2, %): factor(%): elimifac(%); subs(defR2,%);
tmqx:= xcollect(%, sqrt(b^4-b^2*c^2+c^4-a^2*b^2-a^2*c^2+a^4), factor);
ency(%);
fintax(tmqx); parse(%);
defR2; isolate(defR, S);
valom; # latexx(%);
(subs)(valom, cos(omega)^2+sin(omega)^2): collect(%,R);
(factor@simplify)(%, {defR2});
r=S/((a+b+c)/2); (factor@subs)(isolate(defR, S), %); valr;
1522
orthojoin(pX(13)): tmqx:= (factor@subs)(defR2, %[1]); ency(%);
elimifac(tmqx): map(xcollect,%,[R, sqrt(3)]);
a*((-18*a^5*c^4-3*a*c^8+12*a^3*b^6-18*a^5*b^4-3*a*b^8+12*a^7*b^2+12*a^7*c^2+12*a^3*c^6-3*a^9+12*a*b^2*c^6-18*a*b^4*c^4+12*a*b^6*c^2-12*a^3*b^2*c^4-12*a^5*b^2*c^2-12*a^3*b^4*c^2)*R/sqrt(3)+(2*c^9*b+2*b^9*c+12*b^5*c^5-8*b^3*c^7-8*b^7*c^3-7*a^4*b^5*c-9*a^6*b^3*c+5*b^7*a^2*c-5*a^2*b^3*c^5-7*a^4*c^5*b+9*a^8*b*c+14*a^4*b^3*c^3+5*c^7*a^2*b-9*a^6*b*c^3-5*b^5*c^3*a^2))*(a^6*b^2+a^6*c^2-3*a^4*b^4+2*a^4*b^2*c^2-3*a^4*c^4+3*a^2*b^6-2*a^2*b^4*c^2-2*a^2*b^2*c^4+3*a^2*c^6-b^8-b^6*c^2+4*b^4*c^4-b^2*c^6-c^8): tmqx:= map(xcollect,%,[R, sqrt(3)], abctax);
tmpp:= fintax(tmqx); parse(%); ency(%);
www:=62; pX(62); isolate(valR,sqrt(a+b+c));
subs(map(1/id, %),%%); (elimifac@factor)(%[1]): collect(%, R, factor);
tmpp:= fintax(%); parse(%); ency(%);
www;orthojoin(pX(589));
tmqy:= (elimifac@factor@subs)(defR2, %[1]); ency(%);
tmqy: length(convert(%, string));
tmqx:= map(collect,factor(tmqy),R, abctax);length(convert(%, string));
subs(defR2, collect(tmqy, R)): collect(factor(%), R);;length(convert(%, string));
elimifac(factor(subs(defR2, collect(expand(tmqy), R)))): collect(%, R, abctax);
tmpp:= convert(%, string); parse(%); ency(%);length(convert(tmpp, string));
map(xcollect,tmqy,[R, sqrt(3)]):
collect(select(has,%,R)/sqrt(1), R, normal)*remove(has,%,R);
tmqx:= map(xcollect,%,[R, sqrt(3)], abctax);
tmpp:= fintax(%); parse(%); ency(%);
tmqy:=cutt();
map(xcollect,tmqy,[R, sqrt(3)]):
collect(select(has,%,R)/sqrt(3), R, normal)*remove(has,%,R);
tmqx:= map(xcollect,%,[R, sqrt(3)], abctax);
tmpp:= fintax(%); parse(%); ency(%);
1113, 1314
seq(``||u||2=u^2,u=[a,b,c]);
valJ:= J= subs( J = 1/(a*b*c)*(S(6) - S(2,4) + 3*a^2*b^2*c^2)^(1/2), S(6) = a^6 + b^6 + c^6,
S(2,4) = a^2*b4 + a2*c4 + b2*c4 + b2*a4 + c2*a4 + c2*b4,
seq(``||u||4=u^4,u=[a,b,c]),seq(``||u||2=u^2,u=[a,b,c]), J);
substax("(1 - J)* cos (A) - 2* cos (B) *cos (C)"); parse(%);
subs(valJ,%)*sin(A); subs(kashi, %); WW:= op(indets(%, sqrt));
xcollect(%%,WW, factor); elimifac(factor(%)); xcollect(-%,WW, factor@expand);
ency(%);
fac:=a^2*b^2*c^2/R^2/4*(a+b+c); sqrtfac:= a*b*c/R/2*sqrt(a+b+c);
uuu:= (subs)(f(a,b,c) = (d^2 + (4*r - R)*R + sqrt(Q))*SB*SC + (d^2 + 2*r^2 - R^2)*a^2*SA,
Q = 4*d^2*R*(4*r - R) + (d^2 - 3*R^2 + 4*r*(r + R))^2,
d = J*R, valJ, valr,
SA = (b^2 + c^2 - a^2)/2, SB = (-b^2 + c^2 + a^2)/2, SC = (b^2 - c^2 + a^2)/2, Q*fac);
collect(uuu,R, factor): valQ:= QQ= (factor@subs)(valR, %); normal(%-rot(%));
vvv:= (subs)(f(a,b,c) = (d^2 + (4*r - R)*R + sqrt(Q))*SB*SC + (d^2 + 2*r^2 - R^2)*a^2*SA,
sqrt(Q)=sqrt(QQ)/sqrtfac, d = J*R, valJ, valr,
SA = (b^2 + c^2 - a^2)/2, SB = (-b^2 + c^2 + a^2)/2, SC = (b^2 - c^2 + a^2)/2, f(a,b,c));
map(id^2, valR);
collect(vvv, R, factor); vvv1:= collect(subs(map(id^2, valR),map(1/id^2, valR), %), sqrt(QQ),normal);
vvv2:= collect(4*vvv1*a*b*c*(a+b+c), sqrt(QQ), factor);
remove(has,vvv2,QQ); a*b*c*collect(%/a/b/c,a, factor); insidetax(%);
vvv3:= parse(%)+select(has,vvv2,QQ);
subs(valQ,QQ): collect(%,a, factor): insidetax(%):
length(%); valQ:= QQ=parse(%%);
subs(valQ,vvv3); tmpp:=convert(%, string); L=length(%), ency(parse(%));
tmpp := "a*b*c*(2*a^7+(-3*b^2+4*b*c-3*c^2)*a^5+(b+c)*(b-c)^2*a^4-2*b*c*(b-c)^2*a^3-2*(b+c)*(b^2+b*c+c^2)*(b-c)^2*a^2+(b+c)^2*(b-c)^4*a+(b^2+c^2)*(b-c)^2*(b+c)^3)-2*(a^2-b^2+c^2)*(a^2+b^2-c^2)*(a+b+c)^(1/2)*R*(-a^9+(b+c)*a^8+(b^2+c^2)*a^7-(b+c)*(3*c^2-4*b*c+3*b^2)*a^6+(-c^2+2*b^2)*(-2*c^2+b^2)*a^5+(b+c)*(c^2-2*b*c+2*b^2)*(2*c^2-2*b*c+b^2)*a^4-3*(b^2+c^2)*(b^2-c^2)^2*a^3+(b+c)*(c^4+2*b*c^3-2*b^2*c^2+2*b^3*c+b^4)*(b-c)^2*a^2+(b^4+3*b^2*c^2+c^4)*(b^2-c^2)^2*a-(b-c)^2*(b^2-b*c+c^2)^2*(b+c)^3)^(1/2)"; L=length(%), ency(parse(%));
1117....
SubstituteAll(zalo[2],"<sub>",""): SubstituteAll(%,"</sub>","");
par1:= parse("n(a,b,c)=5*(SA)^2*((SB)^2+(SC)^2)-3*a^2*(SA)^3-4*(SB)^2*(SC)^2-SA*SB*SC*(2*SA-SB-SC)");
SubstituteAll(zalo[3],"<sub>",""): SubstituteAll(%,"</sub>","");
par2:= parse("d(a,b,c)=2*a*(4*(SA)^2-b^2*c^2)*(b^2*c^2*(3*a^2-8*SA)+8*(SA)^3)");
(factor@subs)(par1,par2,SA=(b^2+c^2-a^2)/2,SB=(-b^2+c^2+a^2)/2,SC=(b^2-c^2+a^2)/2,
a*n(a,b,c)/d(a,b,c)): elimifac(%);ency(%);
Trigo
3368
r\303\250gles
cos(A-2*pi/5): %=expand(%): rule51:= (eval@subs)(pi=Pi, %);
sin(A-2*pi/5): %=expand(%): rule52:= (eval@subs)(pi=Pi, -%);
sin(A+2*Pi/10): %=subs(rule5, expand(%)):
subs(A=A+Pi/10, rule51, rule52, rule5, %):
xcollect(%, [sin(A),cos(A)], factor):
rule53:= xcollect(%, [sin(A),cos(A),5+sqrt(5)], factor);
cos(A+2*Pi/10): %=subs(rule5, expand(%)):
subs(A=A+Pi/10, rule51, rule52, rule5, %):
xcollect(%, [sin(A),cos(A)], factor):
rule54:= xcollect(%, [sin(A),cos(A),5+sqrt(5)], factor);
sols:=[solve(x^24=1)]: nops(%);
map(evalc@Im, sols); evalf(%); evalf(sin(5*Pi/12));
cos(Pi/12)=(1/4)*sqrt(2)*sqrt(4+sqrt(2)*4^(1/4)*sqrt(3));
rule12:= simplify(%);
cos(5*Pi/12); simplify(%);
rule12 := { cos(1/12*Pi) = (1/4)*sqrt(2)*(1+sqrt(3)),sin(5/12*Pi)=cos(1/12*Pi),
cos(5/12*Pi) = sin(1/12*Pi),sin(1/12*Pi)=(1/4)*sqrt(2)*(-1+sqrt(3)) };
evalf(map(lhs-rhs,%));
www; subs(rule8, kashi, (eval@expand@convert)(parse(fdat[www]), sincos) );
collect( factor(%), R, factor); qqq:= subs(defR2pow,%); subs(rule8, %);
elimifac(numer(%))/elimifac(factor(denom(%)*sqrt(2+sqrt(2))));
tmqx:= xcollect( %, R, factor); ency(%);
subs(rule12, kashi, (expand@convert)(parse(fdat[www]), sincos) );
subs(rule12, %);
collect(%*R^2, R); subs(defR2pow,%); (elimifac@factor)(%):
tmqx:= map(collect, %, R, factor); ency(%);
www; subs(rule51, rule52, rule53, rule54, rule5, kashi, convert(parse(fdat[www]), sincos) ); subs(rule5, kashi, expand(%));
collect(%*R^2, R); subs(defR2pow,%); (elimifac@factor)(%):
tmqx:= map(xcollect, %, [R,4,20], factor); ency(%);
(factor@subs)(apbpc, tmqx): tmpx:= xcollect(elimifac(%), [R,3]); ency(%);
evalf(tmp); evalf(cos(2*Pi/5));
2*sin(Pi/5)*cos(Pi/5)=sin(2*Pi/5); subs(rules5,%);
isolate(%, sin((1/5)*Pi)); rationalize(%);
sin(A+Pi/10); cos(A-Pi/10);
2654
WW1,WW2,WW3,WW4;
subs(kitcircle, tmqx): (elimifac@factor)(%): tmqx:= xcollect(%, [WW1,WW2], (U->`if`(type(U,`*`), map(abctax,U), abctax(U)) )@factor );
2469, 2448
subs(vale, tmqx)*4*a*b*c*sqrt(a^2*b^2+a^2*c^2+b^2*c^2):
tmqx:= xcollect(%, [sqrt(a^6-a^4*b^2-a^4*c^2-a^2*b^4+3*a^2*b^2*c^2-a^2*c^4+b^6-b^4*c^2-b^2*c^4+c^6), sqrt(-a^2*c^2+c^4-b^2*c^2-a^2*b^2+a^4+b^4), A], factor);
ency(tmqx);
www;
tmqx: tmp1:= simplify(%) assuming a>0,b>0,c>0;
indets(tmp1, sqrt):
tmqx:= xcollect(tmp1*4*rotp(sqrt(a)*a)/R, [sqrt(a^3-a^2*b-a^2*c-a*b^2+3*a*b*c-a*c^2+b^3-b^2*c-b*c^2+c^3)], factor); ency(%);