hw06_p03a_newton

# hw06_p03a_newton - f function f_y = fy xold yold f_y =...

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

function [xnew, ynew] = hw06_p03a_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; xnew = xold + delta(1,1); ynew = yold + delta(2,1); y function f_x = fx( xold, yold ) f_x = 2.0*xold;
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: f function f_y = fy( xold, yold ) f_y = 2.0*yold; f function g_x = gx( xold, yold ) g_x = 2.0*xold; g function g_y = gy( xold, yold ) g_y = -2.0*yold; g function f_ = f( xold, yold ) f_ = xold^2 + yold^2 - 4.0; f function g_ = g( xold, yold ) g_ = xold^2 - yold^2 - 1.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