3/29/23 4:15 PM C:\Users\14693\Desktop\...\MECH413hw10.m 1 of 2
%%
%%newton raphson method
clc; %clear command window
clear all; %clear all existing variables
close all; %close all open figures
h = 10;
a = 1;
n = 1;
diff = 1;
while (diff > 0.001)
F = 91 - 40*sqrt(pi*a)*(1.1216+6.5200*(a/h)^2-12.3877*(a/h)^4+89.0554*(a/h)^6-188.6080*
(a/h)^8+207.3870*(a/h)^10-32.0524*(a/h)^12);
F_prime = -40*sqrt(pi)*(1.1216*0.5*a^-0.5+6.5200*2.5*(a^1.5)/(h^2)-12.3877*4.5*(a^3.5)/
(h^4)+89.0554*6.5*(a^5.5)/(h^6)-188.6080*8.5*(a^7.5)/(h^8)+207.3870*10.5*(a^9.5)/(h^10)
-32.0524*12.5*(a^11.5)/(h^12));
anew = a - F/F_prime;
diff = abs(anew - a)/a;
a = anew;
n = n + 1;
end
%n = 3
%a = 1.1366
%%
clc; %clear command window
clear all; %clear all existing variables
close all; %close all open figures
a = 0.1;
b = 1.13658;
alpha = 2.6;
h = 10;
c = 1.6*10^-9;
sigma = 40;
xi1 = 0.0000000000000000;
xi2 = -0.2011940939974345;
xi3 = 0.2011940939974345;
xi4 = -0.3941513470775634;
xi5 = 0.3941513470775634;
xi6 = -0.5709721726085388;
xi7 = 0.5709721726085388;
xi8 = -0.7244177313601701;
xi9 = 0.7244177313601701;
