function f = fja(m,l,gama,b);
  f = tan(m*l)+m/(gama+b);
  return
  
