Advanced Search

Maple 2015 Questions and Posts

These are Posts and Questions associated with the product, Maple 2015
View unanswered posts
Maple 2015 Questions and Posts Feed

print text and Equations...

Asked by: minhhieuh2003 140
2dmath text inline
January 23 2018
1  5

How do you way text below

Maple 2015 

if i use function cat() is not 

 

Sorting Many Tables...

Asked by: Les 190
worksheet print sort
January 22 2018
0  6

The attached worksheet performs two functions:

(1) It lets me print 4 × 6 Index Cards for the short entries in each table.

(2) It allows for easy storage and retrieval of syntax (code).

The worksheet has many tables, each separated by a Page Break.

Questions:

(a) Is there a way to sort all the different tables so they will be arranged in alphebetical order?

(b) When I select one table to print and open the Print Dialog, the "Selection" option is grayed out. (see graphic below).  (1) Is there a way to enable the selection option?  (2) Is there a way to determine what page I am on so I can use the "Pages from...to" option?  If I need to number the pages, will the page numbers reset to parallel a new alphabetic sort order.

Many thanks in advance.  See WC29_4_BY_6_NOTE_CARDS_UNSORTED.mw attached. And see image of Print Dialog below.

Les

Need the numerical solution of system of ODEs...

Asked by: Muhammad Usman 245
ode dsolve nonlinear
January 21 2018
0  4

Dears 

Hope you all are fine. I want to solve the following nonlinear system ODEs numerically

Eq1 := diff(F(eta), eta, eta, eta)-phi1*((diff(F(eta), eta))^2-(diff(F(eta), eta)+diff(G(eta), eta))*(diff(F(eta), eta, eta)))+alpha*(1-phi)^2.5*((diff(F(eta), eta, eta))^2+2*(diff(F(eta), eta))*(diff(F(eta), eta, eta))-(diff(F(eta), eta)+diff(G(eta), eta))*(diff(F(eta), eta, eta, eta, eta))+A*(2*(diff(F(eta), eta, eta, eta))+(1/2)*eta*(diff(F(eta), eta, eta, eta, eta))))-A*phi1*(diff(F(eta), eta)+(1/2)*eta*(diff(F(eta), eta, eta))) = 0; Eq2 := diff(G(eta), eta, eta, eta)-phi1*((diff(G(eta), eta))^2-(diff(F(eta), eta)+diff(G(eta), eta))*(diff(G(eta), eta, eta)))+alpha*(1-phi)^2.5*((diff(G(eta), eta, eta))^2+2*(diff(G(eta), eta))*(diff(G(eta), eta, eta))-(diff(F(eta), eta)+diff(G(eta), eta))*(diff(G(eta), eta, eta, eta, eta))+A*(2*(diff(G(eta), eta, eta, eta))+(1/2)*eta*(diff(G(eta), eta, eta, eta, eta))))-A*phi1*(diff(F(eta), eta)+(1/2)*eta*(diff(F(eta), eta, eta))) = 0;

assoicated with the following BCs

F(0)=0,D(F)(0)=1,G(0)=0,D(G)(0)=p, D(F)(L)=0, D( G(L)=0;

for phi1 := 1.2; alpha := 2; A := 1.5;phi:=0.1;

Please help me to find the solution

Special request to @acer , @Carl Love @Preben Alsholm

conclude not equations...

Asked by: minhhieuh2003 140
solve
January 20 2018
1  9

can result ?

Please help me?

Does anyone can explain the evaluation rules for "...

Asked by: mmcdara 7896
evalf
January 20 2018
0  14

Hi,

I am a little bit surprised by the result of the operation evalf[8](f(y)) in the piece of code that follows.
I was expected the answer to be 2.4494897, not 2.4494898.

Happily the sequence
res := f(y) ; evalf[8](res)
returns the expected result 2.4494897

I suspect the difference comes from some precedence of the operators (f and evalf) but I can't figure out what really happens

Could you enlight me please ?

Thanks in advance

 

> 

restart:
> 

interface(version);

`Standard Worksheet Interface, Maple 2015.2, Mac OS X, December 21 2015 Build ID 1097895`

 (1)
> 

Digits;

10

 (2)
> 

f := x -> sqrt(2.0)*x;

proc (x) options operator, arrow; sqrt(2.0)*x end proc

 (3)
> 

y := sqrt(3.0):
> 

f(y);

2.449489743

 (4)
> 

evalf[9](f(y));  # right

2.44948974

 (5)
> 

evalf[8](f(y));  # ????

2.4494898

 (6)
> 

res := f(y);
evalf[8](res);  # right

2.449489743

  

2.4494897

 (7)
> 

 

 

 

Help limit fraction...

Asked by: minhhieuh2003 140
January 20 2018
1  2

with maple

How can in maple 2015? Help me?

problem plotting double integration expressions...

Asked by: Giulianot 25
plot integration numeric
January 19 2018
0  6

hi everyone,

i have attached a maple worksheet which you can see the issue...azido_displacement.mw
i think tittle says by itself... thanks in advance for taking the time to review and aswer me.

 

theta__o := (1/4)*Pi

(1/4)*Pi

 (1)

omega__o := 0

0

 (2)

tau := 1

1

 (3)

m := 2.28335

2.28335

 (4)

g := 9.80665

9.80665

 (5)

L := .35

.35

 (6)

Iota := 0.9996799726e-1

0.9996799726e-1

 (7)

with(DirectSearch)

[BoundedObjective, CompromiseProgramming, DataFit, ExponentialWeightedSum, GlobalOptima, GlobalSearch, Minimax, ModifiedTchebycheff, Search, SolveEquations, WeightedProduct, WeightedSum]

 (8)

SolveEquations([omega__o+(1/3)*alpha__1*tau-(1/3)*alpha__2*tau+(1/3)*alpha__3*tau = 0, theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(1/3)*(omega__o+(1/3)*alpha__1*tau)*tau-(1/2)*alpha__2*((1/3)*tau)^2+(1/3)*(omega__o+(1/3)*alpha__1*tau-(1/3)*alpha__2*tau)*tau+(1/2)*alpha__3*((1/3)*tau)^2 = (1/2)*Pi, int((m*g*cos(omega__o*t+theta__o+(1/2)*alpha__1*t^2)+alpha__1*(L*m+Iota))/(m*sin(omega__o*t+theta__o+(1/2)*alpha__1*t^2)), t = 0 .. (1/3)*tau)+int((m*g*cos(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau)*t-(1/2)*alpha__2*t^2)-alpha__2*(L*m+Iota))/(m*sin(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau)*t-(1/2)*alpha__2*t^2)), t = 0 .. (1/3)*tau)+int((m*g*cos(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(1/3)*(omega__o+(1/3)*alpha__1*tau)*tau-(1/2)*alpha__2*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau-(1/3)*alpha__2*tau)*t+(1/2)*alpha__3*t^2)+alpha__3*(L*m+Iota))/(m*sin(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(1/3)*(omega__o+(1/3)*alpha__1*tau)*tau-(1/2)*alpha__2*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau-(1/3)*alpha__2*tau)*t+(1/2)*alpha__3*t^2)), t = 0 .. (1/3)*tau) = 0], initialpoint = [alpha__1 = 12.7751705882228, alpha__2 = 18.4817577058678, alpha__3 = 5.70658711764534])

[6.74607137501932*10^(-24), Vector(3, {(1) = HFloat(1.936228954946273e-13), (2) = HFloat(2.027322754116767e-12), (3) = 0.1612e-11}), [`#msub(mi("α",fontstyle = "normal"),mi("1"))` = 14.7208062595154, `#msub(mi("α",fontstyle = "normal"),mi("2"))` = 22.3730290484357, `#msub(mi("α",fontstyle = "normal"),mi("3"))` = 7.65222278892092], 139]

 (9)

alpha__1 := 14.7208062595154

14.7208062595154

 (10)

alpha__2 := 22.3730290484357

22.3730290484357

 (11)

alpha__3 := 7.65222278892092

7.65222278892092

 (12)

x__1 := int(int((m*g*cos(omega__o*t+theta__o+(1/2)*alpha__1*t^2)+alpha__1*(L*m+Iota))/(m*sin(omega__o*t+theta__o+(1/2)*alpha__1*t^2)), t = 0 .. t), t = 0 .. t2)

int(int(.4379530076*(22.39201428*cos((1/4)*Pi+7.360403130*t^2)+13.23607306)/sin((1/4)*Pi+7.360403130*t^2), t = 0 .. t), t = 0 .. t2)

 (13)

x__2 := int(int((m*g*cos(omega__o*t+theta__o+(1/2)*alpha__1*t^2)+alpha__1*(L*m+Iota))/(m*sin(omega__o*t+theta__o+(1/2)*alpha__1*t^2)), t = 0 .. t), t = 0 .. (1/3)*tau)+t2*(int((m*g*cos(omega__o*t+theta__o+(1/2)*alpha__1*t^2)+alpha__1*(L*m+Iota))/(m*sin(omega__o*t+theta__o+(1/2)*alpha__1*t^2)), t = 0 .. (1/3)*tau))+int(int((m*g*cos(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau)*t-(1/2)*alpha__2*t^2)-alpha__2*(L*m+Iota))/(m*sin(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau)*t-(1/2)*alpha__2*t^2)), t = 0 .. t), t = 0 .. t2)

int(int(.4379530076*(22.39201428*cos((1/4)*Pi+7.360403130*t^2)+13.23607306)/sin((1/4)*Pi+7.360403130*t^2), t = 0 .. t), t = 0 .. 1/3)+t2*(int(.4379530076*(22.39201428*cos((1/4)*Pi+7.360403130*t^2)+13.23607306)/sin((1/4)*Pi+7.360403130*t^2), t = 0 .. 1/3))+int(int(-.4379530076*(22.39201428*cos(-1.603220734-4.906935420*t+11.18651452*t^2)-20.11649647)/sin(-1.603220734-4.906935420*t+11.18651452*t^2), t = 0 .. t), t = 0 .. t2)

 (14)

x__3 := int(int((m*g*cos(omega__o*t+theta__o+(1/2)*alpha__1*t^2)+alpha__1*(L*m+Iota))/(m*sin(omega__o*t+theta__o+(1/2)*alpha__1*t^2)), t = 0 .. t), t = 0 .. (1/3)*tau)+(1/3)*tau*(int((m*g*cos(omega__o*t+theta__o+(1/2)*alpha__1*t^2)+alpha__1*(L*m+Iota))/(m*sin(omega__o*t+theta__o+(1/2)*alpha__1*t^2)), t = 0 .. (1/3)*tau))+int(int((m*g*cos(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau)*t-(1/2)*alpha__2*t^2)-alpha__2*(L*m+Iota))/(m*sin(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau)*t-(1/2)*alpha__2*t^2)), t = 0 .. t), t = 0 .. (1/3)*tau)+t2*(int((m*g*cos(omega__o*t+theta__o+(1/2)*alpha__1*t^2)+alpha__1*(L*m+Iota))/(m*sin(omega__o*t+theta__o+(1/2)*alpha__1*t^2)), t = 0 .. (1/3)*tau)+int((m*g*cos(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau)*t-(1/2)*alpha__2*t^2)-alpha__2*(L*m+Iota))/(m*sin(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau)*t-(1/2)*alpha__2*t^2)), t = 0 .. (1/3)*tau))+int(int((m*g*cos(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(1/3)*(omega__o+(1/3)*alpha__1*tau)*tau-(1/2)*alpha__2*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau-(1/3)*alpha__2*tau)*t+(1/2)*alpha__3*t^2)+alpha__3*(L*m+Iota))/(m*sin(theta__o+(1/3)*omega__o*tau+(1/2)*alpha__1*((1/3)*tau)^2+(1/3)*(omega__o+(1/3)*alpha__1*tau)*tau-(1/2)*alpha__2*((1/3)*tau)^2+(omega__o+(1/3)*alpha__1*tau-(1/3)*alpha__2*tau)*t+(1/2)*alpha__3*t^2)), t = 0 .. t), t = 0 .. t2)

int(int(.4379530076*(22.39201428*cos((1/4)*Pi+7.360403130*t^2)+13.23607306)/sin((1/4)*Pi+7.360403130*t^2), t = 0 .. t), t = 0 .. 1/3)+(1/3)*(int(.4379530076*(22.39201428*cos((1/4)*Pi+7.360403130*t^2)+13.23607306)/sin((1/4)*Pi+7.360403130*t^2), t = 0 .. 1/3))+int(int(-.4379530076*(22.39201428*cos(-1.603220734-4.906935420*t+11.18651452*t^2)-20.11649647)/sin(-1.603220734-4.906935420*t+11.18651452*t^2), t = 0 .. t), t = 0 .. 1/3)+t2*(int(.4379530076*(22.39201428*cos((1/4)*Pi+7.360403130*t^2)+13.23607306)/sin((1/4)*Pi+7.360403130*t^2), t = 0 .. 1/3)+int(-.4379530076*(22.39201428*cos(-1.603220734-4.906935420*t+11.18651452*t^2)-20.11649647)/sin(-1.603220734-4.906935420*t+11.18651452*t^2), t = 0 .. 1/3))+int(int(.4379530076*(22.39201428*cos(1.995919816-2.550740930*t+3.826111394*t^2)+6.880423404)/sin(1.995919816-2.550740930*t+3.826111394*t^2), t = 0 .. t), t = 0 .. t2)

 (15)

plot([x__1, x__2, x__3], t2 = 0 .. (1/3)*tau)

  

``


 

Download azido_displacement.mw

 

Help with polarplot and complexplot ...

Asked by: Les 190
complexplot polarplot plotting
January 14 2018
1  2

Could someone help me with the following.  The syntax produces an unfinished graph with a warning.

> with(plots);
> z := polar(1.05, (1/10)*Pi);
                             
> display(polarplot(1, color = grey, axis[radial] = [color = "Blue"]), complexplot(seq(evalc(z))^n,
n = 1 .. 21));

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

the graphic looks like this

The graphic should look like this

Thanks, any help appreciated.

Les

 

Animation of the numerical solution of ODE...

Asked by: Muhammad Usman 245
ode animation
January 09 2018
0  2

Hellow dears!!!

Hope everyone is fine with everything. I want the animation of the numerical solution of ODE i,e., f4 for delta=[0,1.5]. Please see the attachment and fix my problem. I Shall be very thankful to you.

Graph.mw

Special request to  acer 13834, @Carl Love ,Preben Alsholm 10271

Mistake on Maplesoft 2015, Tools -> Math Apps -...

by: Les 190 Maple 2015
January 05 2018
1  3

Irrational numbers: numbers that cannot be represented as a ratio of integers. The decimal form of a rational number is non-repeating and non-terminating.

Change to:

Irrational numbers: numbers that cannot be represented as a ratio of integers. The decimal form of an irrational number is non-repeating and non-terminating. 

or change to

Irrational numbers can be represented by decimal fractions in which the digits go on forever without ever repeating a pattern.  See Downing, Douglas. Dictionary of Mathematics Terms. 2nd ed. Hauppauge, NY: Barron's Ed. Series, Inc., 1995, p. 176).

how to calculate a reduced groebner basis ...

Asked by: Annonymouse 160
December 25 2017
0  2

It is a truth universally acknowledged, that a single man in possession of an algorithm for calculating a Groebner basis, must be in want of an algorthim for calculating a reduced Groebner basis.

It seems odd that i can't find something in the Groebner package - if there isn't something there, I assume that there is a well known piece of code for doing this!

How to Graph a Permutation...

Asked by: Les 190
graphtheory tree
December 23 2017
2  11

Maple 2015

Using with(combinat) the permutation of {a,b,c} is determined.

>restart:
>with(combinat):
>permute({a, b, c})
                  [[a, b, c], [a, c, b], [b, a, c], [b, c, a], [c, a, b], [c, b, a]]

The tree diagram of this permutation is

    

In Maple, using with(combinat) and with(GraphTheory), when I attempt to draw the permutation I get the following error:

>L := permute({a, b, c});
       L := [[a, b, c], [a, c, b], [b, a, c], [b, c, a], [c, a, b], [c, b, a]]
>DrawGraph(L);
  Error, invalid input: GraphTheory:-DrawGraph expects its 1st argument, H, to be of type       {GRAPHLN, list(GRAPHLN), set(GRAPHLN)}, but received [[a, b, c], [a, c, b], [b, a, c], [b, c,      a], [c,   a, b], [c, b, a]]

On Maple, again using with(combinat) and with(GraphTheory) the command permute(3) is used.  The results are manually configured as node-connection lines.  A fair representation of the tree diagram is configured by Maple, although the diagram has numeric instead of alpha configurations, and the a,b,c structure shown above is not easily recognized.

Any suggestions on developing a procedure that will graph (draw) an alpha-labeled permutation welcomed.  Thanks!  WC44_Permutation_Graph.mw

Problem in solving linear system of equations...

Asked by: Muhammad Usman 245
linear solve
December 22 2017
0  5

Dearz

Hope you would be fine with everything. I try to solve the following linear system of equations via fsolve command but the solution doesn't satisfied the system please see and put your valueable comments. Waiting your positive response.

-5.7167551941125971285 d[1, 1] - 0.23520507704562101132 d[1, 2]

   - 4.7759348859301130832 d[1, 3]

   + 82.882747548740738074 d[1, 4]

   + 1.5473302855836067493 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.2926823766365742120 d[2, 3]

   - 22.433527600870893213 d[2, 4]

   - 11.906076336447024126 d[3, 1]

   - 0.48985298599265354856 d[3, 2]

   - 9.9466643924764099316 d[3, 3]

   + 172.61685795222431091 d[3, 4]

   + 153.42462622364681378 d[4, 1]

   - 17.156128463674125233 d[4, 2]

   + 222.04914007834331471 d[4, 3]

   - 2162.1913920527683546 d[4, 4] = 0
-6.3505370802317673052 d[1, 1] - 0.23520507704562101132 d[1, 2]

   - 5.4097167720492832599 d[1, 3]

   + 54.802782951629695640 d[1, 4]

   + 1.7188733853924696263 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.4642254764454370890 d[2, 3]

   - 14.833240696164645293 d[2, 4]

   - 13.226030621801645811 d[3, 1]

   - 0.48985298599265354856 d[3, 2]

   - 11.266618677831031617 d[3, 3]

   + 114.13574573628827681 d[3, 4]

   + 107.19584752215150208 d[4, 1]

   - 17.156128463674125233 d[4, 2]

   + 175.82036137684800302 d[4, 3]

   - 1136.3239123361047712 d[4, 4] = 0
-6.7642088272251297212 d[1, 1] - 0.23520507704562101132 d[1, 2]

   - 5.8233885190426456759 d[1, 3]

   + 34.632657184619275137 d[1, 4]

   + 1.8308401918550417305 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.5761922829080091932 d[2, 3] - 9.373876878125528749 d[2, 4]

   - 14.087569594645296643 d[3, 1]

   - 0.48985298599265354856 d[3, 2]

   - 12.128157650674682449 d[3, 3]

   + 72.128164697121390121 d[3, 4]

   + 77.022155175221117487 d[4, 1]

   - 17.156128463674125233 d[4, 2]

   + 145.64666902991761843 d[4, 3]

   - 601.11088029977885095 d[4, 4] = 0
-1.5473302855836067487 d[1, 1] - 0.06366197723675813430 d[1, 2]

   - 1.2926823766365742115 d[1, 3]

   + 22.433527600870893218 d[1, 4]

   + 1.5473302855836067493 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.2926823766365742120 d[2, 3]

   - 22.433527600870893213 d[2, 4]

   - 7.7366514279180337465 d[3, 1]

   - 0.31830988618379067154 d[3, 2]

   - 6.4634118831828710599 d[3, 3]

   + 112.16763800435446606 d[3, 4]

   + 104.66490008068725185 d[4, 1]

   - 19.162255148264198426 d[4, 2]

   + 181.31392067374404557 d[4, 3]

   - 1455.2623850848598494 d[4, 4] = 0
-1.7188733853924696257 d[1, 1] - 0.06366197723675813430 d[1, 2]

   - 1.4642254764454370885 d[1, 3]

   + 14.833240696164645297 d[1, 4]

   + 1.7188733853924696263 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.4642254764454370890 d[2, 3]

   - 14.833240696164645293 d[2, 4]

   - 8.5943669269623481316 d[3, 1]

   - 0.31830988618379067154 d[3, 2]

   - 7.3211273822271854450 d[3, 3]

   + 74.166203480823226458 d[3, 4]

   + 53.030427038219525869 d[4, 1]

   - 19.162255148264198426 d[4, 2]

   + 129.67944763127631958 d[4, 3]

   - 668.89639482661771723 d[4, 4] = 0
-1.8308401918550417299 d[1, 1] - 0.06366197723675813430 d[1, 2]

   - 1.5761922829080091926 d[1, 3] + 9.373876878125528754 d[1, 4]

   + 1.8308401918550417305 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.5761922829080091932 d[2, 3] - 9.373876878125528749 d[2, 4]

   - 9.1542009592752086523 d[3, 1]

   - 0.31830988618379067154 d[3, 2]

   - 7.8809614145400459657 d[3, 3]

   + 46.869384390627643742 d[3, 4]

   + 19.328418292985322519 d[4, 1]

   - 19.162255148264198426 d[4, 2]

   + 95.977438886042116228 d[4, 3]

   - 305.71973224709969080 d[4, 4] = 0
 7.0561523113686303394 d[1, 1] - 1.9098593171027440292 d[1, 2]

    + 14.695589579779606456 d[1, 3]

    - 96.471127562654332340 d[1, 4]

    - 2.3520507704562101132 d[2, 1]

    + 0.63661977236758134308 d[2, 2]

    - 4.8985298599265354856 d[2, 3]

    + 32.157042520884777447 d[2, 4]

    + 16.464355393193470792 d[3, 1]

    - 4.4563384065730694016 d[3, 2]

    + 34.289709019485748399 d[3, 3]

    - 225.09929764619344213 d[3, 4]

    - 96.434081588704614639 d[4, 1]

    + 26.101410667070835066 d[4, 2]

    - 200.83972425698795490 d[4, 3]

    + 1318.4387433562758754 d[4, 4] = 0
-2.3520507704562101132 d[1, 1] + 0.6366197723675813431 d[1, 2]

   - 4.898529859926535486 d[1, 3] + 32.157042520884777450 d[1, 4]

   - 2.3520507704562101132 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 4.8985298599265354856 d[2, 3]

   + 32.157042520884777447 d[2, 4]

   + 7.0561523113686303394 d[3, 1]

   - 1.9098593171027440293 d[3, 2]

   + 14.695589579779606457 d[3, 3] - 96.47112756265433234 d[3, 4]

   - 11.760253852281050559 d[4, 1] + 3.183098861837906715 d[4, 2]

   - 24.49264929963267742 d[4, 3] + 160.7852126044238874 d[4, 4] = 

  1
1.9098593171027440291 d[1, 1] - 1.9098593171027440292 d[1, 2]

   + 9.5492965855137201456 d[1, 3]

   - 36.287327024952136554 d[1, 4]

   - 0.6366197723675813430 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 3.1830988618379067154 d[2, 3]

   + 12.095775674984045518 d[2, 4]

   + 4.4563384065730694010 d[3, 1]

   - 4.4563384065730694016 d[3, 2]

   + 22.281692032865347008 d[3, 3] - 84.67042972488831863 d[3, 4]

   - 26.101410667070835067 d[4, 1]

   + 26.101410667070835066 d[4, 2]

   - 130.50705333535417533 d[4, 3]

   + 495.92680267434586630 d[4, 4] = 0
-0.6366197723675813431 d[1, 1] + 0.6366197723675813431 d[1, 2]

   - 3.1830988618379067164 d[1, 3]

   + 12.095775674984045516 d[1, 4]

   - 0.6366197723675813430 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 3.1830988618379067154 d[2, 3]

   + 12.095775674984045518 d[2, 4]

   + 1.9098593171027440288 d[3, 1]

   - 1.9098593171027440293 d[3, 2] + 9.549296585513720146 d[3, 3]

   - 36.287327024952136560 d[3, 4] - 3.183098861837906717 d[4, 1]

   + 3.183098861837906715 d[4, 2] - 15.91549430918953358 d[4, 3]

   + 60.47887837492022764 d[4, 4] = 1
-1.4491448767744190950 d[1, 1] - 1.9098593171027440292 d[1, 2]

   + 6.1902923916365570215 d[1, 3]

   - 11.964006709004497915 d[1, 4]

   + 0.4830482922581396984 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 2.0634307972121856740 d[2, 3] + 3.988002236334832639 d[2, 4]

   - 3.381338045806977889 d[3, 1] - 4.4563384065730694016 d[3, 2]

   + 14.444015580485299718 d[3, 3] - 27.91601565434382847 d[3, 4]

   + 19.804979982583727629 d[4, 1]

   + 26.101410667070835066 d[4, 2]

   - 84.600662685699612634 d[4, 3] + 163.5080916897281382 d[4, 4] = 

  0
0.4830482922581396984 d[1, 1] + 0.6366197723675813431 d[1, 2]

   - 2.0634307972121856744 d[1, 3] + 3.988002236334832645 d[1, 4]

   + 0.4830482922581396984 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 2.0634307972121856740 d[2, 3] + 3.988002236334832639 d[2, 4]

   - 1.4491448767744190956 d[3, 1]

   - 1.9098593171027440293 d[3, 2]

   + 6.1902923916365570221 d[3, 3] - 11.96400670900449791 d[3, 4]

   + 2.415241461290698491 d[4, 1] + 3.183098861837906715 d[4, 2]

   - 10.317153986060928369 d[4, 3] + 19.94001118167416332 d[4, 4] = 

  1
 11.581726419330485018 d[1, 1] - 3.8605754731101616728 d[1, 2]

    + 27.024028311771131709 d[1, 3]

    - 158.28359439751662858 d[1, 4]

    - 1.9098593171027440292 d[2, 1]

    + 0.63661977236758134308 d[2, 2]

    - 4.4563384065730694016 d[2, 3]

    + 26.101410667070835066 d[2, 4]

    + 19.221163687741461135 d[3, 1]

    - 6.4070545625804870452 d[3, 2]

    + 44.849381938063409316 d[3, 3]

    - 262.68923706579996884 d[3, 4]

    - 172.31418534244454203 d[4, 1]

    + 57.438061780814847345 d[4, 2]

    - 402.06643246570393142 d[4, 3]

    + 2354.9605330134087411 d[4, 4] = 0
 7.0561523113686303394 d[1, 1] - 2.3520507704562101132 d[1, 2]

    + 16.464355393193470792 d[1, 3]

    - 96.434081588704614639 d[1, 4]

    - 1.9098593171027440292 d[2, 1]

    + 0.63661977236758134308 d[2, 2]

    - 4.4563384065730694016 d[2, 3]

    + 26.101410667070835066 d[2, 4]

    + 14.695589579779606456 d[3, 1]

    - 4.8985298599265354856 d[3, 2]

    + 34.289709019485748399 d[3, 3]

    - 200.83972425698795490 d[3, 4]

    - 96.471127562654332340 d[4, 1]

    + 32.157042520884777447 d[4, 2]

    - 225.09929764619344213 d[4, 3]

    + 1318.4387433562758753 d[4, 4] = 0
 1.9098593171027440291 d[1, 1] - 0.6366197723675813430 d[1, 2]

    + 4.4563384065730694010 d[1, 3]

    - 26.101410667070835067 d[1, 4]

    - 1.9098593171027440292 d[2, 1]

    + 0.63661977236758134308 d[2, 2]

    - 4.4563384065730694016 d[2, 3]

    + 26.101410667070835066 d[2, 4]

    + 9.5492965855137201456 d[3, 1]

    - 3.1830988618379067154 d[3, 2]

    + 22.281692032865347008 d[3, 3]

    - 130.50705333535417533 d[3, 4]

    - 36.287327024952136554 d[4, 1]

    + 12.095775674984045518 d[4, 2]

    - 84.670429724888318626 d[4, 3]

    + 495.92680267434586626 d[4, 4] = 0
-1.4491448767744190950 d[1, 1] + 0.4830482922581396984 d[1, 2]

   - 3.381338045806977889 d[1, 3] + 19.804979982583727629 d[1, 4]

   - 1.9098593171027440292 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 4.4563384065730694016 d[2, 3]

   + 26.101410667070835066 d[2, 4]

   + 6.1902923916365570215 d[3, 1]

   - 2.0634307972121856740 d[3, 2]

   + 14.444015580485299718 d[3, 3]

   - 84.600662685699612634 d[3, 4]

   - 11.964006709004497917 d[4, 1] + 3.988002236334832638 d[4, 2]

   - 27.91601565434382847 d[4, 3] + 163.50809168972813819 d[4, 4] = 

  0
Sols := solve([seq(`$`(F1[l1, l2], l1 = 2 .. 2^K*M-1), l2 = 2 .. 2^K*M), seq(`$`(F2[l1], l1 = 2 .. 2^K*M)), seq(`$`(F3[l1], l1 = 2 .. 2^K*M)), seq(`$`(F4[l1], l1 = 1 .. 2^K*M))], {seq(`$`(d[l1, l2], l1 = 1 .. 2^K*M), l2 = 1 .. 2^K*M)});
map(evalf, subs(Sols, convert(F4, list)));

 

Problem in finding solution of system PDEs...

Asked by: Muhammad Usman 245
numeric differential-equations
December 16 2017
0  3

Dearz!

Hope everyone is fine with everything. I am facing problem to solve the system of PDEs in the attached file. Is there any built-in command to the solve the attached system of PDEs via FEM, FDM, SIMPLER algorithm or some other efficient method? Please try to fix my problem. I am waiting your positive response. Thanks in advance.

PDEs_sol.mw

How to find a norm of polynomial with fractional r...

Asked by: Jameel123 20
polynomial norm
December 15 2017
1  2

problem.mw
 

> 

restart
> 

with(LinearAlgebra):
> 

with(orthopoly):
> 

with(student):
> 

> 

Digits := 32

32

 (1)

``

 (2)
> 

aa := 92.359325848696162368355511315798*x^4+7.1870894227877814492593273580981*x^2-6.8000000000000000000000000000000*10^(-32)*x-34.149715256645549592146037526614*x^3-135.58601962203264776595890928687*x^5+101.25962748430115463423633295054*x^6-30.114080251339988529197699897019*x^7-x^1.25

92.359325848696162368355511315798*x^4+7.1870894227877814492593273580981*x^2-0.68000000000000000000000000000000e-31*x-34.149715256645549592146037526614*x^3-135.58601962203264776595890928687*x^5+101.25962748430115463423633295054*x^6-30.114080251339988529197699897019*x^7-x^1.25

 (3)
> 

norm(aa, infinity)

Error, invalid input: norm expects its 1st argument, p, to be of type {Matrix, Vector, matrix, polynom, vector}, but received 92.359325848696162368355511315798*x^4+7.1870894227877814492593273580981*x^2-0.68000000000000000000000000000000e-31*x-34.149715256645549592146037526614*x^3-135.58601962203264776595890928687*x^5+101.25962748430115463423633295054*x^6-30.114080251339988529197699897019*x^7-x^1.25

  
> 

``


 

Download problem.mw

 

First 22 23 24 25 26 27 28 Last Page 24 of 73