Date: Oct 15, 2013 2:14 PM
Author: Nasser Abbasi
Subject: Re: complex results using ODE45
On 10/15/2013 12:53 PM, Jeffrey wrote:

> I have the two m files given below. "a" and "T" are returned as complex with

>

> Why are values becoming complex?

>

> Please help.

>

> Jeff

> _______________________________

> %this is in file alpha.m

> function dadT = alpha(T,a)

> b = 20

> l= 50.1

> k=exp(l)

> ner = -36108.73

> n=1.269

> dadT = (k/b)*(exp(ner/T)) *((1-a)^n

> ________________________________________

>

You have missing ")" above btw. Last line.

But besides this, which I assume a cut/paste issue only,

look at the expression itself.

What happens for (1-a)^n, and "a" happens to be

larger than one?

EDU>> (1-1.5)^(1.269)

ans =

-0.2754 - 0.3104i