I want to solve two odes with their boundary condition. I wrote the code below:
bcs := phi(0)=0,phi(h)=1,T(0)=0,T(h)=1;
however, this code doesnt get my desired results (the results are complex!). but when I (with hand) integrate Eq3 twice and substitute boundary conditions and replace in Eq2 the answer is easy and straightforward.
How can I change the following algorithm to get my results?
Thanks for your attention in advance