X-Received: by 10.176.84.3 with SMTP id n3mr2216929uaa.28.1478166064352; Thu, 03 Nov 2016 02:41:04 -0700 (PDT) X-Received: by 10.157.8.134 with SMTP id 6mr818276otf.17.1478166064262; Thu, 03 Nov 2016 02:41:04 -0700 (PDT) Path: csiph.com!v102.xanadu-bbs.net!xanadu-bbs.net!news.glorb.com!p16no455434qta.1!news-out.google.com!z32ni979qtc.1!nntp.google.com!n6no455961qtd.0!postnews.google.com!glegroupsg2000goo.googlegroups.com!not-for-mail Newsgroups: comp.graphics.apps.gnuplot Date: Thu, 3 Nov 2016 02:41:03 -0700 (PDT) In-Reply-To: Complaints-To: groups-abuse@google.com Injection-Info: glegroupsg2000goo.googlegroups.com; posting-host=130.56.80.185; posting-account=eO0UFAoAAAAttewk0-_prR-cK1_CGaJg NNTP-Posting-Host: 130.56.80.185 References: User-Agent: G2/1.0 MIME-Version: 1.0 Message-ID: <65025da6-a456-46d5-b987-2b6d8dcbc692@googlegroups.com> Subject: Re: fitting for Bateman equation returns NaN From: apsara@connect.hku.hk Injection-Date: Thu, 03 Nov 2016 09:41:04 +0000 Content-Type: text/plain; charset=UTF-8 Xref: csiph.com comp.graphics.apps.gnuplot:3430 On Thursday, November 3, 2016 at 2:12:59 AM UTC+11, Michael Okuntsov wrote: > > N=7998.68795 > > M=8540.4184 > > A=.04951051 > > > > direct(x)=M*exp(-B*x) > > fed(x)=(A*N/(B-A))*(exp(-A*x)-exp(-B*x)) > > f(x)=fed(x)+direct(x) > > Define initial parameter: > > B=-0.1 > ... Hi Michael, Thank you! it works great now. Why do we define B, though, and is the -0.1 arbitrary? What do I define N or M as, if I want to fix A, B and fit via N or M? Thanks :)