torsion en matlab

ANEXO 1
fprintf ('\t CALCULO EJE IDEAL SOMETIDO A TORSIÓN');
%Variables de entrada: m= material, l=luz, T= Torsor,
%a= sección transversal.
fprintf ('\n Materiales disponibles para el eje');fprintf ('\n 1. Acero.');
fprintf ('\n 2. Aluminio.');
fprintf ('\n 3. Titanio.');
fprintf ('\n 4. Bronce.');
fprintf ('\n 5. Laton');
m= input('\n ¿De que material requiere el eje? ---->');
l=input('\n ¿De cuánto es el valor de luz (m)? ---->');
T= input('\n ¿De cuánto es el valor del torsor (N.m)? ---->');
fprintf ('\n Sección trasversal disponibles para el eje');
fprintf ('\n 1.Circular maciza');
fprintf ('\n 2. Circular con hueco');
fprintf ('\n 3. Cuadrada');
fprintf ('\n 4. Cuadrada con hueco');
a= input('\n ¿Qué sección trasvesal requiere para el eje? ---->');

%ENM=Esfuerzo normal máximo, ECM= Esfuerzo cortante máximo,
%AN= Ángulo máximo de torsión.
ENM=0;ECM=0;AN=0;
%Modulos de elasticidad cortate (N/m2): st=acero estructural A36 , al=aleación aluminio forjado2014-T6
%ti= aleación de titanio Ti-6Al-4V, br=bronce C86100, la=latón rojo C83400
%r1= radio exterior de la sección trasversal
%r2= radio interior de la sección trasversal, t= espesor
%J=Momento polar de inercia
st= 75000000000;
al= 27000000000;
ti= 44000000000;
br= 38000000000;
la= 37000000000;
r1=0;J=0;r2=0;t=0;Ma=0;
%Ma= material en el swith
if m==1
Ma=st;
end
if m==2Ma=al;
end
if m==3
Ma=ti;
end
if m==4
Ma=br;
end
if m==5
Ma=la;
end

switch a
case 1
r1= input('\n ¿Qué radio requiere para la seccióntrasversal circular (m)? ---->');
J= (3.141592/2)*(r1^4);
ECM= ((T*r1)/J);
fprintf ('\n El esfuerzo cortante máximo para el eje es de ----> %2f', ECM);fprintf ('\t N/m2');
AN= (T*l)/(J*Ma);
AN= AN *(57.29577);
fprintf ('\n El ángulo máximo para el eje sometido a torsión es de ---->...