Changeset 50cbdc in git for Singular/LIB/surf.lib
- Timestamp:
- Aug 27, 2001, 4:48:02 PM (23 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- 2567b5a6cb7109be5a83e53eb94abb1c38fb9945
- Parents:
- 3de58c9ca0aeaafdf5cb29f986967bffa405b542
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/surf.lib
r3de58c r50cbdc 1 1 //last change: 02.03.2001 (Eric Westenberger) 2 2 /////////////////////////////////////////////////////////////////////////////// 3 version="$Id: surf.lib,v 1. 19 2001-03-02 13:41:19 westenbExp $";3 version="$Id: surf.lib,v 1.20 2001-08-27 14:48:01 Singular Exp $"; 4 4 category="Visualization"; 5 5 info=" … … 40 40 41 41 proc plot(ideal I,list #) 42 "USAGE: plot(I); I ideal 43 ASSUME: I defines a plane curve or a surface 42 "USAGE: plot(I); I ideal or poly 43 ASSUME: I defines a plane curve or a surface given by one equation 44 44 RETURN: nothing 45 45 NOTE: requires the external program 'surf' to be installed … … 75 75 } 76 76 ideal I=simplify(phi(I),2); 77 if (leadcoef(I[1]) <0) { I[1]=-I[1]; } 77 78 if (ncols(I)==1 and n<=2) // curve 78 79 { … … 131 132 example 132 133 { "EXAMPLE:"; echo =2; 133 134 134 // --------- plane curves ------------ 135 ring rr0 = 0,(x1,x2),dp; 135 136 136 137 137 ideal I = x1^3 - x2^2; 138 plot(I); 138 139 139 140 141 140 ring rr1 = 0,(x,y,z),dp; 141 ideal I(1) = 2x2-1/2x3 +1-y+1; 142 plot(I(1)); 142 143 143 144 // ---- Singular Logo -------------- … … 145 146 plot(logo); 146 147 147 // --------- implicit curves ------------ 148 // implicit curves 149 ideal I(1) = y,-x2; 150 plot(I(1)); 151 ideal I(2) = x2,-y2 +4; 152 plot(I(2)); 148 // Steiner surface 149 ideal J(2) = x^2*y^2+x^2*z^2+y^2*z^2-17*x*y*z; 150 plot(J(2)); 153 151 154 //the lemniscate 155 156 ideal I(3) = x4+2x2y2 + y4, x2-y2; 157 plot(I(3)); 158 159 // ----------- surfaces ------------------- 160 ideal J(1) = 3xy4 + 2xy2, x5y3 + x + y6,10x2; 161 plot(J(1)); 162 163 // Steiner surface 164 165 ideal J(2) = x^2*y^2+x^2*z^2+y^2*z^2-17*x*y*z; 166 plot(J(2)); 167 152 // -------------------- 168 153 plot(x*(x2-y2)+z2); 169 154
Note: See TracChangeset
for help on using the changeset viewer.