Newton rapshon

clc
clear
m=0.8;
i=1;
a_m=[18;28];
delta(2,1)=0;
i_matrix(:,1)=0;
a1_matrix(:,1)=0;
a2_matrix(:,1)=0;
while true
jacob=[-sind(a_m(1,1))-sind(a_m(2,1));...
      -3*sind(3*a_m(1,1))-3*sind(3*a_m(2,1));...
      -5*sind(5*a_m(1,1))-5*sind(5*a_m(2,1));...
      -7*sind(7*a_m(1,1))-7*sind(7*a_m(2,1))]
      
b_func=[2*m;0];
f_func=[cosd(a_m(1,1))+cosd(a_m(2,1));...
     3*cosd(3*a_m(1,1))+3*cosd(3*a_m(2,1));..
     5*cosd(5*a_m(1,1))+5*cosd(5*a_m(2,1));...
     7*cosd(7*a_m(1,1))+7*cosd(7*a_m(2,1))]
delta= jacob/(b_func-f_func);
a_m=acosd(abs(cosd(a_m+delta)));
i_matrix(i,1)=i;
a1_matrix(i,1)=a_m(1,1); 
a2_matrix(i,1)=a_m(2,1);
if max(abs(delta))<0.057
break
end
i=i+1;
end
display(sort(a_m));
plot=(i_matrix,a1_matrix,i_matrix,a2_matrix
a=a_m(1,1);
b=a_m(2,1); 
t=0.02;