plotting - Finding undetermined coefficients with boundary conditions

I have 2 equations which are;

To[x_] := c1 E^(-a1 x)/2 + c2 E^(a1 x)/2 + a2

Tu[x_] := c3 E^(-a3 x)/2 + c4 E^(a3 x)/2

Tu equals almost zero (10^-6) at end of the wires.

To must be equal to Tu at $±L/2$. Also To has to be continous at $±L/2$ by using Drichlet boundary condition which is To'[L/2] == a4 Tu'[L/2].

In order to understand easily, Tu shows ends of figures, and To shows middle of the figure.

All things considered, I have 4 coefficient and 4 boundary conditions which have to satisfy each other. Functions have to be symmetric , therefore I wrote a trigonometric form as a $Cosh[x]$.

How can I find $c1$ and $c2$ to satisfy all boundary conditions?

I wrote a program with Mathematica. But it is wrong. I think that something is wrong on finding $c1$. Can anyone help my code?

Here is my code:

To[x_] := 2 c1 Cosh[x a1] + a2
Tu[x_] := 2 c2 Cosh[x a3] 

ff = Solve[Tu[Le] == 0.000001,c2]

bcs = {To[L/2] == Tu[L/2], To'[L/2] == a4 Tu'[L/2]};
s = First@Solve[bcs[[1]], c1]
Quiet@Solve[bcs[[2]] /. s[[1]], c1]