mpower (Input must be scalar and power matrix)

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

Answer 1

[Torsten]

Set a dot (.) before each ^2-operation.

Best wishes

Torsten.

[FarisM]

Thank you so much Torsten. Best regards Faris