tomleslie

4893 Reputation

15 Badges

9 years, 183 days

MaplePrimes Activity


These are replies submitted by tomleslie

@samiyare 

First of all, I probably should have spotted hat you were multiplying all values in the variable 't_number' by 60 - but I missed this, sorry. The effect of this "scaling" of the time axis, appears to just scale everything else - so one has to adjust ranges over whihc plots are made, search range for the optimizer, etc, etc. Note that both adjsuting ranges and controlling the optimizer range are used essentially to speed-up execution time: if you were pepared to wait quite a long time, then such range restrictions *might* be unnecessary - the relevant minimum might still be obtained, it would just take a loooonnngg time.

For what it is worth, the attached shows a minor variation of my original code implementing the 60*t_number with corresponding changes in plot ranges and optimization search space. This still locates the new "answers" in reasonable execution time


 

  restart;

#
# Parameters and procedures
#
  t_number:=<0, 10, 20, 30, 40>:
  m_number:=<11.50000000, 4.641182547, 1.273311993, 0.361845238, 0.288711649>:
  L:=2:
  getSol:= proc(d)
                local Mx0:= 0.05, cx0:= Mx0/(1-Mx0),
                      Mdb_i:= m_number[1], ct0:= Mdb_i,
                      PDE:= diff(C(x,t),t)=d*diff(C(x,t),x,x),
                      IBC:= {C(x,0)=ct0, C(0,t)=cx0, D[1](C)(L,t)=0};
                return pdsolve(PDE, IBC, numeric);
          end proc :
  Q:= proc(d)
           local pds, i, S:=0:
           if   not type(d, numeric)
           then return 'procname(_passed)'
           end if;
           pds:=getSol(d);
           for i from 1 to 5 do
               pds:-value( t=t_number[i]*60,
                           output=listprocedure
                         );                  # get the solution at the desired time
               eval(C(x,t), %);              # extract the C(x,t) at that time
               int(%, 0..2, numeric) / L;    # compute the average of C(x,t) at that time
               S += (% - m_number[i])^2;     # accumulate the residuals
           end do;
           return sqrt(S);
      end proc:

#
# Rouben's calculation
#
# A few representative values - shows that the
# "optimum" value for the parameter 'd' is
# somewhere around 0.14
#
# Change the following by ~10x to reflect the
# scaling factor
#
#  Q(0.1), Q(0.2), Q(0.5), Q(1.0);
#
  Q(0.001);Q(0.002);Q(0.003);Q(0.004);
  plot(Q(d), d=0.001 .. 0.004, numpoints=5, view=0..2);

HFloat(4.979783246316887)

 

HFloat(1.1695635853944253)

 

HFloat(1.5511287221444159)

 

HFloat(2.6009866748705806)

 

 

#
# Use the Optimization() package to get the value
# of diffusion coefficient which produces the
# minimum overall residual. Compute the solution
# at this "optimal" diffusion coefficient
#
  oSol:= Optimization:-Minimize(Q, 0..0.005);
  pdsOpt:= getSol( oSol[2][1] ):
#
# Using the optimal diffusion coefficient compute
# the average residual at each of the time-values
# in the vector t_number. These values *ought* to
# match somewhat with the corresponding value in
# the vector m_number
#
  resid:= Vector( 5,
                  i-> int
                      ( eval
                        ( C(x,t),
                          pdsOpt:-value
                                  ( t=t_number[i],
                                    output=listprocedure
                                  )
                        ),
                        0..2,
                        numeric
                      )/L
                 );
#
# Get the mean square error in the average C(x,t)-value
# at each of the time-values in the vector t_number
#
  `^`~((resid-m_number),2);

[.8889441947641348, Vector(1, {(1) = 0.2306692069320807e-2})]

 

Vector(5, {(1) = 11.499999999999842, (2) = 10.51283349395565, (3) = 10.108156818332095, (4) = 9.797236955370293, (5) = 9.534923705008406})

 

Vector[column](%id = 18446744074566322166)

(1)

#
# Plot C(x,t) for the optimal value of the diffusion
# coefficient for each of the time-values in t_number
#
  cols:= [red, green, blue, brown, black]:
  plots:-display( [ seq
                    ( pdsOpt:-plot
                      ( t=t_number[j],
                        x=0..L,
                        color=cols[j]
                      ),
                      j=1..numelems(t_number)
                    )
                  ]
                );
#
# Plot the "general 3-D" solution
#
  pdsOpt:-plot3d( x=0..L, t=0..t_number[-1]);

 

 

 


 

Download diffCoeff2.mw

you will have to upload the offending worksheet using the big green up-arrow in the Mapleprimes toolbar.

The TestDimensions() command i used by a lot of commands in the Units() system to verify that the units associated with terms in an expression are "compatible" - in other words making sure that on doesn't add "seconds" to "kilograms". A quick check shows that  it is capable of checking differetn ways of expressing the "same" unit, as shown in the attached "toy" example

restart;
kernelopts(version);
with(Units):
TestDimensions(2.114163508*10^7*Unit(kg/s^2)+2.114163508*10^7*Unit(J/m^2));

`Maple 2019.1, X86 64 WINDOWS, May 21 2019, Build ID 1399874`

 

Automatically loading the Units[Simple] subpackage

 

true

(1)

 


 

Download testUnit.mw

@mehran rajabi 

to the problem I answered in your other post here

https://www.mapleprimes.com/questions/227667-Solve-PDE-With-RungeKutta-Method-And-Plot

Posting the same question multiple times just confuses everyone - don't do it

@maple2015 

if anything- is the problem????

Your PDE is fourth-order in the variable 'x', and therefore needs 4 'x-related' boundary conditions - you have five

Your PDE is second-order in the variable 't', and therefore needs 2 't-related' initial conditions - you have one

The right hand side of this PDE (the piecewise function) looks a bit weird - it has a value at x=L, and absoutely nowhere else. This is probably going to result in a serious discontinuity an any solution (and may make the solution impossible to generate). Is this piecewise definition correct???

@Carl Love 

I have been reading this thread with increasing bewilderment. This confusion exists on two levels

  1. The worksheet  RLC.mw you posted does not generate the displayed answers! When I re-execute it (see attached below) the values of 'Sol' and related plot are different: note the arguments of the sin() and cos() terms should be 3*t, not 7*t. How is this possible? It is relatively easy to demonstrate that this coefficient (aka the the damped frequency) ought to be 3, although this demonstration is further confused by item 2 below
  2. As written, the ODE 'RLC' is incredibly confusing. You are using the symbol 'R' to denote inductance and the symbol 'L' to denote resistance, which is perverse to say the least. However it is the only way that applying Kirchoff's current law to to a series circuit of this type would make any sense. Of course one can use any names one wants, but OP's original post has R = 20 , C = 0.01 , L = 10, whereas  has   R=10, L=20, C= 1/100 so names have been swapped in the equation RLC and values have been swapped back in the params list. Hence the actual effect is nil (other than to confuse old people like me)
  3. In the second attachment below, I have defined the series RLC circuit in the conventional way, and swapped the parameter values around in order to achieve some kind of consistency. This means that both of these worksheets will now generate the same answers (but neither will generate the damped frequency value referred to in (1) above!)

restart:

RLC:= R*diff(x(t),t$2) + L*diff(x(t),t) + x(t)/C = 0;

R*(diff(diff(x(t), t), t))+L*(diff(x(t), t))+x(t)/C = 0

(1)

IC:= x(0)=10, D(x)(0)=0;

x(0) = 10, (D(x))(0) = 0

(2)

params:= [R=10, L=20, C= 1/100]:

Sol:= dsolve(eval({RLC, IC}, params));

x(t) = (10/3)*exp(-t)*sin(3*t)+10*exp(-t)*cos(3*t)

(3)

plot(eval(x(t), Sol), t= 0..5);

 

#Original square wave function:
step:= t-> piecewise(t < 0, 0, 1):
f:= t-> sum((-1)^n*step(t-n), n= 0..50):

#
#It's equivalent to this simpler function:
f:= t-> piecewise(floor(t)::even, 1, 0):

plot(f, 0..5, axes= frame, title= "Square Wave",thickness= 4, gridlines= false);

 

RLC2:= lhs(RLC) = 200*f(t/3);

RLC2 := R*(diff(x(t), t, t))+L*(diff(x(t), t))+x(t)/C = 200*piecewise((floor((1/3)*t))::even, 1, 0)

(4)

Sol2:= dsolve(eval({RLC2, IC}, params), numeric):

plots:-odeplot(Sol2, [[x(t), diff(x(t),t)]], t= 0..25, numpoints= 1000, thickness= 2, title= "Phase portrait", gridlines= false);

 

plots:-odeplot(Sol2, [t, x(t)], t= 0..25, numpoints= 1000, thickness= 2, title= "Solution curve", gridlines= false);

 

 

Download carlsCode.mw

restart:

#
# swapped round R and L
#
RLC:= L*diff(x(t),t$2) + R*diff(x(t),t) + x(t)/C = 0;

L*(diff(diff(x(t), t), t))+R*(diff(x(t), t))+x(t)/C = 0

(1)

IC:= x(0)=10, D(x)(0)=0;

x(0) = 10, (D(x))(0) = 0

(2)

#
# swapped rond values for R and L
#
params:= [R=20, L=10, C= 1/100]:

Sol:= dsolve(eval({RLC, IC}, params));

x(t) = (10/3)*exp(-t)*sin(3*t)+10*exp(-t)*cos(3*t)

(3)

plot(eval(x(t), Sol), t= 0..5);

 

#Original square wave function:
step:= t-> piecewise(t < 0, 0, 1):
f:= t-> sum((-1)^n*step(t-n), n= 0..50):

#
#It's equivalent to this simpler function:
f:= t-> piecewise(floor(t)::even, 1, 0):

plot(f, 0..5, axes= frame, title= "Square Wave",thickness= 4, gridlines= false);

 

RLC2:= lhs(RLC) = 200*f(t/3);

RLC2 := L*(diff(x(t), t, t))+R*(diff(x(t), t))+x(t)/C = 200*piecewise((floor((1/3)*t))::even, 1, 0)

(4)

Sol2:= dsolve(eval({RLC2, IC}, params), numeric):

plots:-odeplot(Sol2, [[x(t), diff(x(t),t)]], t= 0..25, numpoints= 1000, thickness= 2, title= "Phase portrait", gridlines= false);

 

plots:-odeplot(Sol2, [t, x(t)], t= 0..25, numpoints= 1000, thickness= 2, title= "Solution curve", gridlines= false);

 

 


 

Download convRLC.mw

@nm 

because we "froze" two different things. In your code

r    := (x^2-1)/x^2:
Z    := freeze(r):

whereas in mine

Z:=freeze( (x^2-1) ):

 

 

 

@nm 

one of the following more "suitable".

A "trick" I have observed with algsubs() is to (sometimes) do the substitution in a couple of stages


 

  restart;
  expr1:=(x^2-1)*y/x^2+(x^2-1)/x^2;

(x^2-1)*y/x^2+(x^2-1)/x^2

(1)

#
# The 'easy' way
#
  thaw(algsubs((y+1)=freeze((y+1)), expr1));

(y+1)*(x^2-1)/x^2

(2)

#
# Slightly harder - using algsubs and freeze/thaw
#
  Z:=freeze( (x^2-1) ):
  simplify
  ( collect
    ( algsubs
      ( x^2-1=Z,
        expr1
      ),
      Z
    )
  ):
  thaw(%);

(y+1)*(x^2-1)/x^2

(3)

#
# Even harder - using algsubs witout freeze/thaw
#
  subs
  ( Z1=(x^2-1),
    subs
    ( Z2=Z1/x^2,
      collect
      ( algsubs
        ( Z1/x^2=Z2,
          algsubs
          ( (x^2-1)=Z1,
            expr1
          )
        ),
        Z2
      )
    )
  );

(y+1)*(x^2-1)/x^2

(4)

 

 


 

Download frth3.mw

@minhthien2016 

to my previous worksheet to return the radius of the circumcircle and the incircle

triFun2.mw

odetest() only makes very careful simplifications when checking the validity of solutions.

It is sometimes(?) informative to perform the equivalent of odetest() manually - just substitute the solution into the ode, and make your own "simplifications", although I would advise doing this very carefully

For your specific example, the first command in the following evaluates to 'true', while the second evaluates to 'false' - the 'symbolic' option in the simplify() command makes all the difference

evalb(simplify(subs( my_sol, ode), symbolic));
evalb(simplify(subs( my_sol, ode)));

 

@nm
using patmatch! Whilst it is relatively simple to modify my previous response to accommodate your latest counterexample (see the attached), I doubt that this approach could be extended to cover absolutely arbitrary initial  expressions

restart;
expr1:=sin(x)*(x*y)^(1/3);
expr2:=sin(x)+(x*y)^(1/6);
expr3:=sin(x)*(x*y)^(1/3)+3;
f:=ex->`if`(type(ex,`+`),g~([op( ex)]), [op(ex)]);
g:=ex->`if`(type(ex,`*`), op(ex), ex);
getMatch:=expr-> `if`
                 ( ormap
                   ( patmatch,
                     f(expr),
                     `^`(a::name*b::name,c::nonunit(radnum)),
                     'la'
                   ),
                   printf("Pattern Match with c=%a\n", eval(c, la)),
                   printf("No Match\n")
                 ):
getMatch(expr1);
getMatch(expr2);
getMatch(expr3);

sin(x)*(x*y)^(1/3)

 

sin(x)+(x*y)^(1/6)

 

sin(x)*(x*y)^(1/3)+3

 

proc (ex) options operator, arrow; `if`(type(ex, `+`), `~`[g]([op(ex)]), [op(ex)]) end proc

 

proc (ex) options operator, arrow; `if`(type(ex, `*`), op(ex), ex) end proc

 

Pattern Match with c=1/3
Pattern Match with c=1/6
Pattern Match with c=1/3

 

 


 

Download matchPat2.mw

@Carl Love

In addition to the corrections alreadt noted

  1. There are various places on the right hand sides of the OP's expressions where "extra" spaces are probably being interpreted as "implicit multiplications"
  2. Using the same name for the summation index and the loop index is definitely not advisable
  3. In the attached, I have made several corrections based on "guessing" your intent, and this worksheet executes without error. However there is no guarantee that my "guesses are correct, so I suggest you examine the attached very carefully

restart; _local(gamma); m := 10; A := 10; delta := .112; rho := .23; mu := .31; beta := 1.4; alpha := 2.1; gamma := 1.02; q := 2.3; b1 := 50; b2 := 10; b3 := 5; b4 := 20; S(0) := b1; B(0) := b2; V(0) := b3; R(0) := b4; N := 3; for k from 0 to N do S(k+1) := (A*delta*k-(rho+µ)*S(k)-beta*(sum(S(m)*B(j-m), j = 0 .. m)))/(k+1); B(k+1) := -(-(µ+alpha+gamma)*B(k)+beta*(sum(S(m)*B(j-m), j = 0 .. m)))/(k+1); V(k+1) := (rho*S(k)-(1-q)*S(k)-µ*V(k))/(k+1); R(k+1) := (gamma*B(k)-µ*R(k))/(k+1) end do

gamma

 

10

 

10

 

.112

 

.23

 

.31

 

1.4

 

2.1

 

1.02

 

2.3

 

50

 

10

 

5

 

20

 

50

 

10

 

5

 

20

 

3

 

-11.50-50*µ-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10)

 

10*µ+31.20-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10)

 

76.50-5*µ

 

10.20-20*µ

 

.5600000000-(1/2)*(.23+µ)*(-11.50-50*µ-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-.7000000000*S(10)*B(-10)-.7000000000*S(10)*B(-9)-.7000000000*S(10)*B(-8)-.7000000000*S(10)*B(-7)-.7000000000*S(10)*B(-6)-.7000000000*S(10)*B(-5)-.7000000000*S(10)*B(-4)-.7000000000*S(10)*B(-3)-.7000000000*S(10)*B(-2)-.7000000000*S(10)*B(-1)-7.000000000*S(10)

 

(1/2)*(µ+3.12)*(10*µ+31.20-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-.7000000000*S(10)*B(-10)-.7000000000*S(10)*B(-9)-.7000000000*S(10)*B(-8)-.7000000000*S(10)*B(-7)-.7000000000*S(10)*B(-6)-.7000000000*S(10)*B(-5)-.7000000000*S(10)*B(-4)-.7000000000*S(10)*B(-3)-.7000000000*S(10)*B(-2)-.7000000000*S(10)*B(-1)-7.000000000*S(10)

 

-(1/2)*µ*(76.50-5*µ)-8.797500000-38.25000000*µ-10.71000000*S(10)-1.071000000*S(10)*B(-9)-1.071000000*S(10)*B(-8)-1.071000000*S(10)*B(-7)-1.071000000*S(10)*B(-6)-1.071000000*S(10)*B(-5)-1.071000000*S(10)*B(-4)-1.071000000*S(10)*B(-3)-1.071000000*S(10)*B(-2)-1.071000000*S(10)*B(-1)-1.071000000*S(10)*B(-10)

 

5.100000000*µ+15.91200000-.7140000000*S(10)*B(-10)-.7140000000*S(10)*B(-9)-.7140000000*S(10)*B(-8)-.7140000000*S(10)*B(-7)-.7140000000*S(10)*B(-6)-.7140000000*S(10)*B(-5)-.7140000000*S(10)*B(-4)-.7140000000*S(10)*B(-3)-.7140000000*S(10)*B(-2)-.7140000000*S(10)*B(-1)-7.140000000*S(10)-(1/2)*µ*(10.20-20*µ)

 

.7466666667-(1/3)*(.23+µ)*(.5600000000-(1/2)*(.23+µ)*(-11.50-50*µ-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-.7000000000*S(10)*B(-10)-.7000000000*S(10)*B(-9)-.7000000000*S(10)*B(-8)-.7000000000*S(10)*B(-7)-.7000000000*S(10)*B(-6)-.7000000000*S(10)*B(-5)-.7000000000*S(10)*B(-4)-.7000000000*S(10)*B(-3)-.7000000000*S(10)*B(-2)-.7000000000*S(10)*B(-1)-7.000000000*S(10))-.4666666667*S(10)*B(-10)-.4666666667*S(10)*B(-9)-.4666666667*S(10)*B(-8)-.4666666667*S(10)*B(-7)-.4666666667*S(10)*B(-6)-.4666666667*S(10)*B(-5)-.4666666667*S(10)*B(-4)-.4666666667*S(10)*B(-3)-.4666666667*S(10)*B(-2)-.4666666667*S(10)*B(-1)-4.666666667*S(10)

 

(1/3)*(µ+3.12)*((1/2)*(µ+3.12)*(10*µ+31.20-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-.7000000000*S(10)*B(-10)-.7000000000*S(10)*B(-9)-.7000000000*S(10)*B(-8)-.7000000000*S(10)*B(-7)-.7000000000*S(10)*B(-6)-.7000000000*S(10)*B(-5)-.7000000000*S(10)*B(-4)-.7000000000*S(10)*B(-3)-.7000000000*S(10)*B(-2)-.7000000000*S(10)*B(-1)-7.000000000*S(10))-.4666666667*S(10)*B(-10)-.4666666667*S(10)*B(-9)-.4666666667*S(10)*B(-8)-.4666666667*S(10)*B(-7)-.4666666667*S(10)*B(-6)-.4666666667*S(10)*B(-5)-.4666666667*S(10)*B(-4)-.4666666667*S(10)*B(-3)-.4666666667*S(10)*B(-2)-.4666666667*S(10)*B(-1)-4.666666667*S(10)

 

-.2550000000*(.23+µ)*(-11.50-50*µ-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-(1/3)*µ*(-(1/2)*µ*(76.50-5*µ)-8.797500000-38.25000000*µ-10.71000000*S(10)-1.071000000*S(10)*B(-9)-1.071000000*S(10)*B(-8)-1.071000000*S(10)*B(-7)-1.071000000*S(10)*B(-6)-1.071000000*S(10)*B(-5)-1.071000000*S(10)*B(-4)-1.071000000*S(10)*B(-3)-1.071000000*S(10)*B(-2)-1.071000000*S(10)*B(-1)-1.071000000*S(10)*B(-10))+.2856000000-3.570000000*S(10)-.3570000000*S(10)*B(-9)-.3570000000*S(10)*B(-8)-.3570000000*S(10)*B(-7)-.3570000000*S(10)*B(-6)-.3570000000*S(10)*B(-5)-.3570000000*S(10)*B(-4)-.3570000000*S(10)*B(-3)-.3570000000*S(10)*B(-2)-.3570000000*S(10)*B(-1)-.3570000000*S(10)*B(-10)

 

.1700000000*(µ+3.12)*(10*µ+31.20-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-.2380000000*S(10)*B(-10)-.2380000000*S(10)*B(-9)-.2380000000*S(10)*B(-8)-.2380000000*S(10)*B(-7)-.2380000000*S(10)*B(-6)-.2380000000*S(10)*B(-5)-.2380000000*S(10)*B(-4)-.2380000000*S(10)*B(-3)-.2380000000*S(10)*B(-2)-.2380000000*S(10)*B(-1)-2.380000000*S(10)-(1/3)*µ*(5.100000000*µ+15.91200000-.7140000000*S(10)*B(-10)-.7140000000*S(10)*B(-9)-.7140000000*S(10)*B(-8)-.7140000000*S(10)*B(-7)-.7140000000*S(10)*B(-6)-.7140000000*S(10)*B(-5)-.7140000000*S(10)*B(-4)-.7140000000*S(10)*B(-3)-.7140000000*S(10)*B(-2)-.7140000000*S(10)*B(-1)-7.140000000*S(10)-(1/2)*µ*(10.20-20*µ))

 

.8400000000-(1/4)*(.23+µ)*(.7466666667-(1/3)*(.23+µ)*(.5600000000-(1/2)*(.23+µ)*(-11.50-50*µ-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-.7000000000*S(10)*B(-10)-.7000000000*S(10)*B(-9)-.7000000000*S(10)*B(-8)-.7000000000*S(10)*B(-7)-.7000000000*S(10)*B(-6)-.7000000000*S(10)*B(-5)-.7000000000*S(10)*B(-4)-.7000000000*S(10)*B(-3)-.7000000000*S(10)*B(-2)-.7000000000*S(10)*B(-1)-7.000000000*S(10))-.4666666667*S(10)*B(-10)-.4666666667*S(10)*B(-9)-.4666666667*S(10)*B(-8)-.4666666667*S(10)*B(-7)-.4666666667*S(10)*B(-6)-.4666666667*S(10)*B(-5)-.4666666667*S(10)*B(-4)-.4666666667*S(10)*B(-3)-.4666666667*S(10)*B(-2)-.4666666667*S(10)*B(-1)-4.666666667*S(10))-.3500000000*S(10)*B(-10)-.3500000000*S(10)*B(-9)-.3500000000*S(10)*B(-8)-.3500000000*S(10)*B(-7)-.3500000000*S(10)*B(-6)-.3500000000*S(10)*B(-5)-.3500000000*S(10)*B(-4)-.3500000000*S(10)*B(-3)-.3500000000*S(10)*B(-2)-.3500000000*S(10)*B(-1)-3.500000000*S(10)

 

(1/4)*(µ+3.12)*((1/3)*(µ+3.12)*((1/2)*(µ+3.12)*(10*µ+31.20-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-.7000000000*S(10)*B(-10)-.7000000000*S(10)*B(-9)-.7000000000*S(10)*B(-8)-.7000000000*S(10)*B(-7)-.7000000000*S(10)*B(-6)-.7000000000*S(10)*B(-5)-.7000000000*S(10)*B(-4)-.7000000000*S(10)*B(-3)-.7000000000*S(10)*B(-2)-.7000000000*S(10)*B(-1)-7.000000000*S(10))-.4666666667*S(10)*B(-10)-.4666666667*S(10)*B(-9)-.4666666667*S(10)*B(-8)-.4666666667*S(10)*B(-7)-.4666666667*S(10)*B(-6)-.4666666667*S(10)*B(-5)-.4666666667*S(10)*B(-4)-.4666666667*S(10)*B(-3)-.4666666667*S(10)*B(-2)-.4666666667*S(10)*B(-1)-4.666666667*S(10))-.3500000000*S(10)*B(-10)-.3500000000*S(10)*B(-9)-.3500000000*S(10)*B(-8)-.3500000000*S(10)*B(-7)-.3500000000*S(10)*B(-6)-.3500000000*S(10)*B(-5)-.3500000000*S(10)*B(-4)-.3500000000*S(10)*B(-3)-.3500000000*S(10)*B(-2)-.3500000000*S(10)*B(-1)-3.500000000*S(10)

 

-.1275000000*(.23+µ)*(.5600000000-(1/2)*(.23+µ)*(-11.50-50*µ-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-.7000000000*S(10)*B(-10)-.7000000000*S(10)*B(-9)-.7000000000*S(10)*B(-8)-.7000000000*S(10)*B(-7)-.7000000000*S(10)*B(-6)-.7000000000*S(10)*B(-5)-.7000000000*S(10)*B(-4)-.7000000000*S(10)*B(-3)-.7000000000*S(10)*B(-2)-.7000000000*S(10)*B(-1)-7.000000000*S(10))-1.785000000*S(10)-.1785000000*S(10)*B(-9)-.1785000000*S(10)*B(-8)-.1785000000*S(10)*B(-7)-.1785000000*S(10)*B(-6)-.1785000000*S(10)*B(-5)-.1785000000*S(10)*B(-4)-.1785000000*S(10)*B(-3)-.1785000000*S(10)*B(-2)-.1785000000*S(10)*B(-1)-.1785000000*S(10)*B(-10)-(1/4)*µ*(-.2550000000*(.23+µ)*(-11.50-50*µ-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-(1/3)*µ*(-(1/2)*µ*(76.50-5*µ)-8.797500000-38.25000000*µ-10.71000000*S(10)-1.071000000*S(10)*B(-9)-1.071000000*S(10)*B(-8)-1.071000000*S(10)*B(-7)-1.071000000*S(10)*B(-6)-1.071000000*S(10)*B(-5)-1.071000000*S(10)*B(-4)-1.071000000*S(10)*B(-3)-1.071000000*S(10)*B(-2)-1.071000000*S(10)*B(-1)-1.071000000*S(10)*B(-10))+.2856000000-3.570000000*S(10)-.3570000000*S(10)*B(-9)-.3570000000*S(10)*B(-8)-.3570000000*S(10)*B(-7)-.3570000000*S(10)*B(-6)-.3570000000*S(10)*B(-5)-.3570000000*S(10)*B(-4)-.3570000000*S(10)*B(-3)-.3570000000*S(10)*B(-2)-.3570000000*S(10)*B(-1)-.3570000000*S(10)*B(-10))+.2856000000

 

0.8500000000e-1*(µ+3.12)*((1/2)*(µ+3.12)*(10*µ+31.20-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-.7000000000*S(10)*B(-10)-.7000000000*S(10)*B(-9)-.7000000000*S(10)*B(-8)-.7000000000*S(10)*B(-7)-.7000000000*S(10)*B(-6)-.7000000000*S(10)*B(-5)-.7000000000*S(10)*B(-4)-.7000000000*S(10)*B(-3)-.7000000000*S(10)*B(-2)-.7000000000*S(10)*B(-1)-7.000000000*S(10))-.1190000000*S(10)*B(-10)-.1190000000*S(10)*B(-9)-.1190000000*S(10)*B(-8)-.1190000000*S(10)*B(-7)-.1190000000*S(10)*B(-6)-.1190000000*S(10)*B(-5)-.1190000000*S(10)*B(-4)-.1190000000*S(10)*B(-3)-.1190000000*S(10)*B(-2)-.1190000000*S(10)*B(-1)-1.190000000*S(10)-(1/4)*µ*(.1700000000*(µ+3.12)*(10*µ+31.20-1.4*S(10)*B(-10)-1.4*S(10)*B(-9)-1.4*S(10)*B(-8)-1.4*S(10)*B(-7)-1.4*S(10)*B(-6)-1.4*S(10)*B(-5)-1.4*S(10)*B(-4)-1.4*S(10)*B(-3)-1.4*S(10)*B(-2)-1.4*S(10)*B(-1)-14.0*S(10))-.2380000000*S(10)*B(-10)-.2380000000*S(10)*B(-9)-.2380000000*S(10)*B(-8)-.2380000000*S(10)*B(-7)-.2380000000*S(10)*B(-6)-.2380000000*S(10)*B(-5)-.2380000000*S(10)*B(-4)-.2380000000*S(10)*B(-3)-.2380000000*S(10)*B(-2)-.2380000000*S(10)*B(-1)-2.380000000*S(10)-(1/3)*µ*(5.100000000*µ+15.91200000-.7140000000*S(10)*B(-10)-.7140000000*S(10)*B(-9)-.7140000000*S(10)*B(-8)-.7140000000*S(10)*B(-7)-.7140000000*S(10)*B(-6)-.7140000000*S(10)*B(-5)-.7140000000*S(10)*B(-4)-.7140000000*S(10)*B(-3)-.7140000000*S(10)*B(-2)-.7140000000*S(10)*B(-1)-7.140000000*S(10)-(1/2)*µ*(10.20-20*µ)))

(1)

NULL


 

Download badSums2.mw

 

in the Mapleprimes toolbar to upload a worksheet with your calculation

That way the problem can be sensibly invetigated

@Carl Love 

for the wrong question! It was intended for

https://www.mapleprimes.com/questions/227415-How-To-Convert-The-Tekplot-File?reply=answer

because I only have the last six major versions of Maple installed on this machine.

Code I supplied works fine on all of them, but that only goes back to Maple 18!

 

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