waseem

65 Reputation

3 Badges

6 years, 154 days

MaplePrimes Activity


These are questions asked by waseem

Hi, all i am unable to plot the graphs ,can any one help me to overcome the error in plotting the graphs.I am using the maple 13. I am attaching the codes

restart:
with(plots):
with(IntegrationTools):
d1:=0.2:L1:=0.2:L2:=0.2:B1:=0.7:B:=1:beta:=0.01:
d2:=0.6:m:=0.1:k:=0.1: 

h:=z->piecewise( z<=d1,    1,
                z<=d1+L1,   1-(gamma1/(2))*(1 + cos(2*(Pi/L1)*(z-d1-L1/2))), 
                z<=B1-L2/2,  1 ,          
                z<=B1,  1-(gamma2/(2))*(1 + cos(2*(Pi/L2)*(z - B1))),
                z<=B1+L2/2,  1-(gamma2/(2))*(1 + cos(2*(Pi/L2)*(z - B1))),
                 z<=B,    1):
                
A:=(-m^2/4)-(1/(4*k)):
S1:=(h(z)^2)/(4*A)-ln(A*h(z)^2+1)*(1+h(z)^2)/(4*A):
b1:=evalf((1/S1)):               
c1:=evalf(Int(b1,z=0..1)):

plot([seq(eval(c1,gamma2=j),j in[0,0.02,0.06])],gamma1=0.02..0.1,legend = [gamma2 = 0.0, gamma2 = 0.02,gamma2 =0.04],linestyle = [solid,dash,dot],color = [black, black,black],axes=boxed); 
 

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}

 

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

 

 

 

 


 

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)):

 

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);
 

1 2 3 4 5 6 7 Last Page 1 of 11