Advanced Search

Question: Ploting real and complex part of a function

Question:Ploting real and complex part of a function

gaurav_rs 65 Maple 2017
complex plot
December 05 2018
0

Dear users,

I have an issue with finding real part of a complex variable function. In calculating the real part I see two arguments and the plot is not smooth. How to get real part correct. The worksheet is attached.
 

``

 

 
> 

##Toya complex variable method

> 

``

> 

restart;

> 

stress_c:=-(1+1/nu_c)*k*p2*zeta_c/2;

-(1/2)*(1+1/nu_c)*k*p2*zeta_c

(1.1)
> 

p2:=(c0_c-d_1c/k)*(z-a*(cos(alpha)+2*lambda*sin(alpha)))+(1-k)/k*a*(N_infty-T_infty)*exp(2*I*phi_c+2*lambda*(alpha-Pi))*((a*(cos(alpha)-2*lambda*sin(alpha)))/z-a^2/z^2)

(c0_c-d_1c/k)*(z-a*(cos(alpha)+2*lambda*sin(alpha)))+(1-k)*a*(N_infty-T_infty)*exp((2*I)*phi_c+2*lambda*(alpha-Pi))*(a*(cos(alpha)-2*lambda*sin(alpha))/z-a^2/z^2)/k

 (1.2)
> 

``

> 

z := exp(I*theta)

exp(I*theta)

(1.3)
> 

``

> 

k := beta_c/(1+nu_c)

beta_c/(1+nu_c)

(1.4)
> 

nu_c := (kappa2*mu+mu2)/(kappa*mu2+mu)

(kappa2*mu+mu2)/(kappa*mu2+mu)

(1.5)
> 

d_1c := (N_infty+T_infty)*(1/2)

(1/2)*N_infty+(1/2)*T_infty

(1.6)
> 

lambda := -evalf(ln(nu_c)/(2*Pi))

-.1591549430*ln((kappa2*mu+mu2)/(kappa*mu2+mu))

(1.7)
> 

``

> 

beta_c := mu*(1+kappa2)/(kappa*mu2+mu)

mu*(1+kappa2)/(kappa*mu2+mu)

(1.8)
> 

zeta_c := ((z-a*exp(I*alpha))/(z-a*exp(-I*alpha)))^(I*lambda)/((z-a*exp(I*alpha))^.5*(z-a*exp(-I*alpha))^.5)

((exp(I*theta)-a*exp(I*alpha))/(exp(I*theta)-a*exp(-I*alpha)))^(-(.1591549430*I)*ln((kappa2*mu+mu2)/(kappa*mu2+mu)))/((exp(I*theta)-a*exp(I*alpha))^.5*(exp(I*theta)-a*exp(-I*alpha))^.5)

(1.9)
> 

``

> 

c0_c := G_c+I*H_c

G_c+I*H_c

(1.10)
> 

G_c:=(0.5*(T_infty+N_infty)*(1-(cos(alpha)+2*lambda*sin(alpha))*exp(2*lambda*(evalf(Pi)-alpha)))-0.5*(1-k)*(1+4*lambda^2)*(N_infty-T_infty)*(sin(alpha))^2*cos(2*phi_c))/(2-k-k*(cos(alpha)+2*lambda*sin(alpha))*exp(evalf(2*lambda*(Pi-alpha))));

(.5*(N_infty+T_infty)*(1-(cos(alpha)-.3183098860*ln((kappa2*mu+mu2)/(kappa*mu2+mu))*sin(alpha))*exp(-.3183098860*ln((kappa2*mu+mu2)/(kappa*mu2+mu))*(3.141592654-alpha)))-.5*(1-mu*(1+kappa2)/((kappa*mu2+mu)*(1+(kappa2*mu+mu2)/(kappa*mu2+mu))))*(.1013211835*ln((kappa2*mu+mu2)/(kappa*mu2+mu))^2+1)*(N_infty-T_infty)*sin(alpha)^2*cos(2*phi_c))/(2-mu*(1+kappa2)/((kappa*mu2+mu)*(1+(kappa2*mu+mu2)/(kappa*mu2+mu)))-mu*(1+kappa2)*(cos(alpha)-.3183098860*ln((kappa2*mu+mu2)/(kappa*mu2+mu))*sin(alpha))*exp(-.3183098860*ln((kappa2*mu+mu2)/(kappa*mu2+mu))*(3.141592654-1.*alpha))/((kappa*mu2+mu)*(1+(kappa2*mu+mu2)/(kappa*mu2+mu))))

 (1.11)
> 

H_c:=0.5*(1-k)*(1+4*lambda^2)*(-T_infty+N_infty)*(sin(alpha))^2*sin(2*phi_c)/(k*(1+(cos(alpha)+2*lambda*sin(alpha))*exp(2*lambda*(evalf(Pi)-alpha))));

.5*(1-mu*(1+kappa2)/((kappa*mu2+mu)*(1+(kappa2*mu+mu2)/(kappa*mu2+mu))))*(.1013211835*ln((kappa2*mu+mu2)/(kappa*mu2+mu))^2+1)*(N_infty-T_infty)*sin(alpha)^2*sin(2*phi_c)*(kappa*mu2+mu)*(1+(kappa2*mu+mu2)/(kappa*mu2+mu))/(mu*(1+kappa2)*(1+(cos(alpha)-.3183098860*ln((kappa2*mu+mu2)/(kappa*mu2+mu))*sin(alpha))*exp(-.3183098860*ln((kappa2*mu+mu2)/(kappa*mu2+mu))*(3.141592654-alpha))))

 (1.12)
> 

##Input

> 

alpha:=evalf(Pi/6)

.5235987758

(1.13)
> 

phi_c:=alpha;

.5235987758

(1.14)
> 

N_infty:=0;

0

(1.15)
> 

T_infty:=1;

1

(1.16)
> 

a:=1;nu2:=22/100;kappa2:=3-4*nu2;nu:=35/100;kappa:=3-4*nu;mu:=239/100;mu2:=442/10;

1

 

11/50

 

53/25

 

7/20

 

8/5

 

239/100

 

221/5

(1.17)
> 

``

> 

stress_c

-(9321/123167)*(((.5586916801-.5*(.8660254037-.1591549431*ln(123167/182775))*exp(-.8333333329*ln(123167/182775))+0.5946710490e-2*ln(123167/182775)^2)/(22817/11767-(717/11767)*(.8660254037-.1591549431*ln(123167/182775))*exp(-.8333333329*ln(123167/182775)))-(1.668336947*I)*(.1013211835*ln(123167/182775)^2+1)/(1+(.8660254037-.1591549431*ln(123167/182775))*exp(-.8333333329*ln(123167/182775)))-11767/1434)*(exp(I*theta)-.8660254037+.1591549431*ln(123167/182775))-(11050/717)*exp(1.047197552*I+.8333333328*ln(123167/182775))*((.8660254037+.1591549431*ln(123167/182775))/exp(I*theta)-1/(exp(I*theta))^2))*((exp(I*theta)+(-.8660254037-.5000000002*I))/(exp(I*theta)+(-.8660254037+.5000000002*I)))^(-(.1591549430*I)*ln(123167/182775))/((exp(I*theta)+(-.8660254037-.5000000002*I))^.5*(exp(I*theta)+(-.8660254037+.5000000002*I))^.5)

 (1.18)
> 

assume((1/6)*Pi < theta, theta < 2*Pi-(1/6)*Pi)

> 

simplify(evalc(Re(stress_c)))

-0.8815855810e-10*((((1.000000000*cos(theta)^7+(0.5294827753e-2+.5671599115*sin(theta))*cos(theta)^6-4.533186669*cos(theta)^5+(-11.80630620+4.886343937*sin(theta))*cos(theta)^4+3.402782742*cos(theta)^3+(9213180122.+0.9866808100e-1*sin(theta))*cos(theta)^2+(-0.1055437876e11+0.1595769608e11*sin(theta))*cos(theta)-5794103792.*sin(theta)+1760041721.)*cos(.5*arctan(sin(theta)+.5000000002, cos(theta)-.8660254037))+(-.5600908440*cos(theta)^7+(0.6523625301e-2+1.134319823*sin(theta))*cos(theta)^6+4.644568297*cos(theta)^5+(-0.2905669688e-1+10.20004207*sin(theta))*cos(theta)^4-0.1774243515e-1*cos(theta)^3+(0.1595769609e11-9.082306669*sin(theta))*cos(theta)^2+(-7023191163.-9213180109.*sin(theta))*cos(theta)-3154310102.*sin(theta)-7408031461.)*sin(.5*arctan(sin(theta)+.5000000002, cos(theta)-.8660254037)))*cos(.5*arctan(sin(theta)-.5000000002, cos(theta)-.8660254037))+(-.5600908440*cos(theta)^7+(1.134319823*sin(theta)+0.4756356038e-2)*cos(theta)^6+4.644568284*cos(theta)^5+(11.37920491*sin(theta)-0.2640575516e-1)*cos(theta)^4-0.1774243890e-1*cos(theta)^3+(-11.39571957*sin(theta)+0.1595769607e11)*cos(theta)^2+(-9213180108.*sin(theta)-7023191160.)*cos(theta)-7408031458.-3154310086.*sin(theta))*sin(.5*arctan(sin(theta)-.5000000002, cos(theta)-.8660254037))*cos(.5*arctan(sin(theta)+.5000000002, cos(theta)-.8660254037))+(-1.000000000*cos(theta)^7+(-.5671599115*sin(theta)-0.5294826902e-2)*cos(theta)^6+4.531921682*cos(theta)^5+(-4.886343941*sin(theta)+11.76153292)*cos(theta)^4-3.358186195*cos(theta)^3+(-0.9866807692e-1*sin(theta)-9213180122.)*cos(theta)^2+(-0.1595769609e11*sin(theta)+0.1055437877e11)*cos(theta)-1760041726.+5794103798.*sin(theta))*sin(.5*arctan(sin(theta)+.5000000002, cos(theta)-.8660254037))*sin(.5*arctan(sin(theta)-.5000000002, cos(theta)-.8660254037)))*cos(0.314104002e-1*ln(1492820323.-1292820323.*cos(theta)+746410161.*sin(theta))-0.314104002e-1*ln(-1292820322.*cos(theta)-746410161.4*sin(theta)+1492820322.))+(((-.5600908440*cos(theta)^7+(1.134319823*sin(theta)+0.4756356038e-2)*cos(theta)^6+4.626658979*cos(theta)^5+(-0.2905667760e-1+10.24488508*sin(theta))*cos(theta)^4-.1341529536*cos(theta)^3+(0.1595769608e11-9.127079936*sin(theta))*cos(theta)^2+(-7023191161.-9213180109.*sin(theta))*cos(theta)-3154310089.*sin(theta)-7408031435.)*cos(.5*arctan(sin(theta)+.5000000002, cos(theta)-.8660254037))+(-1.134319823*cos(theta)^7-.5671599115*sin(theta)*cos(theta)^6+4.531921682*cos(theta)^5+(11.80860365-4.107288978*sin(theta))*cos(theta)^4-3.402959469*cos(theta)^3+(-9213180123.+0.1774243833e-1*sin(theta))*cos(theta)^2+(0.1055437876e11-0.1595769608e11*sin(theta))*cos(theta)+5794103807.*sin(theta)-1760041748.)*sin(.5*arctan(sin(theta)+.5000000002, cos(theta)-.8660254037)))*cos(.5*arctan(sin(theta)-.5000000002, cos(theta)-.8660254037))+(-1.000000000*cos(theta)^7-.5671599115*sin(theta)*cos(theta)^6+4.537223485*cos(theta)^5+(-4.886343950*sin(theta)+11.80860366)*cos(theta)^4-3.358186195*cos(theta)^3+(-0.9866807250e-1*sin(theta)-9213180123.)*cos(theta)^2+(0.1055437876e11-0.1595769608e11*sin(theta))*cos(theta)-1760041739.+5794103821.*sin(theta))*sin(.5*arctan(sin(theta)-.5000000002, cos(theta)-.8660254037))*cos(.5*arctan(sin(theta)+.5000000002, cos(theta)-.8660254037))+(.5600908440*cos(theta)^7+(-1.134319823*sin(theta)-0.4756356038e-2)*cos(theta)^6-4.644554360*cos(theta)^5+(-10.21771474*sin(theta)+0.2905668928e-1)*cos(theta)^4+0.1774243685e-1*cos(theta)^3+(9.082306650*sin(theta)-0.1595769608e11)*cos(theta)^2+(9213180109.*sin(theta)+7023191165.)*cos(theta)+7408031453.+3154310085.*sin(theta))*sin(.5*arctan(sin(theta)+.5000000002, cos(theta)-.8660254037))*sin(.5*arctan(sin(theta)-.5000000002, cos(theta)-.8660254037)))*sin(0.314104002e-1*ln(1492820323.-1292820323.*cos(theta)+746410161.*sin(theta))-0.314104002e-1*ln(-1292820322.*cos(theta)-746410161.4*sin(theta)+1492820322.)))/((-sin(theta)+2.-1.732050807*cos(theta))^(1/4)*(sin(theta)+2.-1.732050807*cos(theta))^(1/4))

 (1.19)
> 

plot(%, theta = (1/6)*Pi .. 2*Pi-(1/6)*Pi)

 


 

Download Toya_complexPlot2.mw
Please Wait...