clear all
close all

f = @(x) x.^2 - 9;
opt = optimset('Display', 'iter');

fzero(f, 0, opt)