[petrellove]

% Test image: the usual suspect...
im = imread('lena512.bmp');
im = double(im)/256 ;

% Generate noisy image.
rho=3;
sig = std(im();
sigma = sig / rho;
nim = im + sigma * randn(size(im));