Calcution-A

with(VectorCalculus):
with(LinearAlgebra):
f := [x3*(1+x2),x1,x2*(1+x1)]:
g := [0,1+x2,-x3]:
fv := <x3*(1+x2),x1,x2*(1+x1)>:
gv := <0,1+x2,-x3>:
ad_f_g :=  MatrixAdd(MatrixVectorMultiply(Jacobian(g,[x1,x2,x3]),fv),-MatrixVectorMultiply(Jacobian(f,[x1,x2,x3]),gv)):
b := <Diff(w(x1,x2,x3),x1),Diff(w(x1,x2,x3),x2),Diff(w(x1,x2,x3),x3)>:
gv := <0,1+x2,-x3>:
c1 := DotProduct(b, gv):
c2 := DotProduct(b, ad_f_g):

Is there anyway to disable the superscript so it works exactly like in Classic Maple, so the whole equation stays inline? Having it go to where it should be makes marginally easier for me to read it, but having to press the right arrow to leave the superscript is extremely annoying and time consuming.

Edit: example: I want it to stay as 10^5 instead of 10⁵, at least while I'm writing.

Hello,

I'm new to Maple and am using it to solve linear and non-linear differential equations.

I can solve and plot linear DEs using a couple of methods. Im trying to solve a non-linear DE and i could find the series solution using dsolve. But given certain initial conditions, is there any way to plot or see the nature of the behaviour of the DE? I just want to see how the DE will look like maybe for certain points up to a certain level of accuracy. Is there any way to do that?

i have 2 matrices :    |1 2 3|     |1 2 3|
                               |4 5 6|     |4 5 6|
i have to multiply there 2 matrices like below

a= (1*6)+(1*5)+(1*4)+(1*3)+(1*2)+(1*1)=21
b= (2*6)+(2*5)+(2*4)+(2*3)+(2*2)+(2*1)=42
c= (3*6)+(3*5)+(3*4)+(3*3)+(3*2)+(3*1)=63

I need this quite often, but never found a nice way of getting the index of minimal or maximal value of a Vector (or a list or whatever). Is there an easier way than the following?

N := 15:
L := LinearAlgebra[RandomVector](N):
m := min(L):
for i to N do if L[i] = m then print(i); break end if end do;

I have a data set that represents some parameters at different locations in a (linear) system. The way the data are organized is as 1-d vectors which in turn are combined into an Array. I.e. Array[1] is the first vector; Array[2] the second and so on. Each vector represents a data subset, all subsets are against a common independent variable, call it Knob. I plot one of these subsets like this:

plot(knob,vector[1],style=point);

That works fine & I can replace...

Given Maple expressions f(g,x) and g(x) how would you replace g in f  to show a function in x.

c:=12+10x^4+7x^3-89x^2+65x+4sinx

how do i plot a graph of c,

calculate the roots and check the correctness of all found roots.

Thanks for your kind solutions in advance.

i have this unfinished procedure to multiply two matrices 

multmatrix:= proc(A:matrix)

local m,n,i,j,r,c,t,a;

m:=linalg[rowdim](A);

n:= linalg[coldim](A);

i:=linalg[rowdim](B);

j:=linalg[coldim](B);

K:=array(1..m,1..n);

H:=array(1..i,1..j);

I need help with my design project. My project is about purification of americium from stock pile plutonium. I need you help in desgning a double pipe horizontal condenser with a special E-shell (3nozzles): 1:1 countercurrent.

Waste gases from two furnaces are run through a condenser to reduce the amount of water that enters the NOx abatement stage. The non-condensables and a small amount of water that remains in the vapour are in the vapour are then passed onto the...

Hi all,

During simulation I need to use a parameter as a feedback, which is previously determined in a Modelica Custom Component. In this case, can we use Zero Order Hold (ZOH) block/function to hold the previous value? If yes, how should we choose the parameters of ZOH?

Thanks

Hi,

I'm developing an application in Maple 14. I'm using the code below to print a percentage:

# Printing a percentage.
x:=0.7569:
printf("The value of x is %.1f%%%", x*100);

Is this the right way to print a percentage?

I also want to group digits (as in 1,959), but don't know how to do it.

Thanks in advance,

Carlos Mallen

Hi, these are what I have, I'm trying tot get a set of answers.

elenum:=16:

nodes:=elenum*2+1:

L:=evalf((Pi*2)/(elenum*2)):

interface(rtablesize=infinity):

E0:=Matrix([[L/3,2*L/3,L/3)],[(2*L/3,11*L/15,2*L/3)],[(L/11,2*L/3,L/110]]):

for i from 1 to elenum do
node:=(2*(i-1))+1;
E0M[node..node+2,node..node+2]:=E0M[node..node+2,node..node+2]+E0
end do:

E0M:

What should I do if i want a different outcome of...

Hey there,

i´ve a numerical solved differential equation and a "normal" equation. How can i calculate the differenz of both?

>S:=0.79*cos(t)+0.3*sin(t);

> DGL1 := 20*(D(D(x)))(t)+10*(D(x))(t)+2500*x(t)+300 = S;
> init := x(1.2) = 0, (D(x))(1.2) = 10.63081252;
> F := dsolve({DGL1, init}, numeric);
>U:=F-S;

> p1 := plot(U, [t, x(t)], 1.2 .. 2, numpoints = 200);
%;
Error, (in plot) two lists...