maple2015

120 Reputation

6 Badges

4 years, 218 days

MaplePrimes Activity


These are questions asked by maple2015

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

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:

Dear Friends

I need to arrange legends from 2 to 1.2

Please hint me

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

 

 

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

Hi dear friends

How do I improve my ability to manipulate code?

Thanks

_________________________________________________________________________________________

 

restart:

Digits := 30: m := 20: p0 := 8.: g := -0.3: nu := 0.3: a := 1:


w := sum(b[n]*r^(1.02*n), n = 1 .. m):


ODE := r^3*(g*r+a)^3*(diff(w, r, r, r, r))+2*(diff(w, r, r, r))*r^2*(4*g*r+a)*(g*r+a)^2+(diff(w, r, r))*r*(3*g^3*nu*r^3+P*a^3*r^2+6*a*g^2*nu*r^2+11*g^3*r^3+3*a^2*g*nu*r+15*a*g^2*r^2+3*a^2*g*r-a^3)+(diff(w, r))*(6*g^3*nu*r^3+P*a^3*r^2+6*a*g^2*nu*r^2-2*g^3*r^3-3*a*g^2*r^2+a^3):


for P from p0 by 0.1e-1 to p0+1 do
W := simplify(subs(solve({seq(evalf(subs(r = j/(m+2), ODE)), j = 2 .. m-1), subs(r = 1, diff(w, r))}, {seq(b[j], j = 1 .. m)}), w/b[m])):
F := (int(((diff(W, `$`(r, 2))+(diff(W, r))/r)^2-(2*(1-nu))*(diff(W, `$`(r, 2)))*(diff(W, r))/r)*r*(1+g*r/a)^3, r = 0 .. a))/(int((diff(W, r))^2*r, r = 0 .. a)):
delta[trunc(100*(P-p0))] := abs(F-P) end do:


plot(W, r = 0 .. a);


Min := min(seq(delta[j], j = 0 .. 100)):


for j from 0 to 100 do if delta[j] = Min then print(p0+j/(100.)) else  end if end do;

What is the best source of learning maple for an abecedarian to become a professional?

First 6 7 8 9 10 11 Page 8 of 11