Relationship between latitude and gravitational acceleration

% 基本参数
a = 6.378*10^6;
b = 6.356*10^6;
M = 5.988*10^24;
G = 6.67*10^-11;
w = 7.27*10^-5;

%纬度转化：度分秒转化为弧度制
theta1 = 35;
theta2 = 56;
theta3 = 18;
theta = theta1+theta2/60+theta3/3600;
angle = pi*theta/180;

% 具体公式
n = (a*cos(angle))^2+(b*sin(angle))^2;
A=(G*M)^2/(n^2);
B=n*((cos(angle))^2)*(w^4);
C=((2*G*M*(w^2)/(sqrt(n)))*(cos(angle)^2));
g=sqrt(A+B-C);
（引流一下---B站up主：学物理的兔八哥）