Please help me to write character in vertical axis of the run figre as in attached photo
clc;clf;clear;
myLegend1 = {};
rr = [ 0.4 0.5 0.6]
for i =1:numel(rr)
x=0:0.0001:9;
alfa= rr(i);
s=1;h=1;
L1=-0.5;L2=0.4;
c=30;
y=0.5;b=1;
k0=386; ce=3.831*10^2; mu=38.6*10^9;alfat=1.78*10^-5; rho=89.54*10^2; lamda=77.6*10^9;taw=0.5;Tnot=2.93*10^2;
c0=sqrt((lamda+2*mu)/(rho)); Betanot=(3*lamda+2*mu)*alfat; a1=mu/(lamda+2*mu);a2=(mu+lamda)/(lamda+2*mu);a3=(Betanot*Tnot)/(lamda+2*mu);omega=(rho*ce)/(k0);
a4=lamda/(lamda+2*mu);a5=(k0*omega*c0^2)/(k0);a6=(rho*ce*c0^2)/(k0);
a7=(Betanot*c0^2)/(k0); a8=a6*taw; a9=a7*taw; a10=rho*ce*taw*omega*c0^4/(k0); a11=Betanot*taw*omega*c0^4/(k0);w=rho*ce/(k0);
b1=s^2+a1*h^2-(x+y+b).^2;b2=a2*s*h;b3=0;b4=a3*s;b5=alfa*b1;b6=alfa*b2;b7=alfa*b3;b8=2*alfa*b4;
B1=b1*L1^2+b2*L2^2;B2=-b3*L1-b4;
B3=-b5*L1^2-b6*L2^2;
%B4=b7*L1+b8;
B4=B3./4;
sec= ((B1-B2)./(3*B4))-((B1-B2)./(3*B4))+c;
T=-((B1-B2)./(3*B4))+(sec).*exp(-(s*x+h*y+b));
% fprintf(‘%12.10f %15.8f rn’,-((B1-B2)./(3*B4)),(c+((B1-B2)./(3*B4))));
%fir=c+sec;
% N2=zeros(1,1);
% N2=(s*L1+a4*h*L2).*-1*(c+((B1-B2)./(3*B4))).*exp(-(s*x+h*y+b)).*(1-alfa.*T)+a3.*T.*(alfa.*T-1);
% N3=zeros(1,1);
% N3=a1*(h*L1+s*L2)*(1-alfa.*T).*(-1*(c+((B1-B2)./(3*B4))).*exp(-(s*x+h*y+b)));
figure(1)
plot(x,T)
grid on,hold on
myLegend1{i}=[‘alfa= ‘,num2str(rr(i))];
i=i+1;
end
figure(1)
legend(myLegend1)
hold onclc;clf;clear;
myLegend1 = {};
rr = [ 0.4 0.5 0.6]
for i =1:numel(rr)
x=0:0.0001:9;
alfa= rr(i);
s=1;h=1;
L1=-0.5;L2=0.4;
c=30;
y=0.5;b=1;
k0=386; ce=3.831*10^2; mu=38.6*10^9;alfat=1.78*10^-5; rho=89.54*10^2; lamda=77.6*10^9;taw=0.5;Tnot=2.93*10^2;
c0=sqrt((lamda+2*mu)/(rho)); Betanot=(3*lamda+2*mu)*alfat; a1=mu/(lamda+2*mu);a2=(mu+lamda)/(lamda+2*mu);a3=(Betanot*Tnot)/(lamda+2*mu);omega=(rho*ce)/(k0);
a4=lamda/(lamda+2*mu);a5=(k0*omega*c0^2)/(k0);a6=(rho*ce*c0^2)/(k0);
a7=(Betanot*c0^2)/(k0); a8=a6*taw; a9=a7*taw; a10=rho*ce*taw*omega*c0^4/(k0); a11=Betanot*taw*omega*c0^4/(k0);w=rho*ce/(k0);
b1=s^2+a1*h^2-(x+y+b).^2;b2=a2*s*h;b3=0;b4=a3*s;b5=alfa*b1;b6=alfa*b2;b7=alfa*b3;b8=2*alfa*b4;
B1=b1*L1^2+b2*L2^2;B2=-b3*L1-b4;
B3=-b5*L1^2-b6*L2^2;
%B4=b7*L1+b8;
B4=B3./4;
sec= ((B1-B2)./(3*B4))-((B1-B2)./(3*B4))+c;
T=-((B1-B2)./(3*B4))+(sec).*exp(-(s*x+h*y+b));
% fprintf(‘%12.10f %15.8f rn’,-((B1-B2)./(3*B4)),(c+((B1-B2)./(3*B4))));
%fir=c+sec;
% N2=zeros(1,1);
% N2=(s*L1+a4*h*L2).*-1*(c+((B1-B2)./(3*B4))).*exp(-(s*x+h*y+b)).*(1-alfa.*T)+a3.*T.*(alfa.*T-1);
% N3=zeros(1,1);
% N3=a1*(h*L1+s*L2)*(1-alfa.*T).*(-1*(c+((B1-B2)./(3*B4))).*exp(-(s*x+h*y+b)));
figure(1)
plot(x,T)
grid on,hold on
myLegend1{i}=[‘alfa= ‘,num2str(rr(i))];
i=i+1;
end
figure(1)
legend(myLegend1)
hold on clc;clf;clear;
myLegend1 = {};
rr = [ 0.4 0.5 0.6]
for i =1:numel(rr)
x=0:0.0001:9;
alfa= rr(i);
s=1;h=1;
L1=-0.5;L2=0.4;
c=30;
y=0.5;b=1;
k0=386; ce=3.831*10^2; mu=38.6*10^9;alfat=1.78*10^-5; rho=89.54*10^2; lamda=77.6*10^9;taw=0.5;Tnot=2.93*10^2;
c0=sqrt((lamda+2*mu)/(rho)); Betanot=(3*lamda+2*mu)*alfat; a1=mu/(lamda+2*mu);a2=(mu+lamda)/(lamda+2*mu);a3=(Betanot*Tnot)/(lamda+2*mu);omega=(rho*ce)/(k0);
a4=lamda/(lamda+2*mu);a5=(k0*omega*c0^2)/(k0);a6=(rho*ce*c0^2)/(k0);
a7=(Betanot*c0^2)/(k0); a8=a6*taw; a9=a7*taw; a10=rho*ce*taw*omega*c0^4/(k0); a11=Betanot*taw*omega*c0^4/(k0);w=rho*ce/(k0);
b1=s^2+a1*h^2-(x+y+b).^2;b2=a2*s*h;b3=0;b4=a3*s;b5=alfa*b1;b6=alfa*b2;b7=alfa*b3;b8=2*alfa*b4;
B1=b1*L1^2+b2*L2^2;B2=-b3*L1-b4;
B3=-b5*L1^2-b6*L2^2;
%B4=b7*L1+b8;
B4=B3./4;
sec= ((B1-B2)./(3*B4))-((B1-B2)./(3*B4))+c;
T=-((B1-B2)./(3*B4))+(sec).*exp(-(s*x+h*y+b));
% fprintf(‘%12.10f %15.8f rn’,-((B1-B2)./(3*B4)),(c+((B1-B2)./(3*B4))));
%fir=c+sec;
% N2=zeros(1,1);
% N2=(s*L1+a4*h*L2).*-1*(c+((B1-B2)./(3*B4))).*exp(-(s*x+h*y+b)).*(1-alfa.*T)+a3.*T.*(alfa.*T-1);
% N3=zeros(1,1);
% N3=a1*(h*L1+s*L2)*(1-alfa.*T).*(-1*(c+((B1-B2)./(3*B4))).*exp(-(s*x+h*y+b)));
figure(1)
plot(x,T)
grid on,hold on
myLegend1{i}=[‘alfa= ‘,num2str(rr(i))];
i=i+1;
end
figure(1)
legend(myLegend1)
hold on figure MATLAB Answers — New Questions









