## 13 Badges

11 years, 87 days

## your question is not as clear enough !...

when u assign sth to a symbol it gets its value and when u execute that,it will return the value not itself ! use eval instead !

 >
 >
 (1)
 >
 (2)
 >

Download eval.mws

## please insert your data more precisley...

here is one way :

 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >

Download matrix.mw

## do u mean this ?...

 > restart:q:=x+y;
 (1)
 > f := unapply(q, x, y);c:=10^(-1):
 (2)
 > G(x,y,xi,eta,t):=(1+2*(sum((exp(-Pi^(2)*n^(2)*c*t))*cos(n*Pi*x)*cos(n*Pi*xi),n=1..10)))*(1+2*(sum((exp(-Pi^(2)*m^(2)*c*t))*cos(m*Pi*y)*cos(m*Pi*eta),m=1..10)));
 (3)
 > w:=int(int(f(xi,eta)*G(x,y,xi,eta,t), eta=0..1), xi=0..1);W:=unapply(w,x,y,t);
 (4)
 > plot(W(.51086, .49427, t), t = 0 .. 1);
 > plottools:-getdata(%); M := %[-1]; #ExportMatrix("new2.dat", M);
 (5)
 > M(1..5,1..2);
 (6)
 > for i to 200 do X||i := M(i,1);Y||i:=M(i,2);od;
 (7)
 >

Download loop.mw

## use simplify...

 > restart: f:=simplify(cosh(t)/(cosh((17/15)*t)+cosh(t)));
 (1)
 > evalf(Int(f, t = 0 .. infinity));
 >
 >
 > evalf(int(f, t = 0 .. infinity));
 >
 >
 (2)
 >

Download int.mw

## two ways...

 > restart:
 > sort(x+y+z, order = tdeg(z, y, x));
 (1)
 > sort( x+y+z, order = plex(z, y, x));
 (2)
 >

Download one_way.mw

## one way...

 > restart;
 > with(Physics):
 >
 > Setup(mathematicalnotation = true):
 > Setup(quantumop = {Q,P}, algebrarule = {%AntiCommutator(Q[i], P[k]) = 2*KroneckerDelta[i, k], %Commutator(Q[j], P[k]) = 2*I*(add(LeviCivita[j, k, l]*Q[l], l = 1 .. 3))}):
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >

Download corrected.mw

## see ?userinfo...

from maple help page :

userinfo - print useful information to the user

Calling Sequence
userinfo(lev, fn, e1, e2 ... );

 The user must assign a non-negative integer to some of the entries in the global table infolevel before invoking the procedure. If the entry infolevel[all] is a non-negative integer then every userinfo call will print if its level is less than or equal to infolevel[all] .

 Throughout the Maple library userinfo statements have been used with the following conventions:
Level 1: reserved for information that the user must be told.
Level 2,3: general information, including technique or algorithm being used.
Level 4,5: more detailed information about how the problem is being solved

The item infolevel[hints] is initialized to a value 1. Maple sometimes returns unevaluated answers when it does not have enough information to produce an explicit answer (because, for example, such an answer would not be correct over all complex numbers). This facility is intended to report hints as to which further information, given through assume(), would enable Maple return an explicit answer.

## as preben says,every thing is ok...

actucallyA is, A=0.5+sqrt(3)/2*I .and all (x+y) can simplify from Numerator and denominator of your A and actually A is independent of x and y.

 A := (1/2)*(-x-y+(-3*(x+y)^2)^(1/2))/(x+y)

 x := 10+I y := 10-I

 -1/2+(1/40)*(-1200)^(1/2) -.5000000000+.8660254038*I

 x := 10+10*I y := 10-I (10/481-(9/962)*I)*((-20-9*I)+(-957-1080*I)^(1/2)) -.5000000000-.8660254038*I

 x := 10+5*2^(1/2)+(5*I)*2^(1/2) y := 10-I (1/2)*((-20+I)-5*2^(1/2)-(5*I)*2^(1/2)+(-3*(20-I+5*2^(1/2)+(5*I)*2^(1/2))^2)^(1/2))/(20-I+5*2^(1/2)+(5*I)*2^(1/2)) -.5000000000-.8660254035*I

Download correct.mw

## make local in proc...

in printf u should write everyhting in string , "Ghana Chocolates" . and make local varibales in the procedure. good luck

 Factory Simulation: Ghana Chocolates Chocolate Factory Simulation Project .     W    O(t,t+1)        B(t)        T(t)        F(t)        R(t)     ---------------------------------------------------------------------------     1         50.         50.        250.        250.        150.        2        100.        100.        200.        200.        200.        3          0.          0.        300.        300.         50.        4         50.         50.          0.          0.        400.        5         50.        100.        125.        125.        325.        6         50.         50.        120.        120.        280.        7         50.         50.        117.        117.        283.        8         50.         50.        114.        114.        286.        9         50.         50.        112.        112.        288.       10         50.         50.        111.        111.        289.     ---------------------------------------------------------------------------     W          50      875/18          50       ---------------------------------------------------------------------------

Download suppress_warnins.mw

## use two loops...

i write a procedure , which M is input for the procedure.

 > restart:
 > mat:=proc(M) local A,i,j;A:=Matrix(M,M): for i to M do for j to M do A(i,j):=(psi||1||(i-1))((2*j-1)/(2*m)); od;od;end proc;
 (1)
 > mat(3);
 (2)
 >
 >

Download Mat.mw

## i have answered exactly the same questio...

http://www.mapleprimes.com/questions/201345-Geom3d-Seq-Of-Point--In-R3

 > restart:with(plots):
 >
 > Tab := [[0, 0, 0], [1/5, 0, 0], [2/5, 0, 0], [0, 1/5, 0], [1/5, 1/5, 0], [2/5, 1/5, 0], [0, 2/5, 0], [1/5, 2/5, 0], [2/5, 2/5, 0], [0, 0, 1/5], [1/5, 0, 1/5], [2/5, 0, 1/5], [0, 1/5, 1/5], [1/5, 1/5, 1/5], [2/5, 1/5, 1/5], [0, 2/5, 1/5], [1/5, 2/5, 1/5], [2/5, 2/5, 1/5], [0, 0, 2/5], [1/5, 0, 2/5], [2/5, 0, 2/5], [0, 1/5, 2/5], [1/5, 1/5, 2/5], [2/5, 1/5, 2/5], [0, 2/5, 2/5], [1/5, 2/5, 2/5], [2/5, 2/5, 2/5]];
 (1)
 > pointplot3d(Tab,color = magenta,axes = normal, symbol = box);
 >

Download plot3d.mw

## use pointplot3d...

 > restart:with(plots):
 > x := [seq]( (1/5)*i,i=0..12);  #  x[i] the x-coordinate
 > y := [seq]( (1/5)*j,j=0..12); # y[j] the y-coordinate
 > t := [seq]( (1/5)*k,k=0..12);  #  t[k] the t-coordinate
 >
 (1)
 > listplot3d([x,y,t]);
 > pointplot3d([seq]([op(i,x),op(i,y),op(i,t)],i=1..nops(x)),color = magenta,axes = normal, symbol = box);
 >

Download plot.mw

## remove print from your function...

 Fun := proc (x) options operator, arrow; x^2 end proc

 A := [1, 2, 3, 4]

 [1, 4, 9, 16] 1, 4, 9, 16

Download map.mw

## stats[describe, quadraticmean](data)...

it is for 15 number of data.

maple.mws

testfile.txt

## your boundary condition are incorrect...

az u can see from the wroksheet i uploaded, u can solve your problem without your boundary conditions, but with them u can not .
i changed your way of typing , and also changed theta[p] to sigma .
u have write G(10)=-f(10) while you have note assigned f(10):=0 , so never use this . also u have written L:=[0.2]; and i do not know where u have used it in your equations.

dsolve without your boundary take alot of memory and time, so i terminated the process. but your problem,is with your boundary conditions,and actually i do not know ho to correct them. maybe preben alsholm can help. good luck

 > restart:
 > fixedparameter := {M = .5, B = .5, theta[r] = -10, L0 = 1, s = .1, Pr = 1}:
 > Eq1 := eval((1-theta(eta)/theta[r])*(diff(f(eta), eta, eta, eta))+(diff(f(eta), eta, eta))*(diff(theta(eta), eta))/theta[r]+(1-theta(eta)/theta[r])^2*(f(eta)*(diff(f(eta), eta, eta))-(diff(f(eta), eta))^2-M*(diff(f(eta), eta))+B*H(eta)*(F(eta)-(diff(f(eta), eta)))),fixedparameter);
 (1)
 > Eq2 := eval(G(eta)*(diff(F(eta), eta))+F(eta)^2+B*(F(eta)-(diff(f(eta), eta))) ,fixedparameter);
 (2)
 > Eq3 := eval(G(eta)*(diff(G(eta), eta))+B*(f(eta)+G(eta)) ,fixedparameter);
 (3)
 > Eq4 := eval(G(eta)*(diff(H(eta), eta))+H(eta)*(diff(G(eta), eta))+F(eta)*H(eta),fixedparameter);
 (4)
 > Eq5 := eval((1+s*theta(eta))*(diff(theta(eta), eta, eta))+(diff(theta(eta), eta))^2*s+Pr*(f(eta)*(diff(theta(eta), eta))-(diff(f(eta), eta))*theta(eta))+(2/3)*B*H(eta)*(sgima(eta)-theta(eta)),fixedparameter);
 (5)
 > Eq6 := eval(2*F(eta)*sgima(eta)+G(eta)*(diff(sgima(eta), eta))+L0*B*(sgima(eta)-theta(eta)) ,fixedparameter);
 (6)
 > bcs := {f(0) = 0, (D(f))(0) = 1, (D(f))(10) = 0,F(10) = 0,G(10) = 0,H(10) = n,theta(10) = 0,sgima(10) = 0};
 (7)
 > L := [0.2]:
 >
 > dsolve({seq(Eq||i,i=1..6)},bcs);dsolve({seq(Eq||i,i=1..6)},bcs,numeric);
 > dsolve({seq(Eq||i,i=1..6)});
 > #for k to 1 do R := dsolve(eval({Eq10, Eq11, Eq4, Eq7, Eq8, Eq9, bcs1, bcs2, bcs3, bcs4, bcs5, bcs6}, n = L[k]), [f(eta), F(eta), G(eta), H(eta), theta(eta), sgima(eta)], numeric, output = listprocedure); Y || k := rhs(R[5]); YP || k := rhs(R[6]); YJ || k := rhs(R[7]); YS || k := rhs(R[2]) end do
 >

Download dsolve.mw

 4 5 6 7 8 9 Page 6 of 9
﻿