                   Grupy punktowe w krystalografii: przekształcenia położenia punktu: Zadanie 1
                   © W.J. 2008-2016 

Procedury i definicje systemowe:

kill(all)$

load(draw)$

number: true ; ratprint: false;

sa: float(sin(%pi/3))$ ca: float(cos(%pi/3))$ sa2: float(sin(2*%pi/3))$ ca2: float(cos(2*%pi/3))$

parametry1: [
line_width=1,
user_preamble="set size ratio 1",
color=black,
head_type = 'nofilled, head_length = .05, 
vector([-1,0],[2.1,0]), text:label(["y",1.1,-0.1]),
vector([0,1],[0,-2.1]), text:label(["x",-0.1,-1.1]),
implicit( x^2+y^2=1, x,-1.20,1.2, y,-1.2,1.2 ),

line_width=2,
implicit( (x-1)^2+(y-0)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),
implicit( (x-0)^2+(y-1)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),
implicit( (x-0)^2+(y+1)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),
implicit( (x+1)^2+(y-0)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),

implicit( (x-sa)^2+(y-ca)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),
implicit( (x-sa2)^2+(y-ca2)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),
implicit( (x-ca)^2+(y+sa)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),
implicit( (x-ca2)^2+(y+sa2)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),

implicit( (x+sa)^2+(y-ca)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),
implicit( (x+sa2)^2+(y-ca2)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),
implicit( (x+ca)^2+(y-sa)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),
implicit( (x+ca2)^2+(y-sa2)^2=0.0005, x,-1.2,1.2, y,-1.2,1.2 ),

line_width=1
] $

Definicje macierzy elementów symetrii:

L1:      matrix( [ 1, 0, 0], [ 0, 1, 0], [ 0, 0, 1]) $
C:       matrix( [-1, 0, 0], [ 0,-1, 0], [ 0, 0,-1]) $
P001:    matrix( [ 1, 0, 0], [ 0, 1, 0], [ 0, 0,-1]) $
P010:    matrix( [ 1, 0, 0], [ 0,-1, 0], [ 0, 0, 1]) $
P100:    matrix( [-1, 0, 0], [ 0, 1, 0], [ 0, 0, 1]) $
P110:    matrix( [ 0,-1, 0], [-1, 0, 0], [ 0, 0, 1]) $
P011:    matrix( [ 1, 0, 0], [ 0, 0,-1], [ 0,-1, 0]) $
P101:    matrix( [ 0, 0,-1], [ 0, 1, 0], [-1, 0, 0]) $
Pm110:   matrix( [ 0, 1, 0], [ 1, 0, 0], [ 0, 0, 1]) $
Pm101:   matrix( [ 0, 0, 1], [ 0, 1, 0], [ 1, 0, 0]) $
P0m11:   matrix( [ 1, 0, 0], [ 0, 0, 1], [ 0, 1, 0]) $
L2001:   P001.C $
L2010:   P010.C $
L2100:   P100.C $
L2110:   P110.C $
L2011:   P011.C $
L2101:   P101.C $
L2m110:  Pm110.C $
L2m101:  Pm101.C $ 
L20m11:  P0m11.C $
L4100:   matrix( [ 1, 0, 0], [ 0, 0,-1], [ 0, 1, 0]) $
L4010:   matrix( [ 0, 0, 1], [ 0, 1, 0], [-1, 0, 0]) $
L4001:   matrix( [ 0,-1, 0], [ 1, 0, 0], [ 0, 0, 1]) $
A4001:   L4001.C $
A4010:   L4010.C $
A4100:   L4100.C $

 
Zadanie 1: Oblicz rachunkiem macierzowym wynik działania kolejno trzech elementów symetrii na punkt P1.
 

Pkt1: matrix([0.7], [0.2], [0.4]) $

I teraz 10 przykładów działania na punkt.
Uwaga, kolejność działań jest nieprzemienna, 
tu najpierw działa oś lub płaszczyzna diagonalna, drugi to płaszczyzna prostopadła do osi "z" a trzeci to centrum. 
Przykłady różnią się pierwszym elementem.

Najpierw same wyniki liczbowe:

[ C.P001.L2m110.Pkt1, 
  C.P001.L2011.Pkt1, 
  C.P001.L2101.Pkt1, 
  C.P001.L2m101.Pkt1, 
  C.P001.L20m11.Pkt1 ];

[ C.P001.Pm110.Pkt1, 
  C.P001.P011.Pkt1, 
  C.P001.P101.Pkt1, 
  C.P001.Pm101.Pkt1, 
  C.P001.P0m11.Pkt1 ];

Osie zwyke

1. L2m110

[ NPkt1: Pkt1,
  NPkt2: L2m110.NPkt1,
  NPkt3: P001.L2m110.NPkt1,
  NPkt4: C.P001.L2m110.NPkt1 ];

wxdraw2d(

parametry1,

if NPkt1[3,1]>0 then color= red else color=blue,
if NPkt1[3,1]>0 then 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),

if NPkt2[3,1]>0 then color= red else color=blue,
if NPkt2[3,1]>0 then 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),

if NPkt3[3,1]>0 then color= red else color=blue,
if NPkt3[3,1]>0 then 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),

if NPkt4[3,1]>0 then color= red else color=blue,
if NPkt4[3,1]>0 then 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 )

);

2. L2011

[ NPkt1: Pkt1,
  NPkt2: L2011.NPkt1,
  NPkt3: P001.L2011.NPkt1,
  NPkt4: C.P001.L2011.NPkt1 ];

wxdraw2d(
parametry1,
if NPkt1[3,1]>0 then color= red else color=blue,
if NPkt1[3,1]>0 then 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt2[3,1]>0 then color= red else color=blue,
if NPkt2[3,1]>0 then 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt3[3,1]>0 then color= red else color=blue,
if NPkt3[3,1]>0 then 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt4[3,1]>0 then color= red else color=blue,
if NPkt4[3,1]>0 then 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 )
);

3. L2101

[ NPkt1: Pkt1,
  NPkt2: L2101.NPkt1,
  NPkt3: P001.L2101.NPkt1,
  NPkt4: C.P001.L2101.NPkt1 ];

wxdraw2d(
parametry1,
if NPkt1[3,1]>0 then color= red else color=blue,
if NPkt1[3,1]>0 then 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt2[3,1]>0 then color= red else color=blue,
if NPkt2[3,1]>0 then 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt3[3,1]>0 then color= red else color=blue,
if NPkt3[3,1]>0 then 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt4[3,1]>0 then color= red else color=blue,
if NPkt4[3,1]>0 then 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 )
);

4. L2m101

[ NPkt1: Pkt1,
  NPkt2: L2m101.NPkt1,
  NPkt3: P001.L2m101.NPkt1,
  NPkt4: C.P001.L2m101.NPkt1 ];

wxdraw2d(
parametry1,
if NPkt1[3,1]>0 then color= red else color=blue,
if NPkt1[3,1]>0 then 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt2[3,1]>0 then color= red else color=blue,
if NPkt2[3,1]>0 then 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt3[3,1]>0 then color= red else color=blue,
if NPkt3[3,1]>0 then 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt4[3,1]>0 then color= red else color=blue,
if NPkt4[3,1]>0 then 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 )
);

5. L20m11

[ NPkt1: Pkt1,
  NPkt2: L20m11.NPkt1,
  NPkt3: P001.L20m11.NPkt1,
  NPkt4: C.P001.L20m11.NPkt1 ];

wxdraw2d(
parametry1,
if NPkt1[3,1]>0 then color= red else color=blue,
if NPkt1[3,1]>0 then 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt2[3,1]>0 then color= red else color=blue,
if NPkt2[3,1]>0 then 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt3[3,1]>0 then color= red else color=blue,
if NPkt3[3,1]>0 then 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt4[3,1]>0 then color= red else color=blue,
if NPkt4[3,1]>0 then 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 )
);

Płaszczyzny

1. Pm110

[ NPkt1: Pkt1,
  NPkt2: Pm110.NPkt1,
  NPkt3: P001.Pm110.NPkt1,
  NPkt4: C.P001.Pm110.NPkt1 ];

wxdraw2d(
parametry1,
if NPkt1[3,1]>0 then color= red else color=blue,
if NPkt1[3,1]>0 then 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt2[3,1]>0 then color= red else color=blue,
if NPkt2[3,1]>0 then 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt3[3,1]>0 then color= red else color=blue,
if NPkt3[3,1]>0 then 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt4[3,1]>0 then color= red else color=blue,
if NPkt4[3,1]>0 then 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 )
);

2. P011

[ NPkt1: Pkt1,
  NPkt2: P011.NPkt1,
  NPkt3: P001.P011.NPkt1,
  NPkt4: C.P001.P011.NPkt1 ];

wxdraw2d(
parametry1,
if NPkt1[3,1]>0 then color= red else color=blue,
if NPkt1[3,1]>0 then 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt2[3,1]>0 then color= red else color=blue,
if NPkt2[3,1]>0 then 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt3[3,1]>0 then color= red else color=blue,
if NPkt3[3,1]>0 then 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt4[3,1]>0 then color= red else color=blue,
if NPkt4[3,1]>0 then 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 )
);

3. P101

[ NPkt1: Pkt1,
  NPkt2: P101.NPkt1,
  NPkt3: P001.P101.NPkt1,
  NPkt4: C.P001.P101.NPkt1 ];

wxdraw2d(
parametry1,
if NPkt1[3,1]>0 then color= red else color=blue,
if NPkt1[3,1]>0 then 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt2[3,1]>0 then color= red else color=blue,
if NPkt2[3,1]>0 then 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt3[3,1]>0 then color= red else color=blue,
if NPkt3[3,1]>0 then 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt4[3,1]>0 then color= red else color=blue,
if NPkt4[3,1]>0 then 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 )
);

4. Pm101

[ NPkt1: Pkt1,
  NPkt2: Pm101.NPkt1,
  NPkt3: P001.Pm101.NPkt1,
  NPkt4: C.P001.Pm101.NPkt1 ];

wxdraw2d(
parametry1,
if NPkt1[3,1]>0 then color= red else color=blue,
if NPkt1[3,1]>0 then 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt2[3,1]>0 then color= red else color=blue,
if NPkt2[3,1]>0 then 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt3[3,1]>0 then color= red else color=blue,
if NPkt3[3,1]>0 then 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt4[3,1]>0 then color= red else color=blue,
if NPkt4[3,1]>0 then 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 )
);

5. P0m11

[ NPkt1: Pkt1,
  NPkt2: P0m11.NPkt1,
  NPkt3: P001.P0m11.NPkt1,
  NPkt4: C.P001.P0m11.NPkt1 ];

wxdraw2d(
parametry1,
if NPkt1[3,1]>0 then color= red else color=blue,
if NPkt1[3,1]>0 then 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt1[2,1])^2 + (y-(-NPkt1[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt2[3,1]>0 then color= red else color=blue,
if NPkt2[3,1]>0 then 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt2[2,1])^2 + (y-(-NPkt2[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt3[3,1]>0 then color= red else color=blue,
if NPkt3[3,1]>0 then 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt3[2,1])^2 + (y-(-NPkt3[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 ),
if NPkt4[3,1]>0 then color= red else color=blue,
if NPkt4[3,1]>0 then 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.005, x,-1.2,1.2, y,-1.2,1.2 ) 
 else 
  implicit( (x-NPkt4[2,1])^2 + (y-(-NPkt4[1,1]))^2=0.002, x,-1.2,1.2, y,-1.2,1.2 )
);