q, a vector of 5 complex coefficients of a quartic
xt, a vector of 3 x grid values, xlo, xinc and xhi
yt, a vector of 3 y grid values, ylo, yinc and yhi
maxiter, the max number of iterations
and then paints the portion of the complex plane (where xlo< x< xhi and
ylo< y< yhi) with points colored by root. Your qnewt must follow the
meshgrid and find procedure sketched on the lecture page.
Your qnewt must also make use of your
homegrown polyder function. Your qnewt must automatically label its
portrait with the user defined polynomial, as in the portrait below, achieved via
qnewt([1 0 -0.84 0 -0.16],[.45 .0001 .55],[-.05 .0001 .05],20)
Submit your fully documented code along with 3 portraits of 3 distinct quartics - 1 of which looks like the one at right.
First draft of qnewt.m
2 pts for proper use of matlab's roots function
2 pts for sketch of myownpolyder subfunction
2 pts for sketch of qlab
Final draft of qnewt.m
12 pts for headers
please model your headers after
seurat. Most importantly, define
only those variables that go in or
out of the function in question. All other variables
can be documented at time of first usage.
5 pts for further comments in code
5 pts for indentation
14 pts for correct myownpolyder
12 pts for correct qlab (no 1s or 0s or +-)
6 pts for three portraits, titled with the user specified quartic.
The first plot should be EXACTLY like that above.