- 积分
- 1
- 注册时间
- 2008-2-27
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2015-1-30 19:31:47
|
显示全部楼层
来自 陕西西安
本帖最后由 高贵的点点滴滴 于 2015-1-30 19:44 编辑
- Clear["Global`*"]
- v = 0;
- r = 1;
- d = 10;
- (*f[t0_?NumericQ,t1_?NumericQ]:=(1+v)Log[r]+Log[t0]+(1+v)Integrate[ \
- Sin[s],{s,t0,t1}];
- g[t0_?NumericQ,t1_?NumericQ]:=(t1/t0)^(1/(1+v))Sqrt[(1-t0)/(1-t1)]E^-\
- Integrate[1,{s,t1,t0}]-d;*)
- f[x_, y_] := Sin[x y];
- g[x_, y_] := Cos[x^2 y];
- fig = ContourPlot[{f[x, y] == 0, g[x, y] == 0}, {x, 0, Pi}, {y, 0,
- Pi}, PlotPoints -> 10, Contours -> 50];
- data = fig[[1, 1]];
- refdata1 = Select[data, Abs[f @@ #] + Abs[g @@ #] < 10^-1 &];
- refdata2 = Union[refdata1, SameTest -> (Norm[#1 - #2] < 10^-3 &)];
- pts = {x,
- y} /. (FindRoot[{f[x, y] == 0,
- g[x, y] == 0}, {{x, #[[1]]}, {y, #[[2]]}},
- AccuracyGoal -> 4] & /@ refdata2)
- ContourPlot[{f[x, y] == 0, g[x, y] == 0}, {x, 0, Pi}, {y, 0, Pi},
- Epilog -> {AbsolutePointSize[5], Red, Point /@ pts}]
复制代码
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|