function [ theta, E, D, FL, EXT, E2, D2, ER, IR] = FindAll(AngleMax,necktodiameter,alpha,beta) global a b a = 52; b = 10;
theta = (AngleMax)-(2*asind(1./(necktodiameter))); E = cosd(theta/2) – sind(alpha)*cosd(beta)*cosd(a)*cosd(b); D = ((1-((sind(alpha)).^2*(cosd(beta)^2).^2))*(1-(cosd(a)).^2*(cosd(b).^2))) – E.^2; FL = acosd(((((-sind(beta))*(cosd(a))*(sind(b)) + (cosd(alpha))*(cosd(beta))*(sind(a)))*E)… – (((cosd(alpha))*(cosd(beta))*(cosd(a))*(sind(b))+(sind(beta))*(sind(a)))*sqrt(D)))… /((1-sind(alpha)^2*cosd(beta)^2)*(1-cosd(a)^2*cosd(b)^2))); EXT = acosd(((((-sind(beta))*(cosd(a))*(sind(b)) + (cosd(alpha))*(cosd(beta))*(sind(a)))*E)… + (((cosd(alpha))*(cosd(beta))*(cosd(a))*(sind(b))+(sind(beta))*(sind(a)))*sqrt(D)))… /((1-sind(alpha)^2*cosd(beta)^2)*(1-cosd(a)^2*cosd(b)^2))); E2 = cosd(theta/2) – (cosd(alpha)*cosd(beta)*sind(a)); D2 = ((1 – cosd(alpha).^2*cosd(beta).^2)*cosd(a).^2) – E2.^2; ER = acosd (((-sind(beta)*sind(b) + sind(alpha)*cosd(beta)*cosd(b))*E2 … + (sind(alpha)*cosd(beta)*sind(b) + sind(beta)*cosd(b))*sqrt(D2))… /(cosd(a)*(1-(cosd(alpha)^2*cosd(beta)^2)))); IR = acosd (((-sind(beta)*sind(b) + sind(alpha)*cosd(beta)*cosd(b))*E2 … – (sind(alpha)*cosd(beta)*sind(b) + sind(beta)*cosd(b))*sqrt(D2))… /(cosd(a)*(1-(cosd(alpha)^2*cosd(beta)^2))));
Dear all, I try to use many values of 'beta', with others input remain constant, but error of mpower occurred. (input must be scalar and power matrix)
For example: >> [ theta, E, D, FL, EXT, E2, D2, ER, IR] = FindAll(180,2,35,[10,15,20]) ??? Error using ==> mpower Inputs must be a scalar and a square matrix. To compute elementwise POWER, use POWER (.^) instead.
Error in ==> FindAll at 24 D = ((1-((sind(alpha)).^2*(cosd(beta)^2).^2))*(1-(cosd(a)).^2*(cosd(b).^2))) – E.^2;
Please help me
Best Answer