    function [y]=f1(x,c);
%
% Program to compute function y=x^2-c
%
% INPUT
% x      = vector x points
% c      = constant c
% OUTPUT
% y 
%

    y=x.^2-c;
