SimWe仿真论坛's Archiver

hamit 发表于 2008-6-1 23:58

想画2D的图形

我想画 involute(t)  函数的曲线。

c :=t->(t, sin(t)) :
mintt := 0;
maxtt := 4*Pi;
t0 := 0;
sol := dsolve({diff(l(t), t) = eval(sqrt(add((diff([c(tt)], tt))[i]*(diff([c(tt)], tt))[i], i = 1 .. 2)), tt = t), l(t0) = 0}, numeric, method = dverk78, output = listprocedure, maxstep = 0.1e-1);
lgth := subs(sol, l(t));
tangvct := t->evalf(eval((diff([c(tt)], tt))/sqrt(add((diff([c(tt)], tt))[i]*(diff([c(tt)], tt))[i], i = 1 .. 2)), tt = t));
tangvct1 := t->*op(1, expand(tangvct(t))), op(2, expand(tangvct(t))));
involute :=t->evalf(c(t)+(t0-lgth(t))*tangvct1(t));

involute(Pi);

plot([involute(t), t = 0 .. 4*Pi]);

出现错误,请帮个忙。

页: [1]
 

Powered by Discuz! Archiver 6.1.0  © 2001-2007 Comsenz Inc.