i want to solve these two coupled eqaut with finite boundary conditionsions. Can some one help me

eq1:=diff(f(eta),eta,eta,eta,eta)+2*(epsilon/(1+epsilon*theta(eta)))^2*diff(f(eta),eta,eta)*((diff(theta(eta),eta))^2)-(epsilon/(1+epsilon*theta(eta)))*(diff(theta(eta),eta,eta))*(diff(f(eta),eta,eta))=0;

eq2:=diff(theta(eta),eta,eta)+Pr*Re*f(eta)*diff(theta(eta),eta)=0;

Re:=1:Pr:=1:epsilon:=0.25:

bc:=f(1)=0,D(f)(-1)=0,D(f)(1)=1,D(f)(-1)=1,theta(-1)=0,theta(1)=0;