Dear Users,

I would like to know how to convert the following bessel product into a polynomial form.

I have an experession as follows:

BesselI(1,p*x)*BesselK(0,p*x)+BesselI(0,p*x)*BesselK(1,p*x)    #where p is arbitrary real positive constant and x is a variable

How to get the result in maple by simplify, convert or any other functoins.

Thanks for your attention.

i'am using maplet shortcut in maple 

but it did not show in table form, it only gives final answer.

use DocumentTools in 
with(LinearAlgebra):with(VectorCalculus):with(linalg):

# F1,F2,F3 are my equations 
F1:=Do(%f1); 
F2:=Do(%f2);
F3:=Do(%f3);

# these is my initial point
X1:=Do(%x1);
X2:=Do(%x2);
X3:=Do(%x3);

# the whole operation
F := Matrix(1 .. 3, 1 .. 1, [[F1], [F2], [F3]]);
J := Jacobian([F1, F2, F3], [x1, x2, x3]);

I have the following expression 

lambda1-2.9881355172134833*lambda2+16.044692204509636*lambda1^2-114.31727386075951*lambda1*lambda2+202.36051939112247*lambda2^2

I want to get the expression in a form in which all the coeffs are fixed to 7 decimal places. 

I have tried with 

evalf[7](lambda1-2.9881355172134833*lambda2+16.044692204509636*lambda1^2-114.31727386075951*lambda1*lambda2+202.36051939112247*lambda2^2)

Hi guys. I need to do a Runge-Kutta method to solve two equations and I'm having trouble at using functions. My program is as shown below:

restart; A := Matrix();

dt := 0.1e-2;

diff(x(t), t) := v;

diff(v, t) := a(x(t), v);

xo := 0; vo := 0;

for i to 1000 do

xn := xo+v.dt;

vn := vo+a.dt;

dx1 := vn.dt; dv1 := a.dt;

dx2 := (vn+(1/2)*dv1)*dt; dv2 := a(x(t)+(1/2)*dx1, (1/2)*dv1)*dt;

dx3 := (vn+(1/2...

of the integral int(cos(t^3)/(t+2), t = -infinity .. x) as x tends to -infinity with Maple? I have in mind a better asymptotics than 0.

can someone assist me to animate the following:

 A ladder is leaning against a wall. The base of the ladder starts sliding away from the wall, with the top of the ladder sliding down the wall. As the ladder slides, you watch the red point in the middle of the ladder. What figure does the red point trace? (a circle)

[24+24*beta[1]*x+24*beta[2]*x^2+24*beta[3]*x^3+24*beta[4]*x^4]+(2*[[beta[1]+2*beta[2]*x+3*beta[3]*x^2+4*beta[4]*x^3]])*(-24+24*x)+[2*beta[2]+6*beta[3]*x+12*beta[4]*x^2]*(12-24*x+12*x^2)-L^4*m[0]*Omega^2*[-(1/2)*x^2-R*x+R+1/2]*(12-24*x+12*x^2)/EI[0]-F*L^2*(12-24*x+12*x^2)/EI[0]-L^4*m[0]*Omega^2*[-x-R]*(12*x-12*x^2+4*x^3)/EI[0]-L^4*m[0]*omega^2*(6*x^2-4*x^3+x^4)/EI[0]

This is the result of a mathematical procedure. I need to be able to extract the series...

I have already solved my ODEs and used odeplot to analyse the change over time (0 to 10).

Below is my ODE,

d(a(t))/dt = x - y*(a(t))

d(b(t))/dt = m*a - n*(b(t))

My query here is,

I need to substitute a(t) and b(t) in a seperate equation,

h(t) = 60/0.98 * (a(t))*(b(t))

and i need to plot 'h' over time (0 to 10)

How do I go about doing this?

Please help.

Thanks.

Hi,

i am trying to simulate a costum component, but i get this error when starting to compile initial values i get this error:

(in dsolve/numeric/RunSimulation) keyword event_blindeye not recognized
Simulation failed.

Why? This costum component was tested on another PC and worked.

Thanks

Bear with, I'm a noob. This is my code so far:

    > ball := proc (x, y) plots[pointplot]([[x, y]], color = blue, symbol = solidcircle, symbolsize = 20) end proc;

    > animate(ball, [1+sin(t), 1-sin(t)], t = 0 .. 2*Pi, scaling = constrained, frames = 100);

I like the output I have but my aim is to have, say, 8 of these balls, all with their own trajectories, and animated on the same axis. All the trajectories would depend only on one parameter, t.

Hello,

could you pls help me with this code, i dont understand why i get the error.

---------------------------------------------------------------------------------------------

NewtonZweiVariablen := proc (funktion1, funktion2, startwertX, startwertY)

local Ableitung, StartwertNeuVektor, DurchgangsZaehler, Differenz, StartwertVektor, FunktionsVektor, Jacobi, g, h;

Jacobi := linalg:-inverse(Jacobian([funktion1, funktion2], [x, y]));

Hi!

How to solve attached equation in Maple?

Is there a possiblity to solve using units?

wzel

PS

Solution was done in Mathcad

Hello Community,

Wondering if anyone has run into the problem outlined in my code below. I am looking to evalute a commutator involving Pauli spin matrices. One of the matrices is squared. Such a commutator should evaluate to zero. However, the way I implemented it below doesn't give such a simplified result. See Result 7 below. Any suggestions as to how to get Maple to evaluate the exponential before the Commutator? Or to otherwise get to a final, simplified result?

Hi,

I am trying to solve the following system using the function "solve", but the message "warning solutions may have been lost" appears.

This is my script (a,b,c,d,e,f, Resp, Malcyt,Pyrcyt are parameters). The unknown variables are Citcyt, Citmt,Malmt,Pyrmt:

e1 := 0 = a.Pyrcyt-a.Pyrmt+(b-c).Malmt;
 e2 := 0 = d.Malcyt+e.Citmt-(d+b+c).Malmt;
 e3 := 0 = f.Citcyt+c.Malmt-(e+f).Citmt;

 e4 := 0 = 2.*e.Citmt+(b+c).Malmt-Resp;

Hello every one,

I have a question about summing data from an ode.

ode:=diff(T(x),x$2)-T(x)^2=0;

bc:=T(0)=1,T(n)=0;

n:=3:

Thanks