vpasolve用法,vpasolve函数怎么用
求助，引用方程组的符号解
符号方程求解：solve求解代数方程，方程组用linsolve, 数值解用vpasolve,格式[x1,s3,x3,x4]=vpasolve(eqns,vars,int_guess),分别是方程组，变量，初始值，具体百度一下用法就知道啦
如图，解方程组！！！
解：这是一个轮换对称函数，因此，x=y=z;
所以：x^4+x^2-5x^2=x^2-4x+4=(x^2-2)^2=[(x+√2)(x-√2)]^2=0;
x=y=z=+/-√2/2; 是方程组的解。
matlab使用for循环绘制曲线图时，如何自定义曲线颜色？
matlab使用for循环绘制曲线图，自定义每个曲线图的颜色，可以使用plot(x,y,'--','color',rand(1,3));这样的代码实现，因为代码的外层for循环只跑4次，画4次plot，在开头加一个clf擦除旧图就可。
效果图如下
代码修改列出如下：
clc;clear all;
clf;
syms rp r A0 B0 py a pw gama p0
b delta a0 E phi;
c=2.02E3;
b=0.466;
E=2.4E6 ;
a=0.0015;
a0=0;
delta=0;
gama=1025*9.8;
p0=0;
phi=0.00000000001:0.5:5.00000000001;
figure(1);
hold on;
for d=0.03:0.03:0.12;
n=zeros (size (phi));
m=zeros (size (phi));
for i=1: length (phi);
A0= (1+sin(phi(i)*pi/180))/ (1-sin(phi(i)*pi/180));
B0= (2*c*cos (phi(i) *pi/180))/ (1-sin(phi(i)*pi/180));
pw =-(6*rp.^2*(2*A0 - 1)*(B0 - p0 + A0*p0))/ (A0*a^3* (rp/a) ^ (2/A0)*(A0 + 2)*((2*A0)/(a*(A0 - 2)) - 1/(a*(A0 - 2)) + (8*(A0 - 2*b)*(rp/a)*((40 - 2)/A0)*(A0 - 1))/(a*(b - 1)*(A0^2 - 4)) - (4*(A0 - 2*b)*(rp/a)^ ((40 - 2)/A0)* (A0 - 1))/(A0*a*(b - 1)*(A0^2 - 4)) + 1));
eqn=((18.*(A0.*b+b-1).*p0 +6.*B0.*(b+1))./ (2+A0)).*rp.^3+ 6.*(1+b).*a.*pw.*((A0.*b+b-1)./((2+A0).*
(b-1))).*rp.^2 -2.*E.*a.^3==0;
m=vpasolve(eqn, rp,[0, 3]);
n(i)= (m-a)/(d/2);
end;
plot(phi, 100*n,'--','color',rand(1,3));
end
;
axis([O 5 10 100]);
hold off;
代码局部截图如丐
f(x)=coth(x)-1/x如何用Matlab求出x
此方程可以直接用vpasolve（）函数命令求出。
syms x
x=vpasolve(coth(x)-1/x)
x =-0.00000000000000000000000000000011308787783892655728288492802931
验算：
coth(x)-1/x
ans =0.0
说明：vpasolve（）函数适用高版本，对于低版本，可用solve（）函数来计算，用法同上。