## 120 Reputation

5 years, 175 days

## Solving nonlinear equations...

Hi

I need amount of P.

## Plots in tripartite coordinate system...

Maple

Hi

How  can we draw diagrams in tripartite coordinate system?

## Numerical Amounts of Definite Integrals...

------------------------------------------------------------------------------------------------------------------------------

restart:

Digits :=30: m := 4: g :=0.3: nu := 0.4: a := 1:k:=0.2:

W:=cos(n*2*evalf(Pi)*( r-k*a)/((1-k)*a)):

ODE:=((2*(diff(W, r, r, r, r))+2*(diff(W, r, r, r))/r-4*(diff(W, r, r))/r^2+4*(diff(W, r))/r^3-(2*(1-nu))*(diff(W, r, r, r))/r+(4*(1-nu))*(diff(W, r, r))/r^2-(4*(1-nu))*(diff(W, r))/r^3)*r*(1+g*r/a)^3+(2*(2*(diff(W, r, r, r))+2*(diff(W, r, r))/r-2*(diff(W, r))/r^2-(2*(1-nu))*(diff(W, r, r))/r+(2*(1-nu))*(diff(W, r))/r^2))*(1+g*r/a)^3+(6*(2*(diff(W, r, r, r))+2*(diff(W, r, r))/r-2*(diff(W, r))/r^2-(2*(1-nu))*(diff(W, r, r))/r+(2*(1-nu))*(diff(W, r))/r^2))*r*(1+g*r/a)^2*g/a+(6*(2*(diff(W, r, r))+2*(diff(W, r))/r-(2*(1-nu))*(diff(W, r))/r))*(1+g*r/a)^2*g/a+(6*(2*(diff(W, r, r))+2*(diff(W, r))/r-(2*(1-nu))*(diff(W, r))/r))*r*(1+g*r/a)*g^2/a^2-((2*(diff(W, r, r, r)+(diff(W, r, r))/r-(diff(W, r))/r^2))/r-(2*(diff(W, r, r)+(diff(W, r))/r))/r^2-(2*(1-nu))*(diff(W, r, r, r))/r+(2*(1-nu))*(diff(W, r, r))/r^2)*r*(1+g*r/a)^3-((2*(diff(W, r, r)+(diff(W, r))/r))/r-(2*(1-nu))*(diff(W, r, r))/r)*(1+g*r/a)^3-(3*((2*(diff(W, r, r)+(diff(W, r))/r))/r-(2*(1-nu))*(diff(W, r, r))/r))*r*(1+g*r/a)^2*g/a)/(2*(diff(W, r, r))*r+2*(diff(W, r))) :

for i from 1 to m do  for j from 1 to m do
C[i,j]:=int(subs(n=i,ODE)*subs(n=j,W),r=k*a..a)
end do:end do:

## Adjusting the Order of Legends ...

Maple

Dear Friends

I need to arrange legends from 2 to 1.2

Thanks

restart:
e := x*t*(E[c]-E[m])/((2*(x+2))*(x*E[c]+E[m])):
EZ := (E[m]-E[c])*((t-2*z-2*e)/(2*t))^x+E[c]:
assume(x > 0):
d := int(EZ*z^2/(-nu^2+1), z = -(1/2)*t-e .. (1/2)*t-e):
R := simplify(subs(E[c] = k*E[m], d/subs(x = 0, d))):
LS[1.2] := solid:LS[1.4] := longdash:LS[1.6] := dash:LS[1.8] := spacedash:LS[2.0] := dashdot:
for k from 1.2 by .2 to 2.0 do
p[k] := plot(R, x = 0 .. 5, color = black, axes = framed, font = [Times, 13], size = [650, 550], resolution = 1200, legendstyle = [location = right], legend = ['E[c]/E[m]' = k], linestyle = LS[k], thickness = 2, labels = ['n', 'P[FG]/P[Metallic]'], labelfont = [Times, 14], tickmarks = [10, 15])
end do:
with(plots):
display({`\$`(p[1+.2*m], m = 1 .. 5)});

## Numerical Solving of System...

Maple
`Hi dear friends Is there an interactive package management utility or a way for solving following problem?for m>=4 It dosent work!restart:Digits :=30: m := 3: g :=0.3: nu := 0.2: a := 1:w := sum(b[n]*cos(n*r), n = 1 .. m):W := simplify( subs( solve( { subs(r = 1, diff(w, r\$2)+nu*diff(w,r))},{ b[1] }),w)):d1:=diff(W,r):d2:=diff(d1,r):F:= int( ((d2+d1/r)^2-(2*(1-nu))*d2*d1/r)*r*(1+g*r/a)^3,r = 0 .. a) /int( d1^2*r,r = 0 .. a):FW := simplify( subs(solve( {seq(subs(r=n/m,diff(F,b[n])),n=2..m)},{ seq(b[j], j = 2 .. m) }),F));`
 First 6 7 8 9 10 11 12 Page 8 of 12
﻿