Jeroen Nijhof wrote:
> Hello,
>
> Ci(t) fails when t = [0:0.1:10] : it complains about an array size
> mismatch because of this line in 
> /trunk/octave-forge/main/specfun/inst/Ci.m (r4585):
>
>    y(real(z)>0) = -0.5.*(expint_E1(i.*y(real(z)>=0)
> )+expint_E1(-i.*y(real(z)>=0) ));
>
> the z=0 element _is_ included on the left hand side, but not on the
> left. That should probably be replaced by
>
>    y(real(z)>=0) = -0.5.*(expint_E1(i.*y(real(z)>=0)
> )+expint_E1(-i.*y(real(z)>=0) ));
>
> Best regards,
>
> Jeroen Nijhof
>   

Ok, I committed this fix to the SVN

D.


-- 
David Bateman                                [EMAIL PROTECTED]
Motorola Labs - Paris                        +33 1 69 35 48 04 (Ph) 
Parc Les Algorithmes, Commune de St Aubin    +33 6 72 01 06 33 (Mob) 
91193 Gif-Sur-Yvette FRANCE                  +33 1 69 35 77 01 (Fax) 

The information contained in this communication has been classified as: 

[x] General Business Information 
[ ] Motorola Internal Use Only 
[ ] Motorola Confidential Proprietary


-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Octave-dev mailing list
Octave-dev@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/octave-dev

Reply via email to