Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

hi..

is correct this answer for differential equations??

i think order of result should be in (10^6 or 10^9 or higher) range....

please check it

thanks

hpp.mw
 

restart

L := 100*10^(-9):

Eq1 := {-(1017/1600000000000000000000000000000000000000000)*(diff(w(x), x, x, x, x, x, x))+(26169/40000000000000000000000000)*(diff(w(x), x, x, x, x))-0.8325000000e-4*omega^2+1.560937500*10^(-21)*omega^2*(diff(w(x), x, x)), w(0) = 0, w(1/10000000) = 0, (D(w))(0) = 0, (D(w))(1/10000000) = 0, ((D@@2)(w))(0) = 0, ((D@@2)(w))(1/10000000) = 0}:

sys := subs(omega^2 = omega2, Eq1);

{-(1017/1600000000000000000000000000000000000000000)*(diff(diff(diff(diff(diff(diff(w(x), x), x), x), x), x), x))+(26169/40000000000000000000000000)*(diff(diff(diff(diff(w(x), x), x), x), x))-0.8325000000e-4*omega2+0.1560937500e-20*omega2*(diff(diff(w(x), x), x)), w(0) = 0, w(1/10000000) = 0, (D(w))(0) = 0, (D(w))(1/10000000) = 0, ((D@@2)(w))(0) = 0, ((D@@2)(w))(1/10000000) = 0}

 

{-(5085/8)*(diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y))+6542250*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.8325000000e-4*omega2+0.1560937500e-6*omega2*(diff(diff(g1(y), y), y)), 10000000*(D(g1))(0) = 0, 10000000*(D(g1))(1) = 0, 100000000000000*((D@@2)(g1))(0) = 0, 100000000000000*((D@@2)(g1))(1) = 0, g1(0) = 0, g1(1) = 0}

 

{-(5085/8)*(diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y))+6542250*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.8325000000e-4*omega2+0.1560937500e-6*omega2*(diff(diff(g1(y), y), y))}, {10000000*(D(g1))(0) = 0, 10000000*(D(g1))(1) = 0, 100000000000000*((D@@2)(g1))(0) = 0, 100000000000000*((D@@2)(g1))(1) = 0, g1(0) = 0, g1(1) = 0}

 

{g1(0) = 0, g1(1) = 0, (D(g1))(0) = 0, (D(g1))(1) = 0, ((D@@2)(g1))(0) = 0, ((D@@2)(g1))(1) = 0}

 

{diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y), diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), diff(diff(diff(diff(g1(y), y), y), y), y), diff(diff(diff(g1(y), y), y), y), diff(diff(g1(y), y), y), diff(g1(y), y)}

 

{-(1017/1600000000000000000000000000000000000000000)*(diff(diff(diff(diff(diff(diff(w(x), x), x), x), x), x), x))+(26169/40000000000000000000000000)*(diff(diff(diff(diff(w(x), x), x), x), x))-0.8325000000e-4*omega2+0.1560937500e-20*omega2*(diff(diff(w(x), x), x))}

 

{-(5085/8)*(diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y))+6542250*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.8325000000e-4*omega2+0.1560937500e-6*omega2*(diff(diff(g1(y), y), y))}

 

{diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y) = 10292.62537*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.1309734513e-6*omega2+0.2455752212e-9*omega2*(diff(diff(g1(y), y), y))}

 

{diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y) = 10292.62537*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.1309734513e13*omega3+2455752212.*omega3*(diff(diff(g1(y), y), y))}

 

{diff(diff(diff(diff(diff(diff(g1(y), y), y), y), y), y), y) = 10292.62537*(diff(diff(diff(diff(g1(y), y), y), y), y))-0.1309734513e13*omega3+2455752212.*omega3*(diff(diff(g1(y), y), y)), g1(0) = 0, g1(1) = 0, (D(g1))(0) = 0, (D(g1))(1) = 0, ((D@@2)(g1))(0) = 0, ((D@@2)(g1))(1) = 0}

 

{((D@@3)(g1))(0), ((D@@3)(g1))(1), ((D@@4)(g1))(0), ((D@@4)(g1))(1), ((D@@5)(g1))(0), ((D@@5)(g1))(1)}

 

((D@@3)(g1))(0)

 

((D@@3)(g1))(1)

 

((D@@4)(g1))(0)

 

((D@@4)(g1))(1)

 

((D@@5)(g1))(0)

 

((D@@5)(g1))(1)

 

((D@@5)(g1))(1), ((D@@4)(g1))(0), ((D@@5)(g1))(0), ((D@@3)(g1))(1), ((D@@4)(g1))(1), ((D@@3)(g1))(0)

 

HFloat(-8.852947665097804e-24), HFloat(-8.991820290300328e-22), HFloat(8.852947665097804e-24), HFloat(-9.672787782157173e-20), HFloat(-8.991820290300328e-22), HFloat(9.672787782157165e-20)

(1)

sqrt(8.85294766509780*10^(-21)*10^19);

.2975390338

(2)

NULL


 

Download hpp.mw

 

staganation_point11.mw
 

``

restart

l := 1:

1

 

1.5

 

.5

 

[blue, green, red, yellow]

(1)

``

for j to nops(A) do R1 := 2*n/(n+1); R2 := 2*p/(n+1); R3 := 2/(n+1); sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2)-M*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr*f(eta)*(diff(theta(eta), eta))-R2*pr*(diff(f(eta), eta))*theta(eta)+R3*(A[j]*(diff(f(eta), eta))+B*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L+b*((D@@2)(f))(0), (D(f))(20) = 1, theta(0) = 1+s*(D(theta))(0), theta(20) = 0], numeric, method = bvp); plots[odeplot](sol1, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol1, [eta, diff(f(eta), eta)], color = K[j], axes = boxed); tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = K[j], axes = normal) end do:

 

 

``

l := 1:

1

 

[0, 1, 1.5]

 

.5

 

[blue, green, red, yellow]

(2)

for j to nops(n1) do R4 := 2*n1[j]/(n1[j]+1); R5 := 2*p1/(n1[j]+1); R6 := 2/(n1[j]+1); sol2 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R4*(1-(diff(f(eta), eta))^2)-M1*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr1*f(eta)*(diff(theta(eta), eta))-R5*pr1*(diff(f(eta), eta))*theta(eta)+R6*(A1*(diff(f(eta), eta))+B1*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L1+b1*((D@@2)(f))(0), (D(f))(7) = 1, theta(0) = 1+s1*(D(theta))(0), theta(7) = 0], numeric, method = bvp); plots[odeplot](sol2, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol2, sol1, [[eta, diff(f(eta), eta)], [eta, diff(f(eta), eta)]], color = K1[j], axes = boxed); tplt[j] := plots[odeplot](sol2, [[eta, theta(eta)]], color = K1[j], axes = normal) end do; plots:-display([seq(fplt[j], j = 1 .. nops(n1))]); plots:-display([seq(tplt[j], j = 1 .. nops(n1))])

Error, (in plots/odeplot) invalid argument: sol1

 

 

 

``

``

sol1(0);

[eta = 0., f(eta) = 1.00000000000000044, diff(f(eta), eta) = .899599635987511914, diff(diff(f(eta), eta), eta) = -.200800728024976643, theta(eta) = 1.18657688243172332, diff(theta(eta), eta) = .373153764863445037]

(3)

sol1(.1)

[eta = .1, f(eta) = 1.08902313464617162, diff(f(eta), eta) = .881503890693141945, diff(diff(f(eta), eta), eta) = -.162309073227910134, theta(eta) = 1.22040489003745489, diff(theta(eta), eta) = .304232440930656767]

(4)

sol1(.2)

[eta = .2, f(eta) = 1.17641763749368966, diff(f(eta), eta) = .866916683092940898, diff(diff(f(eta), eta), eta) = -.130454301210374102, theta(eta) = 1.24759607023709362, diff(theta(eta), eta) = .240488988787701030]

(5)

sol1(.3)

[eta = .3, f(eta) = 1.19045803452309284, diff(f(eta), eta) = .579367537136023514, diff(diff(f(eta), eta), eta) = -.285675511621370782, theta(eta) = 1.18389221591022696, diff(theta(eta), eta) = .146013974567769960]

(6)

sol1(.4)

[eta = .4, f(eta) = 1.40000000000000034, diff(f(eta), eta) = 1.00000000000000022, diff(diff(f(eta), eta), eta) = -0.243774513041384287e-17, theta(eta) = .625958972186505536, diff(theta(eta), eta) = -.314549395236395634]

(7)

sol1(.5)

[eta = .5, f(eta) = 1.37026161183094430, diff(f(eta), eta) = .874752886901313142, diff(diff(f(eta), eta), eta) = .345911467377074400, theta(eta) = .432494259338694842, diff(theta(eta), eta) = -.382764248064397461]

(8)

sol1(.6)

[eta = .6, f(eta) = 1.36678221814533528, diff(f(eta), eta) = .771028661281065508, diff(diff(f(eta), eta), eta) = .407805382194403932, theta(eta) = .876413930517023876, diff(theta(eta), eta) = -.197648778495384870]

(9)

sol1(2)

[eta = 2., f(eta) = 2.66120522956795602, diff(f(eta), eta) = .991532161353848585, diff(diff(f(eta), eta), eta) = 0.251405465681268682e-1, theta(eta) = .635967939441598018, diff(theta(eta), eta) = -.144641270049362308]

(10)

``

``

``

 

``

``

NULL


 

Download staganation_point11.mw
 

``

restart

l := 1:

1

 

1.5

 

.5

 

[blue, green, red, yellow]

(1)

``

for j to nops(A) do R1 := 2*n/(n+1); R2 := 2*p/(n+1); R3 := 2/(n+1); sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2)-M*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr*f(eta)*(diff(theta(eta), eta))-R2*pr*(diff(f(eta), eta))*theta(eta)+R3*(A[j]*(diff(f(eta), eta))+B*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L+b*((D@@2)(f))(0), (D(f))(20) = 1, theta(0) = 1+s*(D(theta))(0), theta(20) = 0], numeric, method = bvp); plots[odeplot](sol1, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol1, [eta, diff(f(eta), eta)], color = K[j], axes = boxed); tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = K[j], axes = normal) end do:

 

 

``

l := 1:

1

 

[0, 1, 1.5]

 

.5

 

[blue, green, red, yellow]

(2)

for j to nops(n1) do R4 := 2*n1[j]/(n1[j]+1); R5 := 2*p1/(n1[j]+1); R6 := 2/(n1[j]+1); sol2 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R4*(1-(diff(f(eta), eta))^2)-M1*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr1*f(eta)*(diff(theta(eta), eta))-R5*pr1*(diff(f(eta), eta))*theta(eta)+R6*(A1*(diff(f(eta), eta))+B1*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L1+b1*((D@@2)(f))(0), (D(f))(7) = 1, theta(0) = 1+s1*(D(theta))(0), theta(7) = 0], numeric, method = bvp); plots[odeplot](sol2, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol2, sol1, [[eta, diff(f(eta), eta)], [eta, diff(f(eta), eta)]], color = K1[j], axes = boxed); tplt[j] := plots[odeplot](sol2, [[eta, theta(eta)]], color = K1[j], axes = normal) end do; plots:-display([seq(fplt[j], j = 1 .. nops(n1))]); plots:-display([seq(tplt[j], j = 1 .. nops(n1))])

Error, (in plots/odeplot) invalid argument: sol1

 

 

 

``

``

sol1(0);

[eta = 0., f(eta) = 1.00000000000000044, diff(f(eta), eta) = .899599635987511914, diff(diff(f(eta), eta), eta) = -.200800728024976643, theta(eta) = 1.18657688243172332, diff(theta(eta), eta) = .373153764863445037]

(3)

sol1(.1)

[eta = .1, f(eta) = 1.08902313464617162, diff(f(eta), eta) = .881503890693141945, diff(diff(f(eta), eta), eta) = -.162309073227910134, theta(eta) = 1.22040489003745489, diff(theta(eta), eta) = .304232440930656767]

(4)

sol1(.2)

[eta = .2, f(eta) = 1.17641763749368966, diff(f(eta), eta) = .866916683092940898, diff(diff(f(eta), eta), eta) = -.130454301210374102, theta(eta) = 1.24759607023709362, diff(theta(eta), eta) = .240488988787701030]

(5)

sol1(.3)

[eta = .3, f(eta) = 1.19045803452309284, diff(f(eta), eta) = .579367537136023514, diff(diff(f(eta), eta), eta) = -.285675511621370782, theta(eta) = 1.18389221591022696, diff(theta(eta), eta) = .146013974567769960]

(6)

sol1(.4)

[eta = .4, f(eta) = 1.40000000000000034, diff(f(eta), eta) = 1.00000000000000022, diff(diff(f(eta), eta), eta) = -0.243774513041384287e-17, theta(eta) = .625958972186505536, diff(theta(eta), eta) = -.314549395236395634]

(7)

sol1(.5)

[eta = .5, f(eta) = 1.37026161183094430, diff(f(eta), eta) = .874752886901313142, diff(diff(f(eta), eta), eta) = .345911467377074400, theta(eta) = .432494259338694842, diff(theta(eta), eta) = -.382764248064397461]

(8)

sol1(.6)

[eta = .6, f(eta) = 1.36678221814533528, diff(f(eta), eta) = .771028661281065508, diff(diff(f(eta), eta), eta) = .407805382194403932, theta(eta) = .876413930517023876, diff(theta(eta), eta) = -.197648778495384870]

(9)

sol1(2)

[eta = 2., f(eta) = 2.66120522956795602, diff(f(eta), eta) = .991532161353848585, diff(diff(f(eta), eta), eta) = 0.251405465681268682e-1, theta(eta) = .635967939441598018, diff(theta(eta), eta) = -.144641270049362308]

(10)

``

``

``

 

``

``

NULL


 

Download staganation_point11.mw
 

``

restart

l := 1:

1

 

1.5

 

.5

 

[blue, green, red, yellow]

(1)

``

for j to nops(A) do R1 := 2*n/(n+1); R2 := 2*p/(n+1); R3 := 2/(n+1); sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2)-M*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr*f(eta)*(diff(theta(eta), eta))-R2*pr*(diff(f(eta), eta))*theta(eta)+R3*(A[j]*(diff(f(eta), eta))+B*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L+b*((D@@2)(f))(0), (D(f))(20) = 1, theta(0) = 1+s*(D(theta))(0), theta(20) = 0], numeric, method = bvp); plots[odeplot](sol1, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol1, [eta, diff(f(eta), eta)], color = K[j], axes = boxed); tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = K[j], axes = normal) end do:

 

 

``

l := 1:

1

 

[0, 1, 1.5]

 

.5

 

[blue, green, red, yellow]

(2)

for j to nops(n1) do R4 := 2*n1[j]/(n1[j]+1); R5 := 2*p1/(n1[j]+1); R6 := 2/(n1[j]+1); sol2 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R4*(1-(diff(f(eta), eta))^2)-M1*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr1*f(eta)*(diff(theta(eta), eta))-R5*pr1*(diff(f(eta), eta))*theta(eta)+R6*(A1*(diff(f(eta), eta))+B1*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L1+b1*((D@@2)(f))(0), (D(f))(7) = 1, theta(0) = 1+s1*(D(theta))(0), theta(7) = 0], numeric, method = bvp); plots[odeplot](sol2, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol2, sol1, [[eta, diff(f(eta), eta)], [eta, diff(f(eta), eta)]], color = K1[j], axes = boxed); tplt[j] := plots[odeplot](sol2, [[eta, theta(eta)]], color = K1[j], axes = normal) end do; plots:-display([seq(fplt[j], j = 1 .. nops(n1))]); plots:-display([seq(tplt[j], j = 1 .. nops(n1))])

Error, (in plots/odeplot) invalid argument: sol1

 

 

 

``

``

sol1(0);

[eta = 0., f(eta) = 1.00000000000000044, diff(f(eta), eta) = .899599635987511914, diff(diff(f(eta), eta), eta) = -.200800728024976643, theta(eta) = 1.18657688243172332, diff(theta(eta), eta) = .373153764863445037]

(3)

sol1(.1)

[eta = .1, f(eta) = 1.08902313464617162, diff(f(eta), eta) = .881503890693141945, diff(diff(f(eta), eta), eta) = -.162309073227910134, theta(eta) = 1.22040489003745489, diff(theta(eta), eta) = .304232440930656767]

(4)

sol1(.2)

[eta = .2, f(eta) = 1.17641763749368966, diff(f(eta), eta) = .866916683092940898, diff(diff(f(eta), eta), eta) = -.130454301210374102, theta(eta) = 1.24759607023709362, diff(theta(eta), eta) = .240488988787701030]

(5)

sol1(.3)

[eta = .3, f(eta) = 1.19045803452309284, diff(f(eta), eta) = .579367537136023514, diff(diff(f(eta), eta), eta) = -.285675511621370782, theta(eta) = 1.18389221591022696, diff(theta(eta), eta) = .146013974567769960]

(6)

sol1(.4)

[eta = .4, f(eta) = 1.40000000000000034, diff(f(eta), eta) = 1.00000000000000022, diff(diff(f(eta), eta), eta) = -0.243774513041384287e-17, theta(eta) = .625958972186505536, diff(theta(eta), eta) = -.314549395236395634]

(7)

sol1(.5)

[eta = .5, f(eta) = 1.37026161183094430, diff(f(eta), eta) = .874752886901313142, diff(diff(f(eta), eta), eta) = .345911467377074400, theta(eta) = .432494259338694842, diff(theta(eta), eta) = -.382764248064397461]

(8)

sol1(.6)

[eta = .6, f(eta) = 1.36678221814533528, diff(f(eta), eta) = .771028661281065508, diff(diff(f(eta), eta), eta) = .407805382194403932, theta(eta) = .876413930517023876, diff(theta(eta), eta) = -.197648778495384870]

(9)

sol1(2)

[eta = 2., f(eta) = 2.66120522956795602, diff(f(eta), eta) = .991532161353848585, diff(diff(f(eta), eta), eta) = 0.251405465681268682e-1, theta(eta) = .635967939441598018, diff(theta(eta), eta) = -.144641270049362308]

(10)

``

``

``

 

``

``

NULL


 

Download staganation_point11.mw

 

in this program im trying to combine the result, but it showing some error can help me please

 

 

 

How can I sketch the angles such as Pi/6 , Pi/4 or 7Pi/6 and etc by Maple 18? In general, what is maple commend for sketching angles?

Hi

Here is my question: 

I have been trying to plot this on maple but i only get an ampry box

-arctan((2.m.x)/(1-x^2))

Simple as that... M is between 0 to 10

And x and y are supposed to fmbe from 0 to infinity(i dont know what to write for infinity so i give a large number like 2milions)

Am a post graduate student. am using maple as soft ware.. How can I create a 15 x 15 matrix in maple?

I have a data point set:

x_val:=<250,300,350,397,451,497,547,593,647,691,745,788,840,897>:
y_val:=<0,0.5,2,6.3,23.2,48.7,71.2,83.4,90.1,92.8,94.7,95.7,96.9,97.8>:

I want to make a least square fit using this difficult function:
 

function:=x->1-exp(-(k*exp(-(E/(8.314*873.15))*((873.15/x)-1)))*(0.026/350))

but both Statistics[Fit]:
 

with(Statistics):fit_nelog:=Fit(1-exp(-(k*exp(-(E/(8.314*873.15))*((873.15/x)-1)))*(0.026/350)),<x_val|y_val>,x,parameternames=[k,E],output=[parametervector,residualsumofsquares]);

and DirectSearch[DataFit]:

with(DirectSearch):fit_nelog2:=DataFit(1-exp(-(k*exp(-(E/(8.314*873.15))*((873.15/x)-1)))*(0.026/350)),x_val,y_val,x,method=cdos);


give wrong k,E parameters. The correct parameter values were obtained with Excel Solver:

k=27843.3551042397

E=68.4

The approximately correct parameters were fitted when using logarithm form of the function.
How can I obtain correct parameter values in Maple using given form of the function?

int(a(t)*b(t)+2*(diff(a(t), t))*(diff(b(t), t)), a(t));
Error, (in int) integration range or variable must be specified in the second argument, got a(t)
 

do not understand this error message,

how to integrate it?

Hello people in mapleprimes,
I have a question.
I appended two pictures where from the same code, two different orders of
expression appear.
How can I do for this so as not to get error messages?
The cause of this is simplify(%,symbolic) brings different order of term a__0^(-k)*F__D ahead of a parenthesis in a jpg.file and F__D*a__0^(-k) after
that parenthesis in another jpg.fine both in the line above that of  "dairihensu1."

In this case, What I can do?
Please help me.
Best wishes.

taro

my_code.mw

Original code is

e7_4:=F__D*(Omega+1)*beta/(beta-1) = F__I*a__D^(-k)*a__0^k+T^((sigma-k-1)/(-1+sigma))*F__D*phi^(k/(-1+sigma))+F__D;

a1:=beta=k/(sigma-1);
subs_free:=
  proc(a,b,c)
    local b1;
    b1:=isolate(b,c);
    subs(b1,a);
  end proc;
isolate(e7_4,a__D^(-k));simplify(%,symbolic);dairihensu1:=subs_free(%,a1,sigma);e7_5:=applyop(simplify,[2,4,1,3,2],dairihensu1);

A case without error.

A case with a error.

 

Hi everybody,

So today is 10-28-2016 and I explored Leyland Numbers for the first time, on Maple.  Please see my example file and let me know what your impression is.

x_to_the_yth_power_and_y_to_the_xth_power_take_4.mw

x_to_the_yth_power_and_y_to_the_xth_power_take_4.pdf

I have included a .pdf file so that the caual internet observer can also be aware of this information.

Regards,
Matt

 

Dear all,

I am creating an animation, and I was wondering if I can add a multiplier to my equation in a specific range.  So starting from z:=0.2 add a multiplier (z+1). The code I have so far is added. Does anyone know a code for this?

Kind regards

restart; 
with(plots); 
a := -1/2; b := 1/2; c := -2; d := 2; n := 20; 
g := proc (x) options operator, arrow; value(Int(sigma(t), t = 0 .. x)) end proc; 
sigma := proc (z) options operator, arrow; 2*sqrt(2*h^2-4*z^2)*z/h^2 end proc; 
h := i/n; 
for i to n do 
an2[i] := plot(sigma(z), z = -(1/2)*h .. (1/2)*h, view = [a .. b, c .. d], color = AQUAMARINE); 
an3[i] := plot(2*g(x), x = 0 .. (1/2)*h, view = [a .. b, c .. d], color = RED) 
end do; 
p := plots[display]([seq(an2[i], i = 1 .. n)], insequence = true); 
q := plots[display]([seq(an3[i], i = 1 .. n)], insequence = true); display(p, q)

 


 

restart; with(plots); beta := 0.1e-1; Bi := 1; Pr := 3.0; L0 := 1; w = 0.2e-1

Eq1 := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))-(diff(f(eta), eta))^2+beta*H(eta)*(F(eta)-(diff(f(eta), eta))) = 0

diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))-(diff(f(eta), eta))^2+0.1e-1*H(eta)*(F(eta)-(diff(f(eta), eta))) = 0

(1)

Eq2 := G(eta)*(diff(F(eta), eta))+F(eta)^2+beta*(F(eta)-(diff(f(eta), eta))) = 0

G(eta)*(diff(F(eta), eta))+F(eta)^2+0.1e-1*F(eta)-0.1e-1*(diff(f(eta), eta)) = 0

(2)

Eq3 := G(eta)*(diff(G(eta), eta))+beta*(f(eta)+G(eta)) = 0

G(eta)*(diff(G(eta), eta))+0.1e-1*f(eta)+0.1e-1*G(eta) = 0

(3)

Eq4 := H(eta)*F(eta)+H(eta)*(diff(G(eta), eta))+G(eta)*(diff(H(eta), eta)) = 0

H(eta)*F(eta)+H(eta)*(diff(G(eta), eta))+G(eta)*(diff(H(eta), eta)) = 0

(4)

Eq5 := (diff(theta(eta), eta, eta))/Pr+f(eta)*(diff(theta(eta), eta))+(2*beta*H(eta)*(1/3))*(theta[p](eta)-theta(eta)) = 0

.3333333333*(diff(diff(theta(eta), eta), eta))+f(eta)*(diff(theta(eta), eta))+0.6666666667e-2*H(eta)*(theta[p](eta)-theta(eta)) = 0

(5)

Eq6 := G(eta)*(diff(theta[p](eta), eta))+L0*beta*(theta[p](eta)-theta(eta)) = 0

G(eta)*(diff(theta[p](eta), eta))+0.1e-1*theta[p](eta)-0.1e-1*theta(eta) = 0

(6)

bcs1 := f(0) = 0, (D(f))(0) = 1, (D(theta))(0) = -Bi*(1-theta(0)), (D(f))(5) = 0, F(5) = 0, G(5) = -f(5), H(5) = w, theta(5) = 0, theta[p](5) = 0

f(0) = 0, (D(f))(0) = 1, (D(theta))(0) = -1+theta(0), (D(f))(5) = 0, F(5) = 0, G(5) = -f(5), H(5) = w, theta(5) = 0, theta[p](5) = 0

(7)

p := dsolve({Eq1, Eq2, Eq3, Eq4, Eq5, Eq6, bcs1}, numeric)

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

 

odeplot(p, [eta, f(eta)], 0 .. 10);

odeplot(p, [eta, f(eta)], 0 .. 10)

(8)

``

 

 


 

Download from_net.mw

Pressure_loss.mw

Hey all, could someone pls help me with how i can setup the equation for f in my worksheet. It should look like v and Rey with 45 data points. I've tried alot but i can't seem to solve it mysefl. Is it because i solve and map at the same time?

Thanks

I want to run a specific color red outside and yellow inside on my equation here using MAPLE 8.00:

plot3d([(0.5+cos(5*u))*sin(2*v),(0.5+cos(5*u))*cos(2*v),0.5*(cos(5*u)-0.3*cos(15*u)+0.02*cos(25*v))],u=0..2*Pi,v=0..2*Pi,axes=FRAMED);

is there some one here can help me? thanks...here is the example of color I want Eg: 

I got a set like this one:

g:={{3, 4, 5}, {5, 12, 13}, {6, 8, 10}, {7, 24, 25}, {8, 15, 17},

  {9, 12, 15}, {9, 40, 41}, {10, 24, 26}, {12, 16, 20},

  {12, 35, 37}, {14, 48, 50}, {15, 20, 25}, {15, 36, 39},

  {16, 30, 34}, {18, 24, 30}, {20, 21, 29}, {21, 28, 35},

  {24, 32, 40}, {27, 36, 45}}

And i want to count how many time a number apper in this set. Like an example:

3 appear 1

4 appear 1

5 appear 2.

Thank for your reading, :)

Sorry for my bad english :)

i wrote this problem to solve 

Delta= Sum(j=1 to n)SUM(i=j to n)(pi*hj/Ad(t,ij)*Et,ij))

Where n=70,  G= ftj (t)/(4+0.85*t) , where (t =8, 16, 24,…….up to 8*n), hj= 13 for all j except j1 =18

Ad= (Aj+s(mij-1)), where Aj varies

Mij=ES/E(G),          where E(G)= 57sqrt(1000*G)

 

n := 70;

70

(1)

i := seq(1 .. n, 1);

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70

(2)

t := proc (i) options operator, arrow; 8*i end proc;

proc (i) options operator, arrow; 8*i end proc

(3)

j := i;

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70

(4)

F = f(j);

F = f(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70)

(5)

F(1 .. 30) := 8;

8

(6)

F(31 .. 40) := 7;

7

(7)

F(41 .. 70) := 6;

6

(8)

G := proc (F, i) options operator, arrow; F*t/(4+.85*t) end proc;

proc (F, i) options operator, arrow; F*t/(4+.85*t) end proc

(9)

A := f(j);

f(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70)

(10)

A(1 .. 30) := 5184;

5184

(11)

A(31 .. 50) := 3600;

3600

(12)

A(51 .. 62) := 1936;

1936

(13)

A(63 .. 70) := 1024;

1024

(14)

s := f(j);

proc () option remember; table( [( 31 .. 50 ) = 3600, ( 63 .. 70 ) = 1024, ( 1 .. 30 ) = 5184, ( 51 .. 62 ) = 1936, ( 31 .. 40 ) = 3600 ] ) 'procname(args)' end proc

(15)

s(1 .. 10) := 128.0448;

128.0448

(16)

s(11 .. 20) := 63.763;

63.763

(17)

s(21 .. 30) := 79.92;

79.92

(18)

s(31 .. 40) := 64.08;

64.08

(19)

s(41 .. 50) := 47.88:

s(51 .. 62) := 31.944;

31.944

(20)

s(63 .. 70) := 12.49;

12.49

(21)

E := proc (G) options operator, arrow; 57*sqrt(1000*F) end proc;

proc (G) options operator, arrow; 57*sqrt(1000*F) end proc

(22)

Es := 29000;

29000

(23)

m := proc (E) options operator, arrow; Es/E(G) end proc;

proc (E) options operator, arrow; Es/E(G) end proc

(24)

Ad := proc (j, m) options operator, arrow; A+s*(m(E)-1) end proc;

proc (j, m) options operator, arrow; A+s*(m(E)-1) end proc

(25)

P := f(j);

proc () option remember; table( [( 21 .. 30 ) = 79.92, ( 31 .. 50 ) = 3600, ( 41 .. 50 ) = 47.88, ( 63 .. 70 ) = 12.49, ( 1 .. 30 ) = 5184, ( 51 .. 62 ) = 31.944, ( 11 .. 20 ) = 63.763, ( 31 .. 40 ) = 64.08, ( 1 .. 10 ) = 128.0448 ] ) 'procname(args)' end proc

(26)

P(1 .. 68) := 254.7;

254.7

(27)

P(69 .. 70) := 196.8;

196.8

(28)

h := f(j);

proc () option remember; table( [( 21 .. 30 ) = 79.92, ( 31 .. 50 ) = 3600, ( 41 .. 50 ) = 47.88, ( 63 .. 70 ) = 12.49, ( 1 .. 30 ) = 5184, ( 51 .. 62 ) = 31.944, ( 11 .. 20 ) = 63.763, ( 31 .. 40 ) = 64.08, ( 1 .. 10 ) = 128.0448 ] ) 'procname(args)' end proc

(29)

h(1) := 18;

18

(30)

h(2 .. 70) = 13;

h(2 .. 70) = 13

(31)

delta := sum(sum((P.h)/(E(G)*Ad)), i = 1 .. n, j = i)

Error, invalid input: sum uses a 2nd argument, k, which is missing

 

``


 

Download short.mw

First 1064 1065 1066 1067 1068 1069 1070 Last Page 1066 of 2250