flow over a cycle

1.

%FUNCTION

function [dragforce]=force(c_d,a,roh,v)
v=[1:50]
dragforce=roh*a*v.^2*c_d*0.5;
plot(v,dragforce,\'color\',\'g\')
xlabel(\'v\')
ylabel(\'dragforce\')
end

%USING FUNCTION TO CLACULATE DRAG FORCE

[dragforce]=force(0.1,0.1,0.1,50)

2.

%function

function [dragforce]=force(a,roh,v)
c_d=[0.47 0.2 0.27 0.5 0.1]
dragforce=roh*a*v.^2*c_d*0.5;
plot(c_d,dragforce,\'color\',\'g\')
xlabel(\'c_d\')
ylabel(\'dragforce\')
end

%using the function

[dragforce]=force(0.1,0.1,1)

)