function wf()
clc
%%
solve('a*x^2 + b*x + c','x')
solve('a*x^2 + b*x + c','b')
S=solve('x^3 + y^3  = 1','sin(x) - 11*y = 5','x','y')
S.x
S.y
x1=11/1332*(970237+666*2122349^(1/2))^(1/3)-3025/1332/(970237+666*2122349^(1/2))^(1/3)+5/1332
x2=-11/2664*(970237+666*2122349^(1/2))^(1/3)+3025/2664/(970237+666*2122349^(1/2))^(1/3)+5/1332+11/2*i*3^(1/2)*(1/1332*(970237+666*2122349^(1/2))^(1/3)+275/1332/(970237+666*2122349^(1/2))^(1/3))
x3=-11/2664*(970237+666*2122349^(1/2))^(1/3)+3025/2664/(970237+666*2122349^(1/2))^(1/3)+5/1332-11/2*i*3^(1/2)*(1/1332*(970237+666*2122349^(1/2))^(1/3)+275/1332/(970237+666*2122349^(1/2))^(1/3))
y1=1/1332*(970237+666*2122349^(1/2))^(1/3)-275/1332/(970237+666*2122349^(1/2))^(1/3)-605/1332
y2=-1/2664*(970237+666*2122349^(1/2))^(1/3)+275/2664/(970237+666*2122349^(1/2))^(1/3)-605/1332+1/2*i*3^(1/2)*(1/1332*(970237+666*2122349^(1/2))^(1/3)+275/1332/(970237+666*2122349^(1/2))^(1/3))
y3=-1/2664*(970237+666*2122349^(1/2))^(1/3)+275/2664/(970237+666*2122349^(1/2))^(1/3)-605/1332-1/2*i*3^(1/2)*(1/1332*(970237+666*2122349^(1/2))^(1/3)+275/1332/(970237+666*2122349^(1/2))^(1/3))
a=x1^3+y1^3
b=x1-11*y1
%%
s=solve('y*cot(2*x)=-x','y^2=1.05*z^2-x^2','2*z^2-x^2-16=0','x','y','z');
x=s.x
y=s.y
z=s.z
-y*cot(2*x)
sqrt(1.05*z^2-x^2)
2*z^2-x^2-16
%%
x1=3;
x2=4;
x3=5.3;
y1=2.3;
y2=4.5;
y3=5.6;
z1=2.4;
z2=4.5;
z3=3.2;
A=5.6;
B=6.7;
C=8.9;
s=solve('(x1-x0)^2+(y1-y0)^2+(z1-z0)^2=A','(x2-x0)^2+(y2-y0)^2+(z2-z0)^2=B','(x3-x0)^2+(y3-y0)^2+(z3-z0)^2=C','x0','y0','z0');
s
s.x0
s.y0
s.z0