%a
clear, clc, close all
fv=@(v) 1.5*sin(v)-sin(v).^3-0.64;
v=(0:pi/200:pi/2)';
y=fv(v);
plot (v,y)
grid
hold on

%v1,0=0.55 v2,0=1.05
%b
fpv=@(v) 1.5*cos(v)-3*sin(v).^2.*cos(v);
v0=[0.55 1.05];
for i =1:2
    v=v0(i);
    while abs(fv(v))>1e-6
        v=v-fv(v)/fpv(v);
    end
    plot(v,fv(v),'*')
    disp([v fv(v)])
end

%c
clear, figure

xv=@(v) cos(v).^3;
yv=@(v) 1.5*sin(v)-sin(v).^3;

v=0:pi/200:pi/2;
x=xv(v);
y=yv(v);
plot(x,y)