Как задать условие для функции

double Psi(double x)

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
 
int Nt = 51; //Количество расчетных узлов (время)
int Nx = 25; //Количество расчетных узлов (координата)
double A = 5; //Коэффициент в уравнении
double ht = 0.25; //Шаг по времени
double hx = 0.25; //Шаг по координате
double X0 = -2; //Границы по координате
double XN = 4;
double X; //Положение точки
const int MaxNt = 200;
const int MaxNx = 200;
double U[MaxNt][MaxNx]; //Массив для хранения данных
double UPhi[MaxNt][MaxNx]; //Массив для хранения данных
double UPsi[MaxNt][MaxNx]; //Массив для хранения данных
 
double Phi(double x) //Положение точки
{
    return (2*x)/(1+pow(x,2)); //Начальное положение будет задано такой функцией
}
 
double Psi(double x) //Импульс точки
{
    if (x<0)  return 0;
    if (x>2)  return 0;
    if ((x<=0) && (x<=1))  return (x);
    if ((x<=1) && (x<=2))   return (2-x);
    //return cos(0.1*x);//Например, начальный импульс будет задан такой функцией
}
 
double Zero(double x) //Нулевая функция
{
    return 0.0;
}
 
void Print(char * FileName, char * Header, double U[MaxNt][MaxNx])
{
  FILE * F = fopen(FileName,"w");
  fprintf(F,"%10s \n%10s",Header,"T\t");
  //Вывести X
  for (int nx=0; nx<Nx; nx++) fprintf(F,"%10.3lf\t",X0+nx*hx);
  fprintf(F,"\n");
  //Вывести решение
  for (int nt=0; nt<Nt; nt++)
   {
    fprintf(F,"%10.3lf\t",nt*ht); //Время
    for (int nx=0; nx<Nx; nx++)
     fprintf(F,"%10.3lf\t",U[nt][nx]);
    fprintf(F,"\n");
   }
  fclose(F);
}
 
void Print() //Вывод результата на консоль и в файлы
{
    Print("resultPhi.txt","Phi",UPhi); //
    Print("resultPsi.txt","Psi",UPsi);
    Print("result.txt","Полное",U); //Полное
}
 
double Int(double x0,double xn,double F(double),int N)
//Интеграл методом прямоугольников
//x0 - начало интервала
//xn - конец интервала
//F - интегрируемая функция
//N - количество точек
{
    double h = (xn-x0)/(N-1); //Шаг
    double Sum = 0;
    for (int n=0; n<N; n++)
        Sum+=F(x0+h*n);
    return Sum / N;
}
 
double Integral(double x0,double xn,double F(double))
//Любым численным методом найти интеграл заданной функции
//Метод прямоугольников
{
    const double epsilon = 0.0001;//Точность
    int N = 20; //Начальное количество
    double I0 = Int(x0,xn,F,N);
    double I1;
    do {
        N*=2; //Увеличить количество шагов вдвое
        I1 = I0;
        I0 = Int(x0,xn,F,N); //Вычислить интеграл
    } while (fabs(I0-I1)>epsilon); //Если он очень похож на предыдущее значение
    return I0; //То это и есть искомое значение
}
 
double Solve(double Phi(double),double Psi(double),double A, double t,double x)
//Решение уравнения колебания струны
{
    //Решение:
    //U(x,t) = 0.5(Phi(x-at)+Phi(x+at)) + 1/(2*A)Integral(x-at,x+at,Psi);
    return 0.5*(Phi(x-A*t)+Phi(x+A*t)) + 0.5/A * Integral(x-A*t,x+A*t,Psi);
}
 
int main()
{
    system("chcp 1251>nul"); //Переключиться на кирилицу. Не забыть выбрать консольный шрифт!
 
    printf("Введите параметр A :");scanf("%lf",&A);
 //   printf("Сколько точек (по времени)вычислить? :"); scanf("%d",&Nt);
 //   printf("Сколько точек (по координате) вычислить? :"); scanf("%d",&Nx);
    if (Nt>MaxNt) {
        Nt = MaxNt;
        printf("Слишком много шагов по времени. Не больше %d\n",Nt);
    };
    if (Nx>MaxNx) {
        Nx = MaxNx;
        printf("Слишком много по координате. Не больше %d\n",Nx);
    };
 
    //printf("Шаг по времени? :"); scanf("%lf",&ht);
    //printf("Шаг по координате? :"); scanf("%lf",&hx);
    //printf("Координата точки x :");scanf("%lf",&X);
    //Вычисления
     //Для каждого момента времени от 0 с шагом 0.25 вычислить положение точки (полное решение)
     //а также частичные рещения
    printf("\n");
    for (int nx = 0; nx < Nx; nx++) //Для всех координат
    for (int nt = 0; nt < Nt; nt++) //Для всех времен
    {
        printf("\r%d",nt);
        double t = nt * ht; //Время
        X = X0 + nx * hx;
        U[nt][nx] = Solve(&Phi,&Psi,A,t,X); //Общее
        UPhi[nt][nx] = Solve(&Zero,&Psi,A,t,X); //Частное
        UPsi[nt][nx] = Solve(&Phi,&Zero,A,t,X); //Другое частное
    }
    printf("\r");
 
    //Вывод результата
    Print();
 
    system("pause");
}

#include <iostream>
#include <cstdlib>
using std::cin;
using std::cout;
using std::endl;
double myfunc (double x)
{
    double result;
    if (0<=x && x<=1)
    {
        cout << x << endl;
        return 0;
    }
    else if (1<x && x<=2)
    {
        cout << 2-x << endl;
        return 0;
    }
    else
        cout << 0 << endl;
        return 0;
}
int main()
{  
    double x;
    cout << "input X" << endl;
    cin >> x;
    myfunc(x);
    return 0;
}

double Psi(double x)
{
if( x< 0 ) return 0;
else if( x> 2 ) return 0;
else if( x<= 1 ) return x;
else return 2-x
}

double Psi(double x)
{
if (0<=x && x<=1)
        return x;
    else if (1<x && x<=2)
        return (2-x);
    else return 0;
}