## 155 Reputation

11 years, 205 days

## What wrong with the codes Different answ...

Maple

Hai, any one can help me to rectify my error

I have evaluate the codes in maple and got different answer.

restart:
epsilon:=0.2:z:=0.9:m:=10:k:=0.8:
A:=(-m^2/4)-(1/(4*k)):
h(z):=1+epsilon*sin(2*Pi*z):
S1:=(h(z)^2)/(4*A)-ln(A*h(z)^2+1)*(1+h(z)^2)/(4*A):
g1:=evalf((1/S1));

g1 := 8.821345336-11.12386331*I

but in mathematica i am getting  answer as (27.8647 + 35.2042 I)

ClearAll;
\[Epsilon] = 0.2;
h[z_] = 1 + \[Epsilon]*Sin[2*Pi*z];
m = 10;
k = 0.8;
A = -(m^2/4) - (1/4 k);
S1[z_] = (h[z]^2/4*A) - ((1 + h[z]^2)/4*A)*
Log[(A*h[z]^2) + 1] /. {z -> 0.9}

## How to obtain the graph ?...

Maple 18

Hellow, help required to remove the errors

How to  obtain the pressure drop i am unable to get the output. I am uploading the file  and the equations

help_dp.mw

## How to get the desire output for the cod...

Maple 18

Help required to get the derire output  for the differential equation . I am attacing the codes and the sample output.

restart:
with(DETools):
a1:=m^2*((1-N)/(2-N))*(r/2)*Dp:
a2:=((1-N)*m^2/(2-N))*(((Nb-Nt)/64)*(r^5/6-h^2*r/2)-(B[r]/4)*(r^3/4-h^2*r/2)*(Nt/Nb)):
a3:=a1-a2:
a4:=r^2*a3:
DE1:=r^2*diff(v(r),r,r)+r*diff(v(r),r)-(m^2*r^2+1)*v(r)=a4:
b1:=dsolve(DE1,v(r)):

## Evaluate and print values...

Maple 13

how to evaulate the value of R(z) for different values of z=0 to 1 with an interval of 0.1 and print ten values  in one column

R(z):= 1-cos^2*(Pi*z);

## display graphs side by side with multipl...

Greetings, any one help  me to remove the error  in plots display,here is the codes

i want 2 rows and 2 columns graphs with 9 lines in each graphs.

i am attaching the maple questions also  like this i want

https://www.mapleprimes.com/questions/224443-How-Do-I-Get-A-Legend-Only-Three-Curves--In-Maple

https://www.mapleprimes.com/questions/225078-How-To-Remove-Error-In-Plot-Display

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

alpha:=Pi/6:
with(plots):
display
( Vector[row]
([seq [ seq
( plot
( [ seq
( eval(lambda(Gamma2), Nb=j),
j in [0.1,0.2,0.3]
)
d1 in [0.1,0.2,0.5])] ],
delta2=0.02..0.1,
legend=[Nb=0.1,Nb=0.2,Nb=0.3],
labels=[typeset(`&delta;1`), typeset(conjugate(`&Delta;p`))],
title=typeset("Effect of ", ''alpha'', " when ", Gamma,"2=", Gamma2)
),
Gamma2 in [10,20,30,40]
)
]
)
);

 3 4 5 6 7 8 9 Last Page 5 of 15
﻿