With this code:
syms x y z t real
Sa=2*pi*30^2*(1-sin(pi/18))+2*pi*t*(30*cos(pi/18)+(t/2)*tan(pi/18))/cos(pi/18);
Sa = simplify(Sa, 'steps', 10);
Sb=170^2-pi*30^2;
F=-0.3*Sa/(1+x)+((Sa*x^2)/(1+x)+(Sb*y^2)/(1+y))/2+z*(Sa/(1+x)+Sb/(1+y)-(170^2)/1.1);
F = simplify(F, 'steps',10);
[x,y,z,t] = solve((x^2)/2+x+0.3-z==0, (y^2)/2+y-z==0, (1+x)*((170^2)/(1.1)-(Sb)/(1+y))==Sa, F==0, [x, y, z, t])
and after several minutes (I didn’t time it), I got:
x =
-1.0
y =
-0.22540333075851662296414692004352
z =
-0.2
t =
-26.959447166136052326594445120517
Best Answer