how can i solve the problem in line 22

close all
clear all
C=51e-9; % set the capacitance of the beam
A=4.905; % driving acceleration
f=(0:1:73)';
w=2*pi*f; % driving frequency
d=-190e-12; % piezoelectric strain coefficient
e=1.594e-8; % dielectric constant of the piezoelectric material
Y=66e9; % material Young's Modulus
k=0.29; % coupling coefficient
width=12.7e-3; % beam width
lb=21.5e-3;
lm=2e-3;
le=31.5e-3;
tc=0.19e-3;
tsh=0.13e-3;
b=0.16e-3;
R=60000;
s=0.02858; % damping ratio
b1=(3*b*(2*lb+lm-le))/(lb^2*(2*lb+(3*lm/2)));
P=1/(2*w^2)*(A^2*R*C^2*(Y*d*tc*b1/e)^2)/((4*s^2+k^4)*(R*C*w)^2+4*s*k^2*R*C*w+4*s^2);% Power equation

figure 1
plot(w,P)

Answer 1

[Torsten]

f is a vector. Thus you will have to replace all "*" by ".*", all "^" by ".^" and all "/" by "./" in the power equation.

Best wishes

Torsten.