The error is:
Warning: 362 equations in 1 variables. > In /Applications/MATLAB_R2012a.app/toolbox/symbolic/symbolic/symengine.p>symengine at 54 In mupadengine.mupadengine>mupadengine.evalin at 97 In mupadengine.mupadengine>mupadengine.feval at 150 In solve at 160 In Hyper_HW2_P3 at 39 Warning: Explicit solution could not be found. > In solve at 169 In Hyper_HW2_P3 at 39
The code is:
% Given parameters of geoid for Earth U = 6.263685953707e7; % potential [m^2/s^2] GM = 3.986005e14; % gravitational constant [m^3/s^2] J_2 = 1.08263e-3; % Jeffrey constants J_3 = 2.532153e-7; J_4 = 1.6109876e-7; w = 7.292115147e-5; % angular acceleration [rad/s] R_e = 6.378135e6; % equatorial radius [m] R_p = 6.3567506e6; % polar radius [m] phi = [0:180]; % angles of colatitude [deg]
% Relations for Legendre Polynomials for scalar potential equation P_2 = (1/2).*(3.*cosd(phi).^2-1); P_3 = (1/2).*(5.*cosd(phi).^3 - 3.*cosd(phi)); P_4 = (1/8).*(35.*cosd(phi).^3 - 30.*cosd(phi).^2 + 3);
% Determine R_RE from ellipse equation R_RE = (R_e*R_p)./(sqrt(R_p^2.*sind(phi).^2 + R_e^2.*cosd(phi).^2));
% Determine potential from the centrifugal force to include the effects of % rotation for the geoid syms R_GEOID
U_c = -(1/2).*w^2.*R_GEOID^2.*sind(phi).^2;
% Analytically solve for R_GEOID from potential function
P2 = (R_e/R_GEOID)^2*J_2*P_2; P3 = (R_e/R_GEOID)^3*J_3*P_3; P4 = (R_e/R_GEOID)^4*J_4*P_4;
S = solve(-(1/2).*w^2.*R_GEOID.^2.*sind(phi).^2 == U_c,... (GM./R_GEOID).*(1 - (P2 + P3 + P4)) == U )
No products are associated with this question.
When you solve() a matrix, it tries to solve all of the entries simultaneously. It does not consider the elements to be a set of related problems that are to be solved for individually with an array of different results.
What you should do is leave phi as a symbolic parameter instead of making it numeric, and solve() for that, giving you a formula (or set of formulae) that is parameterized on phi. You can then subs() the actual numeric phi values in for the list of solutions.