On Wed, Mar 14, 2012 at 8:34 AM, Laurent wrote:
> #! /usr/bin/sage -python
> # -*- coding: utf8 -*-
>
> from sage.all import *
> from scipy import stats
> import numpy
>
> X=stats.norm.rvs(size=300)
> K=((heaviside(x+1)-heaviside(x-1))).function(x)
>
> f = sum( [ K( (Xi-x) ) for Xi in X ] )
> p
Il 14/03/2012 08:09, Vegard Lima ha scritto:
Hi,
On Wed, Mar 14, 2012 at 7:56 AM, Laurent wrote:
I'm searching for a function
f(x)=1 if x in [a,b], =0 otherwise
how about using the Heaviside step function?
oh yes, of course it works. Thanks.
...
well
My residual problem is that thi