Blame view

Giac_maj/giac-1.4.9/examples/geo/equimorley 1008 Bytes
6663b6c9   adorian   projet complet av...
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
  equimore(m,n):={
  assume(a1=0.3); 
  //[0.00597664509739,0.00303053313965]
  assume(a2=0.4);
  //(a1+a2+a2=pi/3);
  //a1:=0.2;
  //a2:=0.4;
  A:=[0,1+i*texpand(tan(a1)),1+i*texpand(tan(2*a1)),
  1+i*texpand(tan(pi/3+a1)),1+i*texpand(tan(2*a1+2*pi/3)),
  1+i*texpand(tan(a1+2*pi/3)),1+i*texpand(tan(pi/3+2*a1))];
  B:=[1,i*texpand(tan(2*a2)),i*texpand(tan(a2)),
  i*texpand(tan(2*a2+2*pi/3)),i*texpand(tan(a2+pi/3)),
  i*texpand(tan(pi/3+2*a2)),i*texpand(tan(2*pi/3+a2))];
  C0:=texpand(tan(a2*3)/(tan(a1*3)+tan(a2*3))*(1+i*tan(a1*3)));
  C:=[C0,C0+1+i*texpand(tan(pi/3+2*a1-a2)),C0+1+i*texpand(tan(2*pi/3+a1-2*a2)),
  C0+1+i*texpand(tan(2*pi/3+2*a1-a2)),C0+1+i*texpand(tan(pi/3+a1-2*a2)),
  C0+1+i*texpand(tan(2*a1-a2)),C0+1+i*texpand(tan(a1-2*a2))];
  P:=affixe((inter(droite(A[0],A[m]),droite(B[0],B[n])))[0]);
  R:=affixe((inter(droite(A[0],A[n]),droite(C[0],C[m])))[0]);
  Q:=affixe((inter(droite(B[0],B[m]),droite(C[0],C[n])))[0]);
  lpq:=longueur2(P,Q);
  lpr:=longueur2(P,R);
  lqr:=longueur2(Q,R);
  return([normal(lpq-lpr),normal(lpq-lqr)]);
  };