function [y]=f(x)
y=1+x^2
endfunction