ideals.log 1.7 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156
  1. Codemist Standard Lisp 3.54 for DEC Alpha: May 23 1994
  2. Dump file created: Mon May 23 10:39:11 1994
  3. REDUCE 3.5, 15-Oct-93 ...
  4. Memory allocation: 6023424 bytes
  5. +++ About to read file tstlib.red
  6. I_setting(x,y,z);
  7. torder revgradlex;
  8. lex
  9. u := I(x*z-y**2, x**3-y*z);
  10. 2 3
  11. u := i(x*z - y ,x - y*z)
  12. y member I(x,y^2);
  13. 0
  14. x member I(x,y^2);
  15. 1
  16. I(x,y^2) subset I(x,y);
  17. 1
  18. % yes
  19. I(x,y) subset I(x,y^2);
  20. 0
  21. % no
  22. % examples taken from Cox, Little, O'Shea: "Ideals, Varieties and Algorithms"
  23. q1 := u .: I(x);
  24. 3 2 2 2
  25. q1 := i(x - y*z,x *y - z , - x*z + y )
  26. % quotient ideal
  27. q2 := u .+ I(x^2 * y - z^2);
  28. 3 2 2 2
  29. q2 := i(x - y*z,x *y - z , - x*z + y )
  30. % sum ideal
  31. if q1 .= q2 then write "same ideal";
  32. same ideal
  33. % test equality
  34. intersection(u,I(y));
  35. 3 2 2 2 2 3
  36. i(x *y - y *z,x *y - y*z , - x*y*z + y )
  37. % ideal intersection
  38. u .: I(y);
  39. 3 2 2 2
  40. i(x - y*z,x *y - z , - x*z + y )
  41. u .: I(x,y);
  42. 3 2 2 2
  43. i(x - y*z,x *y - z , - x*z + y )
  44. %-----------------------------------------------------
  45. u1 := I(x,y^2);
  46. 2
  47. u1 := i(x,y )
  48. u1u1:= u1 .* u1;
  49. 4 2 2
  50. u1u1 := i(y ,x*y ,x )
  51. % square ideal
  52. u0 :=I(x,y);
  53. u0 := i(x,y)
  54. % test equality/inclusion for u1,u1u1,u0
  55. u1 .= u1u1;
  56. 0
  57. % no
  58. u1 subset u1u1;
  59. 0
  60. % no
  61. u1u1 subset u1;
  62. 1
  63. % yes
  64. u1 .= u0;
  65. 0
  66. % no
  67. u1 subset u0;
  68. 1
  69. % yes
  70. intersection (I(x) , I(x^2,x*y,y^2)) .= intersection(I(x) , I(x^2,y));
  71. 1
  72. end;
  73. (ideals 1100 0)
  74. End of Lisp run after 1.11+0.81 seconds