Maple 13 Questions and Posts

These are Posts and Questions associated with the product, Maple 13

Dear sir, can we get a plot of graph with line style like  line+symbol so on.. 

I am attaching codes and sample graph .

 

h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=(4/h(z)^4)-(sin(alpha)/F)-h(z)^2+Nb*h(z)^4:
lambda:=(F,Nb,delta2)->Int(K1,z=0..1):

L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
A:=plot( [seq(seq(lambda(F,Nb,delta2), Nb=[0.1,0.2,0.3]), F=[0.1,0.2,0.5])],
          delta2=0.02..0.1, linestyle = [solid,longdash,dashdot],'
          thickness = 2',color=[red$3,blue$3,black$3],
          legend=["Nb=0.1","Nb=0.2","Nb=0.3",""$6]):
B:=plots:-textplot([[0.05,-1,F=0.1],[0.05,1.5,F=0.2], [0.05,3,F=0.5]], font=[times, 14]):

plots:-display(A, B);

Dear sir i want to display only three curves 1, curve 2 curve 3  in legends but in the plot it should be nine curve a sampe codes and a sample graph is attaching , please give suggestions. Here is my codes 

h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=(4/h(z)^4)-(sin(alpha)/F)-h(z)^2+Nb*h(z)^4:
lambda:=(F,Nb,delta2)->Int(K1,z=0..1):

L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
A:=plot( [seq(seq(lambda(F,Nb,delta2), Nb=[0.1,0.2,0.3]), F=[0.1,0.2,0.5])], delta2=0.02..0.1, linestyle = [solid,longdash,dashdot],'thickness = 2',color=[red$3,blue$3,black$3]):

B:=plots:-textplot([[0.05,-1,F=0.1],[0.05,1.5,F=0.2], [0.05,3,F=0.5]], font=[times, 14]):

plots:-display(A, B);
 sample graphs 

Hellow sir how to plot stream lines for different values of alpha =6,7,8,9.

 I am attaching the sample graph and codes.
How to get the output like that only.

Thanks in advaned

 

eq1:= alpha*diff(f(eta),eta$3)+f(eta)*diff(f(eta),eta$2)-diff(f(eta),eta)^2+1:
eq2:= diff(h(eta),eta$2)+f(eta)*alpha*diff(h(eta),eta)-diff(f(eta),eta)*h(eta):
eq3:= diff(H(eta),eta) = h(eta):
bc:= f(0)=0, D(f)(0)=0, D(f)(6)=1, h(0)=0, D(h)(6)=1, H(0)=0:
A1:= dsolve(
     {eq||(1..3), bc}, 
     numeric, method= bvp[midrich], abserr = 1e-10, output= operator
):

psi:= unapply(eval(x*f(eta)+H(eta), A1), (x,eta)):
plots:-contourplot((psi(x,eta), x= 0..6, eta= 0..6),alpha=6);

The summation takes too long time. Please help me
 

 

 

 

 

equ.pdf The equation is in the following form:

 

Dear sir I want to plot the graph where the lines are coincides maximize the graph in the same graph.  A sample graph i am attaching .some sample codes are here.

h:=z->1-(gamma2/2)*(1 + cos(2*(Pi/omicron1)*(z - d1 - omicron1))):
K1:=((4/h(z)^4)-(sin(Zeta1)/omega)-h(z)^2+Nb*h(z)^4):
kappa:=Int(K1,z=0..1): 
omicron1:=0.2: omega:=10:d1:=0.2:Zeta1:=Pi/6:

plot( [seq(eval(kappa, Nb=j), j in [0.01,0.02,0.03])], gamma2=0.02..0.1);



 

How can i print maple file? when ever i press ctrl+ p it never responce any thing can we print maple file directly? if yes then please how?

hellow can we draw the figure like this


 

``

restart

``

f := proc (x, z) options operator, arrow; ln(x^z) end proc;

proc (x, z) options operator, arrow; ln(x^z) end proc

(1)

``

``

``

simplify(f(x, 2), ln);

ln(x^2)

(2)

simplify(f(x, 1/2+I*y), ln);

ln(x^(-((1/2)*I)*(I-2*y)))

(3)

``


 

Download simplln.mw

 

with(CodeGeneration)

ua=[ 1 2 3 4];

Matlab(ua);

 

%How can i write the Output of "Matlab(ua)"  on file so I can use the Matlab "load" command
later on to get the result in Matlab?  

 

 

I am looking for a window bat file to run my maple input file.mw

I need to convert this Maple Table to a Matrix

table([(25, 1) = -39, (16, 151) = 32, (33, 1) = -54, (1, 1) = 29, (13, 1) = 32, (31, 101) = -7, (6, 51) = -10, (11, 101) = -1, (28, 151) = -39, (18, 51) = -65, (4, 151) = 29, (8, 151) = -10, (23, 101) = 23, (34, 51) = -54, (40, 151) = 87, (36, 151) = -54, (9, 1) = -1, (37, 1) = 87, (21, 1) = 23, (14, 51) = 32, (22, 51) = 23, (20, 151) = -65, (27, 101) = -39, (3, 101) = 29, (19, 101) = -65, (24, 151) = 23, (32, 151) = -7, (30, 51) = -7, (38, 51) = 87, (7, 101) = -10, (10, 51) = -1, (29, 1) = -7, (35, 101) = -54, (17, 1) = -65, (26, 51) = -39, (15, 101) = 32, (12, 151) = -1, (39, 101) = 87, (5, 1) = -10, (2, 51) = 29]) 

We conjecture that the polynomial h(n) = n^2 + n + 41 is prime for an infinite number of values n.
We furthur conjecture that p(n) = n^2 + 1 is prime an infinite number of times.

I have shown that the set (x,y) with h(y) mod x is congruent to 0 can be written down.  It is p(x,y).  p(x,y) is the set of all divisors of h(n).  See

https://sites.google.com/site/primeproducingpolynomial/

landau.mw

Regards,

Matt

Dear sir,

 i want to plot a graph with  tick marks on all sides , i am attaching codes and  a sample graph, see the attachmnent,

The graph range should be displace only on x and y axis and the vaues should not be on opposite to x and y axis

please see the sample graph

restart:
h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/F)-h(z)^2+Nb*h(z)^4):
lambda:=Int(K1,z=0..1):
F:=0.3:
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
plot( [seq(eval(lambda, Nb=j), j in [0.1,0.2,0.3])], delta2=0.02..0.1);

 

hai,

 i want to plot a graph here is my codes and i want the graph like this, can any body help me

 

restart:
with(DEtools):
with(plots):

h(z):=1-(delta2/2)*(1 + Cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h^4)-(sin(alpha)/F)-h^2+Nb*h^4):
lambda:=int(K1,z=0..1):
F:=0.3:
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
evalf(lambda);

First 7 8 9 10 11 12 13 Last Page 9 of 54