Версия для печати темы

Автор: punkska 28.02.2007 3:16

Задание

x*tg(x)-1/3=0
x принадлежит [0.2;1]
x0=?

% Задача паскалевская. Но теперь её нееобходимо решить в пакете МатЛаб
% Тут тоже есть начальное значение (х0) Но нет эпсилант с каторой мы в паскале сравнивали значение ф-ции
% Но это всё лирика

Изучая приложение (его хелп), смог написать

x0=input('Введите начальное приближение')
x=linspace(x0,1,9)

% берём значения х от значения х0 ( я думаю 0.2) до 1 с 9 шагами (интервалами)

for x=x0:1
fun=x*tg(x)-1/3
fplot('fun',[x0,1])
end

% Вроде даже работает, но график прямая под углом 45!

Grid on
Xlable ('Ось Х')
Ylable ('')
Title('...')



Совсем не получается воспользоваться функцией Fzero для нахождения корней! Возможно я выбрал неудачно способ задания диапозона....
Может вообще всё неверно....
Опытные люди подскажи как правильно делать данные операция в матлабе, какие команды использовать