force_analytic_test.m~ 358 Bytes
Newer Older
ppanchal's avatar
ppanchal committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
% Testing the force formula

% Testing convergence for some alpha values
alpha1 = acosh(10);
alpha2 = acosh(20);

force1 = [];
force2 = [];

for i = 1 : 500
   force1 = [force1 Gsum(i,alpha1)];
   force2 = [force2 Gsum(i,alpha2)];
end

figure;

plot(1:500,force1);
hold on;
plot(1:500,force2);

%alpha = acosh(1+x);

% Next plot

xvals = 0.001:0.001:0.01;