会软件的画个积分图

倪举鹏

某质点初始位置（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

倪举鹏

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 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呢…………