hw06_p04_newton - f_y = 3.0*xold*yold^2; f function g_x =...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
function [xnew, ynew] = hw06_p04_newton( xold, yold) A(1,1) = fx(xold, yold); A(1,2) = fy(xold, yold); A(2,1) = gx(xold, yold); A(2,2) = gy(xold, yold); A b(1,1) = -f(xold, yold); b(2,1) = -g(xold, yold); b delta = A\b; d xnew = xold + delta(1,1); ynew = yold + delta(2,1); y function f_x = fx( xold, yold ) f_x = 2.0*xold + yold^3; f function f_y = fy( xold, yold )
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: f_y = 3.0*xold*yold^2; f function g_x = gx( xold, yold ) g_x = 6.0*xold*yold; g function g_y = gy( xold, yold ) g_y = 3.0*xold^2 - 3.0*yold^2; g function f_ = f( xold, yold ) f_ = xold^2 + xold*yold^3 - 9.0; f function g_ = g( xold, yold ) g_ = 3.0*xold^2*yold - yold^3 - 4.0;...
View Full Document

This note was uploaded on 02/22/2010 for the course CHE 348 taught by Professor Chelikowsky during the Spring '08 term at University of Texas.

Ask a homework question - tutors are online