倪举鹏 发表于 2014-10-12 14:46:48

会软件的画个积分图

某质点初始位置(0,0),初速度0,受到的加速度随时间t变化的矢量为cos(Pi/3*cos(38/41*t))*[-cos(Pi/3*cos(38/41*t)), -sin(Pi/3*cos(38/41*t))]    请画质点的移动轨迹图

BeerRabbit 发表于 2014-10-13 16:44:12

倪举鹏 发表于 2014-10-13 18:16:41

x := int(-sin(theta)*sin(-arccos((1/6)*sqrt(1-cos(theta)^2/(1+2*cos(theta))^2)*(2+cos(theta))*sqrt(6)/(sqrt(1/(1+2*cos(theta)))*(1+cos(theta))))+int((1/36)*sqrt(6)*(2*sqrt(6)*sqrt((9*cos(theta)^3+15*cos(theta)^2+10*cos(theta)+2)/(2*cos(theta)^2+3*cos(theta)+1))*sqrt(3*cos(theta)^2+4*cos(theta)+1)*cos(theta)+sqrt(6)*sqrt((9*cos(theta)^3+15*cos(theta)^2+10*cos(theta)+2)/(2*cos(theta)^2+3*cos(theta)+1))*sqrt(3*cos(theta)^2+4*cos(theta)+1)-sin(theta)*sqrt(-(18*cos(theta)^3+6*cos(theta)^2-12*cos(theta)-12)/(2*cos(theta)^2+3*cos(theta)+1)))/(sqrt(3*cos(theta)^2+4*cos(theta)+1)*sqrt(1/(1+2*cos(theta)))*(2*cos(theta)^2+3*cos(theta)+1)), theta))/(sqrt(3*cos(theta)^2+4*cos(theta)+1)*(1+2*cos(theta))), theta)

y := int(sin(theta)*cos(-arccos((1/6)*sqrt(1-cos(theta)^2/(1+2*cos(theta))^2)*(2+cos(theta))*sqrt(6)/(sqrt(1/(1+2*cos(theta)))*(1+cos(theta))))+int((1/36)*sqrt(6)*(2*sqrt(6)*sqrt((9*cos(theta)^3+15*cos(theta)^2+10*cos(theta)+2)/(2*cos(theta)^2+3*cos(theta)+1))*sqrt(3*cos(theta)^2+4*cos(theta)+1)*cos(theta)+sqrt(6)*sqrt((9*cos(theta)^3+15*cos(theta)^2+10*cos(theta)+2)/(2*cos(theta)^2+3*cos(theta)+1))*sqrt(3*cos(theta)^2+4*cos(theta)+1)-sin(theta)*sqrt(-(18*cos(theta)^3+6*cos(theta)^2-12*cos(theta)-12)/(2*cos(theta)^2+3*cos(theta)+1)))/(sqrt(3*cos(theta)^2+4*cos(theta)+1)*sqrt(1/(1+2*cos(theta)))*(2*cos(theta)^2+3*cos(theta)+1)), theta))/(sqrt(3*cos(theta)^2+4*cos(theta)+1)*(1+2*cos(theta))), theta)

BeerRabbit 发表于 2014-10-13 19:07:15

倪举鹏 发表于 2014-10-13 18:16
x := int(-sin(theta)*sin(-arccos((1/6)*sqrt(1-cos(theta)^2/(1+2*cos(theta))^2)*(2+cos(theta))*sqrt(6 ...

我用的是MMA,没装Maple呢…………
页: [1]
查看完整版本: 会软件的画个积分图