I want the value of 'gama1' to be obtained from the following equation

following is the code that i have written for obtaining the value of the 'gama1' variable. The issue i am facing is that when i run the code, i am gettin an error message that 'gama1' is an unknown parameter while i am actually trying to find out 'gama1' using the given FZERO function.

zeta=linspace(0,2);
    alpha=10;
    k=54;
    c=500;
    rho=7800;
    w=0.5;
    m=0.019;
    n=0.015;
    v=-0.38;
    V=15;
    t1=0.5;
    %Tint(1)=500;
    beta=1;
    Tw=30;
    mu0=3750e6;
    h=0.025e-3;
    gama_p=0.01;
y=500;
Q=zeros(10,1);
 i=[1:1:10];
Tint=zeros(10,1);
for i=1:1:10,
    Tint(i)=y;
    u=i;

     y=y+100;
      Q(i)=i;
  end
  %display(Q)

for i=1:10,
     if Tint(i)<=955,
             mu(i)=1-(3.44e-4)*Tint(i);
     else 
             A=1-((Tint(i)-955)./(1500-955));
             mu(i)=0.68*power(A,1.7);
         end;
end;

      lamda=atan(mu)*180/pi;
      phi = 35 + (0.5*(alpha-(atan(mu)*180/pi)));
      vc=V*sind(phi)/cosd(phi-alpha);
      %lc=t1*sind(phi+lamda-alpha)./(sind(phi)*cosd(lamda));
      lc=t1*(sind(phi-alpha)+(tand(lamda).*cosd(phi-alpha)))/sind(phi);
      %strain=tand(phi-alpha)+(1./tand(phi));

      ref_strain=(V*cosd(alpha))./(h*cosd(phi-alpha));
      ref_stress=mu0*(Tw^v).*power(ref_strain,m);
      B=(beta*ref_stress)/(rho*c*Tw);
      C=(sind(phi).*cosd(phi-alpha))./cosd(alpha);
      D=((rho*(V^2))*sind(phi)*cosd(alpha))./(ref_stress.*cosd(phi-alpha));
      K=k*cosd(phi-alpha)./(rho*c*h*V*cosd(alpha));
      gmma2=(1./tand(phi))+tand(phi-alpha);

      T0= (4.5e-4-(8*Tw*K))./((4.5e-4)-20*K);
      T1=(4*T0-Tw)/3;
      %gama1=0;

      der1=(gama1-gama_p)/(2*h);
  %     AA=m*C.*power(T0,v)*(gama_p^n)*power(der1,(m-1))* (gama1+gama_p)/(h^2);
  %     BB= v*C.*power(T0,v-1)*(gama_p^(n))*power(der1,m).*((T1-Tw)/(2*h));
  %     ZZ=n*C*power(der1,m+1).*power(T0,v)*(gama_p^(n-1))-C.*D*der1;

%     function F=hola(gama1,gama_p,h,m,C,T0,v,n,T1,Tw,D)
    %F=@(gama1) AA+BB+ZZ;
    %syms gama1

      %F=inline('m*C.*power(T0,v)*(gama_p^n)*power(der1,(m-1))* (gama1+gama_p)/(h^2))+(v*C.*power(T0,v-1)*(gama_p^(n))*power(der1,m).*((T1-Tw)/(2*h)))+(n*C*power(der1,m+1).*power(T0,v)*(gama_p^(n-1))-C.*D*der1')
      F=@(gama1) (m*C.*power(T0,v)*(gama_p^n)*power(der1,(m-1))* (gama1+gama_p)/(h^2))+(v*C.*power(T0,v-1)*(gama_p^(n))*power(der1,m).*((T1-Tw)/(2*h)))+(n*C*power(der1,m+1).*power(T0,v)*(gama_p^(n-1))-C.*D*der1);

      gama1=fzero(F,0);
      M=(power(T0,v));
      N=M';
      tou0=C*(((gama1-gama_p)/(2*h))^(m))*N*(gama_p^(n));
      tou=rho*(power(V*sin(phi),2)*gama1)+tou0;

Answer 1

[Purushottama Rao]

gamma1 variable has been used in the equations prior to its declaration. MATLAB does not understand what is gamma1 until it faces gamma1=fzero(F,0)