newtons.m

fem nonlinear source code using weighted residuals

function [c]=Newtons(c,R,x,R1,R2,N,alpha)

for COUNTER = 1:100
    [J]=Jacobian(x,c,alpha,N);
    [R]=Residual(x,c,alpha,N,R,R1,R2);

    delta_c = J\(-R);
    c = c+delta_c;

    normRatio = norm(delta_c)/norm(c);

    if (normRatio)< 1E-1
        disp ('Newtons works');
        break
    end

    if COUNTER ==100
        disp('Newton FAIL');
        break
    end
end

end