function [f]=f(t,w)

%f=w-t^2+1;

f=-w+t+1;