HDTenEightyP/NetNews huggingface dataset:
From:  pat@office.relcom.msk.su (Pataraia A.D.)
Subject: Substitutions in Mma 
Reply-To: pat@office.relcom.msk.su
Organization: RELCOM AG
Date: Fri, 16 Oct 92 06:15:01 +0300 

In Reduce, function SUB(x=y,f) makes
substitution AFTER evaluating f.
Is there a similar function available  in Mma?
( Block[{x=y},f] is not just the same as
Sub(x=y,f) . It makes substitution BEFOR the evaluation of
the expression. For example, Block[{x=5},D[x,x]] gives D[5,5],
while SUB(x=0,DF(x,x)) gives 1 ).
-- 
_______________________________________________________________________
                       Alexander Pataraia 
   ( Institute of Nuclear Research of Russian Academy of Sciences )
E-mail addresses: pat@office.relcom.msk.su and pat@rita.reladmin.msk.su
From: rainbow@ubu.Eng.Sun.COM (Linda Kaplan)
Newsgroups: sci.math.symbolic
Subject: Mathematica
Date: 16 Oct 92 21:43:15 GMT
Organization: Sun Microsystems, Mt. View, Ca.
NNTP-Posting-Host: ubu

Excuse me please.  Is there a newsgroup appropriate for
discussions about Mathematica, the program?  Or is this
(sci.math.symbolic) the closest newsgroup for such discussion?

Thanks,

Rainbow

--
"Anyone who considers arithmatical methods of producing
random digits is, of course, in a state of sin"
			  -Norman
Newsgroups: sci.math.symbolic
From: us291442@mmm.serc.3m.com (Mark B. Mizen)
Subject: Mathematica 2.1 for Windows
Organization: 3M Company
Date: Fri, 16 Oct 92 18:24:18 GMT


I just had a very unpleasant experience with Mathematica and Wolfram
Research and thought I should report my experience to the net.  I am 
using Windows 3.1, DOS 5, and a postscript printer.

Mathematica refuses to completely print cells if either the input or output 
cells exceed one page in length.  The program simply prints multiple copies
of the first output page and will not print subsequent pages.  I called
Mathematica technical support and was told to import the cells and print
them from the word processor!  I told technical support that I did not 
consider this a satisfactory solution for a $1000 program.  In addition
brackets which group cells are not included with the imported cell.  I 
then asked if they could send me a corrected version when available.  I
was then told "We can't promise that."  At that point I hung up in 
disgust.  Wolfram Research seems to have lost sight of their customers.
Am I the only user who expects to be treated better than this?  I will
post any future response from Wolfram.

Mark Mizen
3M Center Bldg 201-3N-04
St. Paul, MN 55101
(612) 733-3925
mbmizen@mmm.com







Newsgroups: sci.math.symbolic
From: taylor@sun.soe.clarkson.edu (Ross Taylor,Peyt312,6652,2655436)
Subject: Thanks for help on Simple (?) Maple question and Thermodynamics
Organization: Clarkson University
Date: Thu, 15 Oct 1992 17:01:23 GMT

My thanks to those of you that helped my question on Maple.

I found that the combine function (as suggested by several) did not 
do exactly what I wanted but the following script that uses combine
did do the trick.

P:=n*R*T/(V-B)-D/(V-alfa*B)/(V-beta*B);
Ar(T,V,n) := int( (P - (n*R*T)/V),V=V0..infinity);
Ar(T,V,n) := subs(V0=V,Ar(T,V,n));
Fr(T,V,n) := Ar(T,V,n) / (R*T);
Fr(T,V,n):=simplify(Fr(T,V,n),power);
Fr(T,V,n):=collect(Fr(T,V,n),[ln(V),ln(V-B)],distributed,factor);
t1:=combine(ln(V-B)-ln(V),ln);
Fr(T,V,n):=subs(ln(V-B)=t1+ln(V),Fr(T,V,n));
t2:=combine(ln(-V+alfa*B)-ln(-V+beta*B),ln);
Fr(T,V,n):=subs(ln(-V+alfa*B)=t2+ln(-V+beta*B),Fr(T,V,n));
Fr(T,V,n):=collect(Fr(T,V,n),ln((V-B)/V),distributed,factor);
Fr(T,V,n):=subs((-V+alfa*B)/(-V+beta*B)=(V-alfa*B)/(V-beta*B),Fr(T,V,n));


Ross Taylor
From: jonesm@nic.cerf.net (Matthew Jones)
Newsgroups: sci.math.symbolic
Subject: How do you display a histogram in Mathematica
Keywords: Mathematica, Histogram
Date: 8 Oct 92 14:24:22 GMT
Distribution: usa
Organization: CERFnet

Subject about sayes it.
How do you graph a histogram in Mathematica? I can't even find a 
way to partition the set of points into bins and then do a ListPlot.
Any help would be appreciated.

Matthew Jones
jonesm@cerf.net
Newsgroups: sci.math.symbolic
From: ma1346631@ntuvax.ntu.ac.sg
Subject: Help required to use Maple `solve'
Organization: Nanyang Technological University - Singapore
Date: Thu, 8 Oct 1992 05:11:32 GMT

I have some problem using the Maple command `solve'. 

For example, if I write
>solve({x^2-4},{x});

it will give me the following
	{x:=2,x:=-2}

My problem is: how do I put these 2 values into separate variable 
names ie. x1:=2 and x2:=-2; so that I can use these 2 values separately
and for different purposes?
From: ee52fas@sdcc15.ucsd.edu (Anand Rao)
Newsgroups: sci.math.symbolic
Subject: mma
Date: 15 Oct 92 23:18:34 GMT

does anybody know if i can take circuit drawings made with
Logicworks (or Spice) and import them into a mathematica document?
(on a macintosh) thanks....
Newsgroups: sci.math.symbolic
From: mcgrant@rascals.stanford.edu (Michael C. Grant)
Subject: Re: 2D convolution, 2D fourier w/ Mma
Organization: Information Systems Laboratory, Stanford University
Date: 13 Oct 92 14:43:41

In article <1992Oct13.175039.6624@bmw.mayo.edu>
lepard@mayo.edu (Kevin O. Lepard) writes:

   I am starting to use Mathematica in a 2D DSP class I'm taking, and I am trying
   to do a 2D convolution.  The Mma book discusses doing this by multiplying the
   fourier transforms of the 2 data sets then taking the inverse fourier
   transform of that.  I understand that.  However, I can't get an answer that
   matches any of the examples in the book (the answers aren't even close to
   each other).

   As an example, let's say I want to convolve

   7 8 9
   4 5 6	(where 1 is located at coordinates (n1,n2) = (0,0), 2 is at (1,0),
   1 2 3	 5 is at (1,1), etc.)

   with

   -1 -1
    1  1 	(where the lower left hand corner is (n1=0, n2=0))

   The answer in the book (Two-Dimensional Signal and Image Processing, by J.S. Lim)
   is:

   -7 -15 -17 -9
    3   6   6  3
    3   6   6  3
    1   3   5  3

   If anyone can help me understand how to do this, I'd be grateful.  I've read
   what it has to say in _Mathematica_A_System_for_doing_Mathematics_by_Computer_
   (several times) but am apparently too dense to understand it.

   I don't have a notebook front end, so if someone could include the syntax
   required to do this from the command line, that would really help a lot.

   Thanks in advance.

Well, I don't know if Mathematica has these functions but this is what I
did in MATLAB:

A = [7,8,9;4,5,6;1,2,3]
B = [-1,-1;1,1]
C = real(ifft2(fft2(A,4,4).*fft2(B,4,4)))

and got the right answer. The form fft2(A,4,4) must be used to pad the data
with zeros out to a 4 by 4 dimension, since that is the final size of the
linear convolution. The ".*" op is element-by-element multiplication, which
I believe is equivalent to "*" in Mathematica for lists/vectors.

If I were to do the padding by hand I'd get

A = [7,8,9,0;4,5,6,0;1,2,3,0;0,0,0,0]
B = [-1,-1,0,0;1,1,0,0;0,0,0,0;0,0,0,0]
C = real(ifft2(fft2(A).*fft2(B)))

Mike
From: lamouche@ERE.UMontreal.CA (Lamouche Guy)
Newsgroups: sci.math.symbolic
Subject: Mathematica - Editor
Date: 15 Oct 92 23:08:23 GMT
Organization: Universite de Montreal


HI.
I am using mathematica 2.0 on a Silicon Graphics machine (Irix4-D).

when I am using the EditIn[] function to reedit previously entered 
expressions, mathematica sends me into vi on my machine. Is there
an environment variable I could set so mathematica would be gentle enough
to call emacs? I haven't found anything in the Mathematica book about
this.
Xref: sparky sci.math:12952 ucb.math:2 sci.math.stat:2097
From: judya@sag4.ssl.berkeley.edu (Judy Acereto)
Newsgroups: sci.math,ucb.math,sci.math.stat
Subject: Visual Math Talk 10/12
Date: 9 Oct 1992 04:46:10 GMT
Organization: Center for EUV Astrophysics
Distribution: ba
NNTP-Posting-Host: sag4.ssl.berkeley.edu

LEONARDO/YLEM LECTURE: MICHELE EMMER ON VISUAL MATHEMATICS

Professor Michele Emmer from Italy will give a talk on Visual 
Mathematics and other aspects of art and mathematics on Monday, 
October 12 at 7:30pm.  Prof. Emmer is the Guest Editor of the recent 
Special Issue on Visual Mathematics of _Leonardo_, the magazine for 
artists using science and technology.  For further information contact 
Dr. Roger Malina : rmalina@cea.berkeley.edu 

The talk will be given at the Center for EUV Astrophysics, University 
of California, Berkeley.  The street address is 2150 Kittredge St., 
Berkeley.
Newsgroups: sci.math
Subject: Re: Array Algebra
From: jeremy@cs.aukuni.ac.nz (Jeremy Gibbons)
Date: Tue, 13 Oct 1992 05:45:29 GMT
Organization: Computer Science Dept. University of Auckland
Cc: jeremy@cs.aukuni.ac.nz

In <1992Oct12.202012.7559@iscsvax.uni.edu> kraai4712@iscsvax.uni.edu writes:

>I would be happy to entertain any discussion about mathematics of arrays with
>an arbitrary number of (integer) dimensions.  Can/has it been done?  Who,  in
>what publication, when?  The obvious would be:
>        scalars   zero dimensions
>        vectors   one dimension
>        matrices  two dimensions

Colin Banger <colin@qucis.queensu.ca> models integer-dimensioned arrays as
an infinite list of elements and a finite list of dimensions; paper in
ATABLE-92 (Montreal; don't know a more precise ref). 

Trenchard More's Array Theory discusses nested heterogeneous arrays;
several papers in APL Quote Quad, implemented in Nial and Q'Nial.

Lenore Mullin's <lenore@uvm-gen.emba.uvm.edu> Mathematics of Arrays (PhD,
Syracuse, 1988) is another APL-based theory, about non-nested arrays.

Richard Bird (PRG, Oxford) discusses two-dimensional arrays constructed
from two "concatenation" operators, "above" and "beside" (Lectures on
Constructive Functional Programming, in M Broy, ed, Constructive Methods in
Computer Science, Springer, 1989).  This is work with Chris Wright, who
hasn't published anything (yet).

>I am very interested in this as I have found some convenient ways to represent
>some combinatorial problems very nicely in this light.  Some relationships
>acquire a very intuitive feel when the spatial relationships become tangible
>through the use of such models.

I'd be interested in hearing your ideas. Have you written anything up yet?

Jeremy

---
Jeremy Gibbons <jeremy@cs.aukuni.ac.nz>         tel: +64 9 373 7599
   Department of Computer Science,              fax: +64 9 373 7453
   University of Auckland, Private Bag 92019, Auckland, New Zealand.
From: stephen@mont.cs.missouri.edu (Stephen Montgomery-Smith)
Newsgroups: sci.math
Subject: Re: Probably you know the probability ;-)
Date: 15 Oct 92 16:52:50 GMT
Organization: University of Missouri

In <1992Oct15.045308.3552@noose.ecn.purdue.edu> kavuri@lips.ecn.purdue.edu (Surya N Kavuri ) writes:

>In article <1992Oct15.025957.11871@galois.mit.edu>, jbaez@riesz.mit.edu (John C. Baez) writes:
>> In article <1992Oct14.193802.25129@noose.ecn.purdue.edu> kavuri@lips1.ecn.purdue.edu (Surya N Kavuri ) writes:
>> >  I am given a "rod" of fixed length l. I took a hammer and
>> >  broke it into three parts.
>> >
>> >   (a) what is the probability that they form a triangle ?
>> >   (b) what is the probability that they form an isoceles triangle ?
>> 
>> There is not enough information to answer unless you tell us the
>> probability distribution of the lengths of the 3 parts.  In real
>> life this would depend in a complicated way on how you broke the
>> rod.
>> 

>  All the points on the rod are equally likely to be breakage points.
>  

It seems to me that one would break the rod into lengths (a,b,c) such
that all admissible triplets (i.e. a+b+c = l) are equally likely.  I
think that if you pick the break points uniformly along the rod that
you get a different distribution.

Well, I haven't done any calculations, so I might be wrong.  But thats the
beauty of nn --- write first, think later !!!!

Am I right?

Stephen
Newsgroups: sci.math.stat
From: thompson@atlas.socsci.umn.edu (T. Scott Thompson)
Subject: Re: Help with a problem...
Reply-To: thompson@atlas.socsci.umn.edu
Organization: Economics Department, University of Minnesota
Date: Mon, 12 Oct 1992 15:49:36 GMT

ytakahashi@eagle.wesleyan.edu writes:

>	Three people work to solve a problem independantly with differing
>chances of solving it.  Given these probabilities, what is the probability
>that the problem is solved?

This looks a lot like homework so I won't give a complete answer:

>Since they are independent, P(AB) = P(A)*P(B)  by definition
>I extended this to... P(ABC) = P(A)*P(B)*P(C)
>	(I'm not sure if that's legal to do)

Yes, this is legal.  It is the definition of independence for 3 events.

Aside: P(AB)=P(A)*P(B), P(BC)=P(B)*P(C) and P(AC)=P(A)*P(C) does _not_
imply that P(ABC)=P(A)*P(B)*P(C).  So although you are using the
correct expression for P(ABC), the reasoning that got you there is not
necessarily correct.

>P(AUBUC) = P(A)+P(B)+P(C) - P(ABC)   by definition

>plug in from above yields...
>P(AUBUC) = P(A)+P(B)+P(C) - P(A)*P(B)*P(C)

But this is wrong.  (If you are not convinced, try P(A)=P(B)=P(C)=1/2
in the formula.  Then the right-hand side is 3/2 - 1/8 = 11/8, which
is not a probability.)  You are incorrectly assuming here that

P(AUBUC) = P(A) + P(B) + P(C) - P(ABC)

I will leave it up to you to figure out why this is not correct, and
what correction is needed.

Hint: Write out the 2*2*2 = 8 different outcomes and count how many
times each one is included on the right-hand side of your formula.  A
Venn diagram of the possibilities will help.
--
T. Scott Thompson              email:  thompson@atlas.socsci.umn.edu
Department of Economics        phone:  (612) 625-0119
University of Minnesota        fax:    (612) 624-0209
From: hougen@vision.csl.uiuc.edu (Darrell Roy Hougen)
Newsgroups: sci.math.stat
Subject: Re: Goodness of Fit with Zero Degrees of Freedom
Date: 8 Oct 92 19:11:34 GMT
Article-I.D.: roundup.1b2156INNnnf
Organization: Center for Reliable and High-Performance Computing, University of Illinois at Urbana-Champaign
NNTP-Posting-Host: vision.csl.uiuc.edu

chaas@coe.drexel.edu (Chuck Haas) writes:


*>Suppose I have a set of categorical data.
...
*> If the model has k parameters and there are m groups, then I know
*>that I could do a goodness of fit test (e.g., Fisher chi-square or
*>likelihood ratio) and compare to the chi-square distribution with m-k
*>degrees of freedom.

I'm not sure I understand your problem, but shouldn't the number of
degrees of freedom be equal to the number of data rather than the
number of categories?

Darrell
Xref: sparky sci.math.stat:2104 sci.math.num-analysis:2983
From: ravi@bobcat.ent.ohiou.edu (Pattalachinti)
Newsgroups: sci.math.stat,sci.math.num-analysis
Subject: Re: Help! !!!!!!!!!!!!!
Date: 11 Oct 92 06:20:34 GMT
Organization: College of Engg. & Tech., Ohio University, Athens, Ohio

Can anybody out there have some very  rubost and 
efficient methods which can be executed on PC.  The language used 
is Fortran 77.

It is very much appreciated if you folk suggest me some books which
deal with non-linear equation solving algorithms which work 
efficiently.

my address::ravi@bobcat.ent.ohiou.edu

Thanking in advance,
ravi

My belief
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
GRAVITY EXISTED BEFORE NEWTON DISCOVERED IT       !!!!!!!!!!!!!  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Newsgroups: sci.math.stat
From: hrubin@pop.stat.purdue.edu (Herman Rubin)
Subject: Re: Rational functions: a linear model?
Organization: Purdue University Statistics Department
Date: Mon, 12 Oct 1992 14:30:50 GMT

In article <1992Oct11.192921.22075@zip.eecs.umich.edu> ayman@eecs.umich.edu (Ayman I. Kayssi) writes:
>hrubin@pop.stat.purdue.edu (Herman Rubin) wrote:
>>In article <1992Oct11.153538.12821@zip.eecs.umich.edu> ayman@eecs.umich.edu (Ayman I. Kayssi) writes:
>>>y = a0 + a1*x + ... + am*x^m + b1*(-y)*x + b2*(-y)*x^2 + ... + bn*(-y)*x^n

>>If you are trying a regression on your transformed equation, the estimates
>>will not even be consistent because of dependent variables on the right side.


>TableCurve (from Jandel) lists the following form under *linear* equations:

>    a + c x + e x^2 + ...
>y = ---------------------
>    1 + b x + d x^2 + ...

>How do they (non-iteratively) calculate the parameters?

I am not familiar with the work; I do not even know if data is being fitted,
or a theoretical function.  A standard numerical analysis situation is the
fitting of approximations of fixed type; as infinite series, the problem 
is hopelessly unidentified.  As I have stated before, there are quite 
substantial numerical problems, although these may not be too crucial for
the degree of fit.  If the fit is to observational data, it may even be a
good idea to use ridge regression or some other Bayesian, or approximate
Bayesian, approacth.

But I cannot see any good reason why it is important, or even desirable, to
calculate non-iteratively, unless this is to be a "field" situation, with
thousands or millions of cases to be computed.
-- 
Herman Rubin, Dept. of Statistics, Purdue Univ., West Lafayette IN47907-1399
Phone: (317)494-6054
hrubin@pop.stat.purdue.edu (Internet, bitnet)  
{purdue,pur-ee}!pop.stat!hrubin(UUCP)
From: freeman@olympus.as.arizona.edu (Jonathan Freeman)
Newsgroups: sci.math.stat
Subject: Need Reference for Iteratively Reweighted LS
Date: 15 Oct 92 18:32:39 GMT
Distribution: usa
Organization: University of Arizona, Tucson, AZ

I have been using iteratively reweighted least squares to fit a non-linear
model with a number of parameters. The errors are unequal so I use
weighted least squares. The problem I have is that I am unsure of the 
final scaling factor. Usually this scaling factor involves a square and 
derivative of the influence function used. I have a number of papers (Huber, Street et al) which discuss this, but they are all a bit fuzzy. Can 
anyone give me a VERY straightforward reference preferably in 
a journal. Again I want a good reference for weighted multi-parameter
IRLS and am unsure of how to obtain the final errors on the parameter
estimates. Thanks! e-mail or post.

                 Jonathan Freeman
                 Steward Observatory 
	         University of Arizona
Newsgroups: sci.math.stat
From: rob@med.umich.edu (Robert Strawderman)
Subject: Re: A question or two re: splines
Date: Wed, 07 Oct 92 11:59:20 EDT
Organization: University of Michigan Health Sciences, Ann Arbor
Distribution: umich
Reply-To: rob@med.umich.edu (Robert Strawderman)
Originator: rob@gaspra


The way you suggest is the one that I am currently using - I had a feeling that the
procedure I posed would not be well-behaved - you've confirmed my suspicions.

In regard to the "knot-selection" procedure you posed, there are two such references
which come to mind:  Pat Smith's paper (a NASA tech report) on using B-splines in 
regression, and a paper by Kooperberg and Stone which appeared recently in one
of the stat journals (Comp Stat and Data Analysis, I believe).

Thanks for your reply...

-------------------------------------------------------------------------------
*******************************************************************************
\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\

Rob Strawderman				Email:  rob@biostat.med.umich.edu
Department of Biostatistics  		Office:	(313) 936 - 1002
University of Michigan

              "Let no man lack for the thirst of real ale...."

///////////////////////////////////////////////////////////////////////////////
*******************************************************************************
-------------------------------------------------------------------------------
Newsgroups: sci.math.stat
From: ecaxron@ariel.lerc.nasa.gov (Ron Graham)
Subject: TQM Mailing List being formed
News-Software: VAX/VMS VNEWS 1.41    
Organization: NASA Lewis Research Center
Date: 13 Oct 1992 12:40 EST  

Based on comment I have seen in sci.engr that (a) a TQM Mailing List
is desirable and (b) none currently exists, I am going to work on forming one.

Topics for such a Mailing List could include, but not be limited to:

o  definitions of TQM, TQL, or whatever you call it;
o  methods for building quality into processes;
o  the works of Deming, Taguchi, and other quality wizards;
o  balancing technical precision with customer relations.

And so on.  If you are interested, send me your e-mail address.  Right now
I have no idea how exactly to pull this off, but I see a digest format as
most likely at this point, because I do not have list server capability 
that I know of.  If you would like to donate list server capability, let me
know.  Otherwise, it'll probably be digest format unless the distribution
gets too large.

RG

"Drive Out Fear."  -  from Deming's 14 Points
From: Barry.John.Good@brunel.ac.uk (Barry John Good)
Newsgroups: sci.math.stat
Subject: Partitioned clustering algorithm
Date: 15 Oct 92 09:00:42 GMT
Organization: Brunel University, West London, UK


I am posting this on behalf of a client.
Can anybody tell me of a statistical package that offers
partitioning methods for cluster analysis.
All the packages I have access to only offer
the hierarchical approach.
 
 Please reply by e-mail; I seldom look at this newsgroup.
  
  Barry Good
Newsgroups: sci.math.stat
From: ssmithap@cc.curtin.edu.au
Subject: Root Solver Program
Organization: Curtin University of Technology
Date: Mon, 12 Oct 1992 05:50:15 GMT


  Hi!

  Could anyone tell me a name of a program for an MSDOS based machine
  that can find the roots of an equation up to the nth order. (ideally
  up to 10th order or better).  Where abouts would I be able to ftp 
  such a program from.

  Could you please e-mail me your replies.

  Your help would be greatly appreciated

  Thanks

   Andrew
   ssmithap@cc.curtin.edu.au
From: jkiley@andy.bgsu.edu (James H. Kiley)
Newsgroups: sci.math.stat
Subject: Trinomials. :)
Keywords: Help, of course. :)
Date: 16 Oct 92 19:02:00 GMT
Organization: Bowling Green State University B.G., Oh.


	For a game I play (yes, it's only a game, I know, why am I troubling 
all of you with it?  Well, somebody out there might know...)  it becomes
helpful if not necessary to know a good expansion formula for (a+b+c)^N. 
Now I'm aware that the expansion 'formula' for (a+b)^N is just Pascal's Triangle,
but I'm wondering if any of you wizards out there know of something similar
for (a+b+c)^N without my having to go through the whole sadistic algebraic
death of actually deriving each and every one for N=1 to 10?
	Thanks in advance, that is if anyone actually knows of one...

	Jim Kiley
	jkiley@andy.bgsu.edu
From: mglacy@lamar.ColoState.EDU
Newsgroups: sci.math.stat
Subject: Another RNG question
Date: 16 Oct 92 19:59:57 GMT
Organization: Colorado State University, Fort Collins, CO  80523

I've been trying out a new RNG code. It appears to yield  a very
nice random u.r.v (0,1).  (At least it looks fine under k-s tests,
which is sufficient for my purposes.)  However,  when I use it  to
generate uniformly distributed random *integers*, it produces lousy
results.  I've just been multiplying each random real by 
some N, and then truncating the result to an integer. To make
the range be 1..N, I've then been adding 1 to the truncated result.

Any suggestions on why might be going wrong here?  Is this to
be expected?

Thanks.
=-=-=-=-=-=-=-=-=-==-=-=-=
Mike Lacy, Sociology Dept.      CSU and I have a deal: We don't take
Colo. State Univ.               responsibility for one another's opinions.
Fort Collins CO 80523           voice (303) 491-6721
From: dahlhaus@nariwork15.ethz.ch (Dirk Dahlhaus,ETF D111,2791,,,)
Newsgroups: sci.math.stat
Subject: search for spectral estimators
Keywords: spectral estimator
Date: 14 Oct 92 11:06:20 GMT
Reply-To: dahlhaus@nariwork15.ethz.ch
Organization: Sun Microsystems, Inc.

I search for a tool in mathematica, maple, matlab etc.
to do nonparametric spectral estimation from a sequence
of N data points. Is there anyone to help me?
Thanks Dirk

=========================================================
     Swiss Federal Institute of Technology, Zurich
             Communication Technology Lab
     Tel.: +41 1 256-2791      Fax: +41 1 262-0943
             email: dahlhaus@nari.ikt.ethz.ch
=========================================================
Xref: sparky sci.math.stat:2119 sci.math.num-analysis:2994