Question: Numerical Solving of System

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));
Please Wait...