Hello,

I look for a way to color the different subsystems of my system in the 2D mode.

It should be to have a clearer presentation of my system.

Is somebody known as it is possible to color a subsystem ?

Thank you for help.

Hello,

I study the model of crankshaft in MapleSim.

For the moment, i manage to plot different parameters function of the time. For exemple, the angle of the crank function of the time or the displacement of the slider function of the time.

But, i like to plot the displacement of the slider funtion of the angle of the crank.

Is it possible and how can i do this?

Thank you for your help.

Hello,

I have 2 questions about the resolution of a ODE for the simple pendulum equations.

1) I need to add 'maxfun'=10^7 in the dsolve function. With this option, dsolve works. Can you explain me in a very simple way the principle of this option ? It's just to feel better how this option runs.

2) gamma(t) is a angular parameter.

With

tgammax:=odeplot(sol,[gamma(t),x(t)],0..10,numpoints=200,color=blue,legend="z"):
display(tgammax);

I manage to plot x(gamma) but with gamma in rad.

How can i do to x(gamma) but with gamma in degree? I should do gamma_deg=360/(2*Pi)*gamma(t) but as gamma(t) is defined numerically, this lign doesn't work.

Thank you for helpe

Here some ligns of my code :

Cinit:=D(psi)(0)=0,psi(0)=-0.93,D(x)(0)=0,x(0)=1.13,D(gamma)(0)=0,gamma(0)=0.79;sol:=dsolve({sys,Cinit},numeric,'maxfun'=10^7);
tx:=odeplot(sol,[t,x(t)],0..10,numpoints=200,color=blue,legend="x"):
display(tx); 

tgamma:=odeplot(sol,[t,gamma(t)],0..10,numpoints=200,color=blue,legend="gamma"):
display(tgamma);

tgammax:=odeplot(sol,[gamma(t),x(t)],0..10,numpoints=200,color=blue,legend="z"):
display(tgammax);

Hello,

I look for solving a DAE system i obtain after having determined the equations of a mechanical system with kinematic closed loops.

For that, thanks to the partitioning method, i could transform my DAE in ODE system.

But now, i don't manage to solve my ODE system.

The first issue was the calculation of big matrix with trigonometric functions. With your help in the post "Resolution of a big product of matrix with trigonometric function" (http://www.mapleprimes.com/questions/200012-Rsolution-Of-A-Big-Product-Of-Matrix) , i could calculate the different matrix involved are calculated.

Now, the resolution with my differential system is very long and never finished.

@Carl Love 3670 gives me good advices. He asked me to try to avoid symbolic calculations of the matrix.

Do not use the output option for dsolve.

Then the following procedure evaluates the Matrix AA:

AA:= proc(S,t)
local Cu_inv:= eval(Cu, S(t))^(-1), tCu_inv:= eval(tCu, S(t))^(1);
     eval(Avv - Avu.Cu_inv.Cv + tCv.tCu_inv.Auu.Cu_inv.Cv - tCv.tCu_inv.Auv, S(t))
end proc:

           To use it, invoke AA(Sol, t) where t is an actual numeric value.

Unfortunately, i don't see how i can use this method for the moment? But, i think that there is a step i don't understand. For me, i have to make the calculation of the big matrix AA and QQ before solving my differentiel system because my differentiel system is composed with AA and QQ.

Here you can find the system i try to solve.

and here you can find the maple file without the steps leading to the setting of the equations

calcul_des_matrices_.mw

Hello,

I try to calculate a big matrix product which is used after for the resolution of a differential system.

In the link below, you have the matrix product which causes me trouble.

http://dl.free.fr/p8fP9XqWn 

The matrix are composed with trigonometric functions.

I think that the problem is due to the fact that Maple tries to determine a analytical solution.

But, for my use, i only need a numerical result that i can use after for solving my differential system.

Have you some ideas in order to make this matrix product? 

Thank a lot