%file: f1.m
 function [out] = f1(x)
   out = (x-0.98543)*(x-pi)*(x-exp(1))*(x+7.0);
 return;

