Code Directory
 Visual Basic & VB.NET
New Code
White-label Grocery Delivery App Solution 2.0
Best Classified Script 5.1
Readymade B2B Script 1.3.1
Uber Clone with Safety Measure Addons 2.0
Equity Investing Software 1.3.2
C# QR Code Generator 2020.6.0.0
The .Net PDF Library 2020.7.1
dbExpress driver for MySQL 7.2
dbForge Documenter for Oracle 1.2
dbForge Studio for Oracle 4.2
Excel .Net Library 2020.6
fsMediaLibrary.NET 2019.11
VaxVoIP SIP Server SDK 5.2.0
Database Workbench Pro 5.7.4
dbForge Data Generator for Oracle 2.2
Top Code
dbExpress driver for MySQL 3.00
Database Workbench Pro 5.7.4
Availability Booking Calendar PHP 1.0
dbExpress driver for SQL Server 4.45
ATN Site Builder 3.0
dbForge Studio for Oracle 3.10
MLM Software ONE 1.5.46
IcrediBB Bulletin Board System 1.0
Invoice Manager by PHPJabbers 3.0
MLM Unilevel Plan Software 1.0.2
ICPennyBid Penny Auction Script 4.0
Paste phpSoftPro 1.4.1
Simple Web Content Management System for Scripts 1.1
ATN Resume Finder 2.0
PHP Review Script 1.0
Top Rated
phpEnter 5.1.
Single Leg MLM 1.2.1
Azizi search engine script PHP 4.1.10
Paste phpSoftPro 1.4.1
Extreme Injector 3.7
Deals and Discounts Website Script 1.0.2
Solid File System OS edition 5.1
Classified Ad Lister 1.0
Aglowsoft SQL Query Tools 8.2
Invoice Manager by PHPJabbers 3.0
ICPennyBid Penny Auction Script 4.0
PHP Review Script 1.0
ATN Resume Finder 2.0
ATN Site Builder 3.0
Availability Booking Calendar PHP 1.0
GCD of Polynomials and Polynomials raised to some Power including Fractional Power 1.0
File ID: 82530

GCD of Polynomials and Polynomials raised to some Power including Fractional Power 1.0
Download GCD of Polynomials and Polynomials raised to some Power including Fractional Power 1.0http://www.mathworks.comReport Error Link
License: Shareware
File Size: 92.2 KB
Downloads: 8
Submit Rating:
GCD of Polynomials and Polynomials raised to some Power including Fractional Power 1.0 Description
Description: 15th July, 2005 : Poly_POWER.m is now corrected !

So, for most reasonable cases of Multiple Roots including Multiple Real Roots, this Programme should now work.

For eg, Poly_POWER works successfully for :
{ x^6 + (-12+18j)*x^5 + (-75-180j)*x^4 + (920+180j)*x^3 + (-1785+1800j)*x^2 + (-732-3582j )*x + (2035+828j) } ^ 0.5
The answer is "approx":
{ ( x^3 + (-6 + 9j) * x^2 + (-15 - 36j) * x + (46 + 9j) }
and for Liouville's Constant based Polys, like :
( { x^6 - 75*x^3 - 190*x + 21 } ^ 3 } ^ 0.3333


Functional Description of Poly_GCD :
If we need to verify the fact that a Polynomial has multiple roots iff (if and only if) it has a common factor with it's derivative, we need two things :

a) A function to compute the GCD of 2 Polynomials. Since I could not find a Standard Matlab function for this, I created this function : Poly_GCD.m :
GCD = Poly_GCD ( sx_poly, rx_poly )
For eg, if :
sx_poly = x^9 - 3x^8 + 0x^7 + 2x^6 + 6x^5 + 0x^4 - 4x^3 - 6x^2 - 3x ? 1
and rx_poly = 9x^8 - 24x^7 + 0x^6 + 12x^5 + 30x^4 + 0x^3 - 12x^2 - 12x ? 3,
then their GCD = x^6 - 2x^5 - x^4 + 3x^2 + 2x + 1

b) A function which can raise a Polynomial to some power in order to simulate an overall polynomial with multiple roots. Since I could not find a Standard Matlab function for this, I created Poly_POWER.m :
P = Poly_POWER ( poly, n )

In addition to computing the GCD of the 2 input Polynomials : sx_Poly and rx_Poly, Poly_GCD.m and Poly_GCD_Main.m also find suitable polynomials :
Nx & nx and cx & dx such that :
check_GCD_Orig = Nx * sx_Poly + nx * rx_Poly
check_GCD_Used = cx * sx_Poly_Used + dx * rx_Poly_Used (Internally used vars)

While finding HCF of 2 Polynomials, as the degrees of the polynomials increase, the accumulated FP errors increase, and can only be "salvaged" to some extent. The trick lies in being able to devise a way to detect the correct Zero limit on the Remainder, and thereby, stop after the correct number of HCF divisions. This requires lots of experiments to find what limits to apply for the definition of 0 (zero) ; this has been done with some complex empirical logics obtained by trial and error.

Convergence also depends upon whether we convert the input Polys to monic, or the intermediate computed Poly also to monic, or we do not convert at all. These 3 combinations create a whole lot of extra logics, almost 30 % of the development effort ! We choose the best combination.
When compared to my very early submission of this zip file, the earlier Poly_GCD.m is now renamed as Poly_GCD_Main.m, and Poly_GCD.m is now a "top level" function.

I think that all these functions : Poly_GCD, Poly_POWER, CMPLX_GCD, Ch_Rem_Thr_Poly.m, Ch_Rem_Thr_Int.m, Gen_Primes_Eq_2_Sqs etc are essential as commonly required functions, and can be placed in the dir : ...matlabspecfun

Should generally work for R14, R13 and R12.

License: Shareware

Related: apply, definition, Complex, Limits, experiments, divisions, requires, empirical, logics, depends

O/S:BSD, Linux, Solaris, Mac OS X

File Size: 92.2 KB

Downloads: 8

More Similar Code

chatFX combines the power of Flash and XML to produce an intuitive, multi-platform, highly customizable chat client anyone can install and use.

XSQL combines the power of XML and SQL to provide a language and database independent means to store and retrieve SQL queries and their results.

Sometimes you get a list of lists and want to swap rows and columns, i.e. transpose the list. Yet, what if the rows have different lengths? Here is some advice you might find useful in such situations.

BTL is a template language that combines power of JSTL and XSLT to produce documents in XML, HTML, XHTML, XSL-FO, PDF or other formats, based on the JavaBean input.

The purpose of project is libarary that is foundation for importers from Rational Rose file formats to other formats and importers to some specific formats like nsuml xmi.

It is application of of Simulink Block to Cardiac PDE VI1 system of 2 non linear coupled equations of PDE and try to design solution with dicretized space and time and issue in 1D.I associate a Doc file to introduction theory.Just see this page...

These algorithms are an implementation of the iterative amplitude adapted Fourier Transform method of Schreiber and Schmitz to make surrogate time series. These surrogate time series or fields were made to have the Power spectrum and the amplitude...

Simios is a collection of libraries and tools to ease development on any CLI framwork, focusing in mono and MS.NET. All libraries are 100% managed and are beeing developed by the community with some support by other Mono communities as...

SVGMap is an open source web tool to help you in the visualisation of experimental data associated to some graphical representation.

More at :

1. Project website

2. Reference :

DSPack is a set of Components and class to write Multimedia Applications using MS Direct Show and DirectX technologies. DSPack 2.2 is designed to work with DirectX 9 on Win9X, ME, 2000, and Windows XP operating systems. DSPack 2.2 is designed to...

User Review for GCD of Polynomials and Polynomials raised to some Power including Fractional Power
- required fields

Please enter text on the image