Форум программистов
 

Восстановите пароль или Зарегистрируйтесь на форуме, о проблемах и с заказом рекламы пишите сюда - alarforum@yandex.ru, проверяйте папку спам!

Вернуться   Форум программистов > IT форум > Помощь студентам
Регистрация

Восстановить пароль
Повторная активизация e-mail

Купить рекламу на форуме - 42 тыс руб за месяц

Ответ
 
Опции темы Поиск в этой теме
Старый 03.02.2016, 17:38   #1
mital25
Пользователь
 
Регистрация: 06.03.2014
Сообщений: 28
По умолчанию по Scilab

доброго времени суток Господа! решаю тут эту хреновину и запарился в конец тангенс не рисует гладким а углами что я сделал не так? и куда то убегает синяя линия графика раз через раз .вот мой код
Код:
clear all
 
global e a m h x
a=1e-9;
m=9.1e-31;
h=1.05e-34;
e=1.6e-19;
 
step = 0.001;
V=5;

function le=aa(E)
    le=tan(a*sqrt(2.*m*E*e)/h) ;
endfunction

function ri=bb(V,E)
    G=V-E;
    G1=sqrt((2.*m*E)/(h^2))//G1=E.*G;
    G2=sqrt((2.*m*(V-E))/(h^2))//G2=(E.*2-V) ;
    ri=2*sqrt(G1)./G2;
endfunction

function y = fun(E)
    y =(1/sqrt((2.*m*E)/(h^2)))*tan((sqrt((2.*m*E)/(h^2)))*a)+(1/(sqrt((2.*m*(V-E))/(h^2)))) //tan(a*sqrt(2.*m*E*e)/h) - 2*sqrt(E.*(V-E))./(E.*2-V);
endfunction

Emin = 0;
Emax = 1;

clf()

E = [Emin*1.01: step: Emax];//(%pi*1.5*h/a)^2/(2*m*e)*0.9];
y1 = aa(E);
y2 = bb(E, V);

plot( E, [y1], 'r')
plot( E, [y2], 'b')
mtlb_axis([min(E) max(E) 0 V])
xgrid

x = fsolve(, fun);
disp(sprintf("E=%f, fun(E)=%f", x, fun(x)))

x = fsolve(, fun);
disp(sprintf("E=%f, fun(E)=%f", x, fun(x)))

x = fsolve(, fun);
disp(sprintf("E=%f, fun(E)=%f", x, fun(x)))

x = fsolve(, fun);
disp(sprintf("E=%f, fun(E)=%f", x, fun(x)))
mital25 вне форума Ответить с цитированием
Ответ


Купить рекламу на форуме - 42 тыс руб за месяц

Опции темы Поиск в этой теме
Поиск в этой теме:

Расширенный поиск


Похожие темы
Тема Автор Раздел Ответов Последнее сообщение
Scilab. While. Егор_Чита Помощь студентам 2 10.11.2015 09:48
задача по Scilab Максим2912 Помощь студентам 1 28.06.2012 12:37
Matlab/Scilab sawyer1204 Помощь студентам 1 24.06.2012 12:28
Scilab задачи. bubu Помощь студентам 1 27.12.2011 20:44
Scilab Igel93 Помощь студентам 0 20.12.2011 22:17