Carl Love

Carl Love

28075 Reputation

25 Badges

13 years, 58 days
Himself
Wayland, Massachusetts, United States
My name was formerly Carl Devore.

MaplePrimes Activity


These are replies submitted by Carl Love

@digerdiga We can compare the residual sums of squares, and the nonlinear fit is better, but only slightly better.

 

restart:

X:= [51, 61, 71, 81, 101, 121, 141, 161, 201]:

Y:= [1.041, 1.034, 1.0296, 1.0259, 1.021, 1.017, 1.015, 1.013, 1.010]:

Sol1:= Statistics:-Fit(
     1+a*n^b, X, Y, n,
     output= [leastsquaresfunction, residualsumofsquares]
);

[1+HFloat(2.1045354792676307)/n^HFloat(1.0016417508071458), 0.462254043633566e-6]

f1, rsos1:= Sol1[]:

X1:= ln~(X):  Y1:= ln~(Y -~ 1):
Sol2:= Statistics:-Fit(a1+b*n, X1, Y1, n, output= [parametervalues, standarderrors]);

Sol2 := [[a1 = .799926272563527, b = -1.01422767794265], Vector[row](2, {(1) = 0.567833227040259e-1, (2) = 0.122755809934050e-1})]

f2:= eval(1+exp(a1)*n^b, Sol2[1]);

1+HFloat(2.22537685111362)/n^HFloat(1.0142276779426467)

rsos2:= add((eval(f2, n= X[k])-Y[k])^2, k= 1..nops(X));

HFloat(5.959559940306572e-7)

plots:-display([
     plot([f1, f2], n= 50..210),
     plots:-pointplot(zip(`[]`, X, Y))
]);

 

Download Linearized_Fit.mw

I updated the procedure substantially to include the eigenvalues, and I added some color to improve the eye-tracking of the circles. I also added a second example, where the disks are more separated.

I updated the procedure substantially to include the eigenvalues, and I added some color to improve the eye-tracking of the circles. I also added a second example, where the disks are more separated.

@DJJerome1976 

I agree with DJ Jerome's skepticism about RealDomain. I am particularly skeptical about RealDomain:-solve. That's why I said to put

with(RealDomain, `^`);

rather than

with(RealDomain);

into your initialization file. The former will overload only the `^` operator, which is enough to accomplish everything discussed in this thread and which doesn't introduce any problems that I know of.

@DJJerome1976 

I agree with DJ Jerome's skepticism about RealDomain. I am particularly skeptical about RealDomain:-solve. That's why I said to put

with(RealDomain, `^`);

rather than

with(RealDomain);

into your initialization file. The former will overload only the `^` operator, which is enough to accomplish everything discussed in this thread and which doesn't introduce any problems that I know of.

@Alina8faq Yes, there is a way that you can use I as a variable in Maple 13. First, choose a symbol to use as the imaginary unit. Let's say that you choose J (although it need not be a single letter). Then give the command

interface(imaginaryunit= J);

Now you are free to use I however you like.

Even though you use the I with subscripts, it is still confusing to algsubs. But all will be okay if you give the command above.

@Alina8faq Yes, there is a way that you can use I as a variable in Maple 13. First, choose a symbol to use as the imaginary unit. Let's say that you choose J (although it need not be a single letter). Then give the command

interface(imaginaryunit= J);

Now you are free to use I however you like.

Even though you use the I with subscripts, it is still confusing to algsubs. But all will be okay if you give the command above.

@ANANDMUNAGALA Please note that I have significantly updated the procedure since you posted the above Reply.

@ANANDMUNAGALA Please note that I have significantly updated the procedure since you posted the above Reply.

@oxy Even using a variety of assumptions, I cannot simplify the expression or convert it to StandardFunctions.

@oxy Even using a variety of assumptions, I cannot simplify the expression or convert it to StandardFunctions.

@spradlig RealDomain does make your plot work the way that you want, and it does make 4^(1/2) evaluate to 2. Did you try it? If you want these behaviours to be the default, then include this line in your initialization file:

with(RealDomain, `^`):

@spradlig RealDomain does make your plot work the way that you want, and it does make 4^(1/2) evaluate to 2. Did you try it? If you want these behaviours to be the default, then include this line in your initialization file:

with(RealDomain, `^`):

It will work for the plain numeric expressions, but not the plots.

It will work for the plain numeric expressions, but not the plots.

First 600 601 602 603 604 605 606 Last Page 602 of 709