This is the part where the error occurs and it reads : "Error using ./ Matrix dimensions must agree."
These are the sizes of F2 and J2.
size(F2)ans = 122 1 size(J2)ans = 120 122
Is this possible to divide a matrix by a vector like this? J2 also has many 0 components. Any help will be appreciated. Thanks!
%20 newtons iteration
for i=20[F1,F2,J1,J2]=thomas_newton(s,tm); delta_s=-F2'./J2'; s=s+delta_s;if norm(delta_s)<tol breakendend
Best Answer