SanzharMukatay

15 Reputation

3 Badges

5 years, 87 days

MaplePrimes Activity


These are questions asked by SanzharMukatay

Need help with write calculations in maple and ploting

Can anyone help with ploting the movement of the body along the segment AB, BC this task?

restart;
alpha:=Pi/6: V[A]:=1: l:=3: f:=0.2: d:=2.5:
# The movement of the body along the segment AB
de:=m*diff(x1(t),t,t)=G*sin(alpha)-F[fr]:
G:=m*g: F[fr]:=m*g*f*cos(alpha): g:=9.81:
ics:=x1(0)=0, D(x1)(0)=V[A]:
de:=evalf(de/m);
B:=rhs(de);
sol:=dsolve({de,ics}, x1(t)); # The solution of de
fsolve(rhs(sol)=l);
tau:=evalf[3](%[2]);
V['B']:=B*tau+V[A];
 

diff(diff(x1(t), t), t) = 3.205858157

 

3.205858157

 

x1(t) = (3205858157/2000000000)*t^2+t

 

-1.715094465, 1.091236545

 

1.09

 

4.494385391

(1)

# The movement of the body along the curve BC
sys:=m*diff(x(t),t,t)=0, m*diff(y(t),t,t)=m*g;
ics_sys:=x(0)=0, y(0)=0, D(x)(0)=V['B']*cos(alpha), D(y)(0)=V['B']*sin(alpha):
Sol:=evalf(dsolve({sys, ics_sys})); # The solution of sys
T:=evalf(solve(eval(x(t),Sol)=d));
h:=evalf(eval(eval(y(t),Sol),t=T));
T:=evalf[2](T);
h:=evalf[3](h);   
 

m*(diff(diff(x(t), t), t)) = 0, m*(diff(diff(y(t), t), t)) = 9.81*m

 

{x(t) = 3.892251925*t, y(t) = 4.905000000*t^2+2.247192696*t}

 

.6423016927

 

3.466940605

 

.64

 

3.47

(2)

 


 

Download sol_Maple.mw

Hi, need help with ploting task. Could someone help with this?

Task:

Hi, I have a problem with solution task. Why eq2 doesn't accept value of x1?


restart;
P := 12; M := 6; q := 2; BC := 2; Q := q*BC; CN := (1/2)*BC; NB := CN; AK := 1; KC := AK*evalf(cot((1/6)*Pi)); KN := KC+CN; BK := NB+KN;
                            P := 12
                             M := 6
                             q := 2
                            BC := 2
                             Q := 4
                            CN := 1
                            NB := 1
                            AK := 1
                       KC := 1.732050808
                       KN := 2.732050808
                       BK := 3.732050808

eq1 := x1 = (P*evalf(cos((1/6)*Pi))-Q*sin((1/6)*Pi))/sin((1/6)*Pi);
                    eq1 := x1 = 16.78460970

eq2 := x2 = AK*P+BK*x1+KN*Q+M;
            eq2 := x2 = 3.732050808 x1 + 28.92820323

solve({eq1, eq2}, {x1, x2});
              {x1 = 16.78460970, x2 = 91.56921942}
NULL;
B := vector(2, [(P*evalf(cos((1/6)*Pi))-Q*sin((1/6)*Pi))/sin((1/6)*Pi), AK*P+BK*x1+KN*Q+M]);
        B := [16.78460970, 3.732050808 x1 + 28.92820323]

A := matrix(2, 2, [[1, 0], [0, 1]]);
                               [1  0]
                          A := [    ]
                               [0  1]

C := evalm(1/A);
                               [1  0]
                          C := [    ]
                               [0  1]

x := evalm(`&*`(C, B));
        x := [16.78460970, 3.732050808 x1 + 28.92820323]

 

hello everyone, I had mistake from maple plot.

Can someone help me with animation my plot


restart;
with(plots);

with(plottools);

K := 48;
R1 := 30; R2 := 20; OA := R1+R2;
Vodilo := PLOT3D(cuboid([-3, -3, 0], [OA+4, 4, 2]));

Koleso1 := PLOT3D(cylinder([0, 0, 5], R1, 18), cylinder([0, 0, -2], 1, 16));
Koleso2 := PLOT3D(cylinder([0, 0, 7], R2, 6), cylinder([0, 0, -2], 1, 16));
for i from 0 to K do  end do;

Amplituda := .3;
phi := sin(6.28*i/K)*Amplituda;

phi1 := phi*(R1+R2)/R2;

P1 := rotate(Koleso2, 0, 0, -phi1);


P2 := rotate(Vodilo, 0, 0, -phi);

x := OA*cos(phi); y := OA*sin(phi);

P[i] := display(Koleso1, translate(P1, x, y, 0), P2);


display(seq(P[i], i = 0 .. K), orientation = [55, 101], insequence = true, scaling = constrained, axes = normal);
Error, (in plots:-display) expecting plot structure but received: P[0]

 

Page 1 of 1