Changeset d406b6 in git
- Timestamp:
- Feb 18, 1999, 8:29:01 PM (25 years ago)
- Branches:
- (u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
- Children:
- d8e380c4561b7282d49138294559898e2ca63423
- Parents:
- 2e6e424cfca4721c43fc5841574211f5da6ef325
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/jordan.lib
r2e6e424 rd406b6 1 1 /////////////////////////////////////////////////////////////////////////////// 2 2 3 version="$Id: jordan.lib,v 1. 8 1999-02-11 11:47:45mschulze Exp $";3 version="$Id: jordan.lib,v 1.9 1999-02-18 19:29:01 mschulze Exp $"; 4 4 info=" 5 5 LIBRARY: jordan.lib PROCEDURES TO COMPUTE THE JORDAN NORMAL FORM … … 13 13 14 14 LIB "ring.lib"; 15 ///////////////////////////////////////////////////////////////////////////////16 17 ///////////////////////////////////////////////////////////////////////////////18 // PROCEDURE: jordan PROCEDURES TO COMPUTE THE JORDAN NORMAL FORM //19 // AUTHOR: Mathias Schulze, email: mschulze@mathematik.uni-kl.de //20 15 /////////////////////////////////////////////////////////////////////////////// 21 16 … … 208 203 } 209 204 example 210 { 211 "EXAMPLE:"; 212 echo=2; 205 { "EXAMPLE:"; echo=2; 213 206 ring R=0,x,dp; 214 207 matrix M[3][3]=3,2,1,0,2,1,0,0,3; … … 218 211 /////////////////////////////////////////////////////////////////////////////// 219 212 220 ///////////////////////////////////////////////////////////////////////////////221 // PROCEDURE: jordanmatrix PROCEDURES TO COMPUTE THE JORDAN NORMAL FORM //222 // AUTHOR: Mathias Schulze, email: mschulze@mathematik.uni-kl.de //223 ///////////////////////////////////////////////////////////////////////////////224 225 213 proc jordanmatrix(list l) 226 214 "USAGE: jordanmatrix(l); l list of ideal and list of intvecs 227 RETURN: the Jordan matrix J with215 RETURN: the Jordan matrix J defined by l: 228 216 l[1] : eigenvalues of J 229 217 l[2][i][j] : size of j-th Jordan block with eigenvalue l[1][i] … … 305 293 } 306 294 example 307 { 308 "EXAMPLE:"; 309 echo=2; 295 { "EXAMPLE:"; echo=2; 310 296 ring R=0,x,dp; 311 print(jordanmatrix(list(ideal(2,3),list(intvec(1),intvec(2))))); 312 } 313 /////////////////////////////////////////////////////////////////////////////// 314 315 /////////////////////////////////////////////////////////////////////////////// 316 // PROCEDURE: jordanform PROCEDURES TO COMPUTE THE JORDAN NORMAL FORM // 317 // AUTHOR: Mathias Schulze, email: mschulze@mathematik.uni-kl.de // 297 list l; 298 l[1]=ideal(2,3); 299 l[2]=list(intvec(1),intvec(2)); 300 print(jordanmatrix(l)); 301 } 318 302 /////////////////////////////////////////////////////////////////////////////// 319 303 … … 329 313 } 330 314 example 331 { 332 "EXAMPLE:"; 333 echo=2; 315 { "EXAMPLE:"; echo=2; 334 316 ring R=0,x,dp; 335 317 matrix M[3][3]=3,2,1,0,2,1,0,0,3; … … 337 319 print(jordanform(M)); 338 320 } 339 ///////////////////////////////////////////////////////////////////////////////340 341 ///////////////////////////////////////////////////////////////////////////////342 // PROCEDURE: inversemat PROCEDURES TO COMPUTE THE JORDAN NORMAL FORM //343 // AUTHOR: Mathias Schulze, email: mschulze@mathematik.uni-kl.de //344 321 /////////////////////////////////////////////////////////////////////////////// 345 322 … … 369 346 } 370 347 example 371 { 372 "EXAMPLE:"; 373 echo=2; 348 { "EXAMPLE:"; echo=2; 374 349 ring R=0,x,dp; 375 350 matrix M[3][3]=3,2,1,0,2,1,0,0,3; … … 378 353 } 379 354 /////////////////////////////////////////////////////////////////////////////// 355
Note: See TracChangeset
for help on using the changeset viewer.