作业帮 > 数学 > 作业

mathematica中解方程组并将结果画出图来

来源:学生作业帮 编辑:作业帮 分类:数学作业 时间:2024/05/08 18:37:45
mathematica中解方程组并将结果画出图来
代码如下
Q = 0
\[Mu] = 0
sol = Solve[
y[a,k,mu]*(a^2 *(1 + r) + r^2 *(-3 + 2 *Q^2/r + r)) ==
a *mu* (1 - r) +
k* Q* (a^2 + Q^2 - r) && (y[a,k,mu]*(r^2 + (a)^2) - mu*(a) -
k*Q*r)^2 == (r^2 - 2*r + a^2 +
Q^2)*(1 + ((a)^2)/2*(1 + mu^2)*y[a,k,mu]^2 -
2*(a)*mu*y[a,k,mu] + (\[Mu]^2)*r^2) &&
r \[GreaterSlantEqual] (1 + (1 - (a^2 + Q^2))^(1/2)),{y[a,k,
mu],r},Reals];
T1 = ParametricPlot[Evaluate[y[1,0,mu] /.sol],{mu,-1,1},
PlotRange -> {0.1,0.5},Frame -> True,AspectRatio -> 1/1]
其中,y=y[a,k,mu].求解出y后画出y随mu变化的图像 a可以取若干个数值,k也可以去若干个数值.但是我的不出图,求大神给看看.
可以出的,


再问: 如果我想依次 画出y[1,0,mu] y[0.99,0,mu] y[0.9,0,mu] y[1,1,mu] y[1,2,mu]怎么弄?另外,要是把您画的这个的直线部分的下部分抹掉怎么弄?如果分开弄的话,我知道怎么弄。另外,您画的这个图中直线部分是r=1的时候,亦即方程组中r>=(1 + (1 - (a^2 + Q^2))^(1/2))中的等号部分。求教。
再答: 不好意思,关于直线部分其实也是一个解,ta在列表首First可以获取,Plot[Evaluate[y[1, 0, mu] /. Take[sol,{2,Length@sol}]可以不画这个解;关于你的Solve求解很慢,我现在没想到好办法把它做成函数。