/* ================================================================  */ 
#include <stdio.h>
#include <stdlib.h>
#include "math.h"
#include "param.c"
/*#include "unix.h"*/

 
double FF1(double u,double up,double t)
{
 double f1;

 f1= up;
 // printf("  up %f %f \n",up,f1);
 return(f1);
}

double FF2(double u,double up,double t)
{
 double f2;
 extern double lbd,k,omega,A0;
 //printf(" l %f %f %f %f u %f up %f %f \n",lbd,k,omega,A0,u,up,t);
  f2=-lbd*up  -  k *  (u)  + A0 * sin(omega*t);
 // f2=-lbd*up  -  k *  sin(u)  +     omega *sin(u) *cos(u);
 return(f2);
}


/*
printf(" Resolution du pendule force                  \n");
printf("  2                                           \n");
printf(" d u        du                                \n");
printf(" --- + lbd  -- + k f(u)   = A0 sin (omega t)  \n");
printf("   2        dt                                \n");
printf(" dt                                           \n");
printf("                                              \n");
printf("   f(u)= u     lineaire !                                \n");*/